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SECTION  1 
INTRODUCTION 


The  WESCOM  code  is  a  digital  computer  program  for  use  in  evaluating 
nuclear  weapon  effects  on  satellite  communications  systems. t  The  code  is 
intended  for  use  in  the  evaluation  of  system  performance  for  a  variety  of 
weapon  and  system  parameters.  The  code  is  designed  for  operation  on  most 
computers  used  for  system  studies.  Operational  information  is  given  in  Volume  1 
(User's  Manual).  In  this  volume  a  description  of  the  code  structure  and 
selected  routines  is  given  for  use  in  conjunction  with  FORTRAN  source  listings. 
Volume  3  provides  summary  descriptions  of  the  environment,  propagation,  and 
signal  processing  models  for  users  and  analysts  requiring  an  understanding  of 
the  level  of  modeling  used  in  the  code. 

The  WESCOM  code  consists  of  a  main  driver  (control)  routine  and  a 
number  of  subroutines  and  functions.  The  FORTRAN  names  and  a  brief  description 
of  the  purpose  of  each  routine  are  given  in  Table  1.  Section  2  of  this  volume 
describes  the  interconnections  between  routines  and  common  blocks.  Intermediate 
level  descriptions  of  selected  routines  are  given  in  Section  3. 


*  WESCOM  is  an  acronym  for  Weapon  Effects  on  Satellite  COMmunications. 
t  Excluding  effects  on  system  materials. 
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Table  1  Description  of  WESCOM  routines. 


Name  of 
Control 
ACCUM 

AM3ABS 

ANLYT2 

WIND 

ARRLIM 

ATMOS F 

ATMOS I 

ATMOS U 

AZEER 

BALFIT 

BB 

BDPSK 

BETAS 


Routine  Description 

Main  driver  routine  for  WESCOM  code 

Integrates  f(x)dx  for  a  specified  interval  in 
x  assuming  either  a  linear  or  exponential 
variation 

Computes  the  absorption  from  ambient  oxygen 
(0^)  and  water  vap”  along  a  specified  path 

Computes  solution  to  differential  equation 
used  in  deionization  above  100  km 

Driver  routine  for  antenna  noise  module 

Determines  debris  location  and  size  from 
debris  marker  particle  location 

Provides  properties  of  the  normal  or  heaved 
atmosphere  above  100  km 

Initializes  atmospheric  properties  for  routine 
ATMOS U 

Determines  properties  of  the  undisturbed 
atmosphere 

Computes  true  azimuth,  elevation,  and  slant 
range  of  one  point  relative  to  another  point 
(points  specified  by  vectors) 

Fits  a  fourth  order  curve  to  ballistic  equation 
solution 

Calculate  maximum  FM  baseband  frequency 

Computes  probability  bit  error  for  a  binary 
differential  psk  system  in  a  fading  environ¬ 
ment 

Computes  absorption  and  backscattcr  for  beta 
sheath 


6 


Table  1 

Name  of 
BP  1 1  Ll) 

BIZ 

BKRIAD 

BLKDTA 

BTRACii 

CARDIN 

CHIiMD 

CHIiMDQ 

CIH-MII- 

Cl  1 1 -MIR 

Clll-MQ 

CUM I  ON 

CHXSPC 

Cl. OS  I 

CMDI.DT 

CMI-'liDT 

COMB 


i 


(Continued) 

Routine  Description 

Computes  properties  of  magnetic  dipole  field 
at  a  point 

Computes  function  used  in  Rooker  interpolation 

Determines  four  point  interpolation  using 
Booker  exponential  formulation 

Defines  parameters  and  constants  for  lablcd 
common  blocks 

Determines  location  of  ionized  air  parcels 
along  a  geomagnetic  field  line 

Reads  input  cards  and  writes  card  imago  output 

Driver  routine  for  D-region  chemistry  module 

Driver  routine  for  steady-state  D-region 
chemistry 

Determines  chemistry  and  ionization  for  E- 
and  l:- regions 

Determines  chemistry  and  ionization  for 
heated  regions 

Determines  steady-state  ionization  above 
100  km 

Determines  species  concentrations  in  I!-  and 
P- regions 


Determines  energy  histogram  for  charge  ex¬ 
change  particles 

Determines  location  of  closest  point  of  ap¬ 
proach  between  two  vectors 

lid  its  burst  list  for  bursts  to  be  used  in 
D-region  chemist ty  calculations 

lid  its  burst  list  for  bursts  to  be  used  in 
Pi-  and  ('-region  chemistry  calculations 

Determines  number  of  ways  of  choosing  I  quan- 
tities  from  a  total  of  M  quantities  without 
regard  to  order 
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Table  1  (Continued) 
Name  of  Routine 


Description 


CONS PL 

COORDV 

COORDX 


CORDTF 

CPSK 

CROSS 

DATE 

DBM  IT 

DBMET 

DBSTAM 

DEBRIS 

DECODE 

DHDEP 

DEEPSK 

1)11  I  NT 


Determines  energy  histogram  for  loss-cone 
part ic les 

Determines  vector  coordinates  from  latitude, 
longitude,  and  altitude 

Transforms  coordinates  of  a  point  in  one 
coordinate  system  to  any  of  the  other  three 
systems;  the  coordinate  systems  are: 

1.  ground  range,  azimuth,  altitude 
2  X  feast),  Y  (north),  Z  (vertical) 

3.  slant  range,  azimuth,  elevation 

4.  elevation,  azimuth,  altitude 

Determines  coordinate  conversions  for  input 
geometry  and  prepares  summary  output 

Computes  probability  of  bit  error  for  an 
M-ary  coherent  psk  system  in  a  noise  environ¬ 
ment  (no  fading) 

Calculates  cross  product  for  two  vectors 

Computes  Gregorian  date  and  zone  time  after 
specified  time  interval 

Determines  time-dependent  debris  region 
quantities  for  high-altitude  bursts 

Determines  t  i mo -dependent  debris  region  quail 
tities  for  low-altitude  bursts 

Driver  routine  for  post -stabi 1 i zat ion  debris 
geomet  ry 

Driver  routine  to  determine  energy  deposition 
by  heavy  particles 

Decodes  four  digit  integer  into  elements 

Determines  delayed  energy  deposition  rate 

Computes  probability  of  bit  error  for  a 
coherent  psk  system  for  slow  fading  environ 
merit 

Evaluates  three  intigrals  associated  with 
diffusion  ionization  profile 
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Table  1 
Name  of 
dish 

DOT 

DPSK 

DQI’SK 

DRATI- 

DSTINT 

DTNP.P 

DTNHQ 

DllSTMI 

P.ANOM 

F.DI.PB 

I'.DRPG 

PDF.PND 

nr 

HLDHN 


(Continued) 

Ro  u  t ine  Des  cription 

Calculates  exponential  term  in  debris  distri¬ 
bution  function 

Calculates  vector  dot  product 

Computes  probability  of  bit  error  for  an 
M-ary  differentially  coherent  psk  system  in 
a  noise  environment  (no  fading) 

Computes  probability  of  bit  error  for  a  dif¬ 
ferentially  coherent  4-ary  psk  system  for 
slow  fading 

Determines  reaction  rate  coefficients  for 
D-region  chemistry  solutions 

Determines  absorption  due  to  dust  clouds 

Computes  l)-rcgion  positive  ion  and  elec¬ 
tron  concentrations  due  to  prompt 
rad iat ion 

Computes  steady-state  D-rcgion  positive 
ion  and  electron  concentrations  due  to 
delayed  radiation 

Determines  time-independent  quantities 
for  dust  regions 

Calculates  eccentric  anomaly  for  satellite 
orbit 

Determines  energy  deposition  rate  due  to 
beta  particles 

Determines  energy  deposition  rate  due  to 
gamma  rays 

Determines  energy  deposition  rate  due  to 
neutron  elastic  collisions  and  capture 
reactions 

Calculates  the  exponential  integral  function 
for  positive  arguments 

Driver  routine  for  chemistry  and  ionization 
calculations 
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Name  of  Routine 


Description 


I  Nil  Driver  routine  lor  i-  .mil  I -region  chemistry 

moilu !  e 

I.NJVIB  Computes  the  concentration  of  \.  in  tin  t'l  rst 

v i hrat iona  1  st  ate 


I  QCI10S  locates  the  air  parcel  in  the  fireball  at 

apogee  that  will  stop  at  the  magnetic  equator 

i  \!  R  Determines  equilibrium  air  species  concent  ra- 

t  i  ons 


’  <11  nn 


Determines  equilibrium  aluminum  and  uranium 
spec i es  concent  rat  ions 


rur 


evaluates  error  function 


i  im: 

I  IIXI'I  i 


i:\ivit 


i:\rivr 


r: 


r\c 


rm:.\ 


I  R I  NT 


I  IUTMI1 


ruMin 


!  BMI  IT 


I  B'll  f 


Ivaluates  complement  of  the  error  function 

Determines  total  l)\  radiation  and  fraction 
radiated  in  five  radiation  groups 

Provides  exponential  interpolation 

Determines  the  integral  of  a  quantity  when 
the  quantity  is  defined  by  a  series  of  exponen¬ 
tials 

Ivaluates  function  used  in  beta  particle 
energy  deposition  routine 

Computes  k  factorial 

Determines  ionization  and  collision  frequen¬ 
cies  at  point  within  fireball 

Driver  routine  for  fireball  environment  and 
propagation  effects  module 

Computes  mass  density  and  mass  density  scale 
height  at  bottom  of  high-altitude  fireballs 

Determines  t ime- independent  fireball  quantities 
for  high-altitude  bursts 

Determines  time-dependent  fireball  quantities 
for  high-altitude  bursts 

Determines  t ime- independent  fireball  quantities 
for  low-altitude  bursts 
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Table  1  (Continued) 


Name  of 

FBMLT 

FBI  1ST 

FNNO 

FOOTPC 
FPC I R 

(TINT 

FPLNK 

FPPLOT 

FQAMB 

FSK'r.'B 

FSKWBQ 

FSKKBY 

FZRT 

GAIN 

GF.OCOR 

GNNO 

GOX 

GWAT 


Routine 


Description 


Determines  time-dependent  fireball  quantities 
for  low-alt itutde  bursts 

Determines  fireball  thermal  emission  in  a 
given  direction 

Determines  solution  of  differential  equation 
used  in  high-altitude  chemistry  model 

Driver  routine  for  fireball  footprint  plots 

Determines  circuit  performance  for  single 
link  circuit  used  for  footprint  plots 

Determines  parametric  values  of  footprint  plot 
quantities  as  a  function  of  ground  terminal 
locat ion 

Determines  link  propagation  quantities  for 
footprint  plots 

Prepares  footprint  printer  page  plots 

Determines  ion-pair  production  rate  below 
100  km  from  electron  density  and  atmospheric 
parameters 

Computer  probability  of  bit  error  for  M-ary 
wideband  FSK  system  in  a  noise  environment 
(no  fading) 

Computes  the  complement  of  the  normal  proba¬ 
bility  distribution 

Computes  a  summation  used  in  routine  FSKWB 
Solves  heavy  particle  range-energy  equation 
Calculates  antenna  gain 

Computes  latitude  and  longitude  of  position 
vector 

Determines  solution  to  set  of  equations  des¬ 
cribing  species  concentrations  for  high- 
altitude  chemistry  model 

Computes  the  0^  incremental  absorption 

Computes  the  incremental  absorption  from  water 
vapor 
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Name  of 
GWPLME 

HDINT 

HEAVE 

HELP 

HERCAN 

HPCHEM 

HTOS 

IDEFLT 

INITAL 

INPUT 

INPUTP 

IONLEK 

IONOSI 

IONOSU 

JAMND 

JULIAN 

LAG  RAN 

LEKSPC 


Routine  Description 

Determines  high-altitude  fireball  motion 
caused  by  gravity  and  atmospheric  winds 

Computes  intergral  used  in  diffusion  ioniza¬ 
tion  model 

Computes  the  vertical  velocity  and  time 
history  of  the  altitude,  expansion  ratio, 
and  density  scale  height  of  a  point  in  a 
vertically  heaving  atmosphere 

Prints  an  error  message  and  can  continue  or 
terminate  the  problem 

Computes  coefficients  used  by  routine  SANCAN 

Partitions  energy  lost  by  heavy  particles 
and  determines  change  in  neutral  and  ion 
species 

Calculates  slant  ranges  along  vector  to 
points  t*at  are  at  specified  altitude 

Prepares  input  default  data 

Determines  D-region  neutral  species  concen¬ 
trations  at  end  of  Phase  2  following  a  burst 

Reads  and  checks  input  data 

Performs  preliminary  processing  of  input 
data  and  prepares  summary  output 

Determines  energy  in  ion  leak  particles  and 
spatial  distribution  paramenter 

Initializes  ionospheric  properties  for 
routine  IONOSU 

Determines  properties  of  the  undisturbed 
ionosphere 

Determines  jammer  noise  density  at  specified 
receiver  for  a  specified  jammer  location 

Computes  Gregorian  calendar  date  to  a  Julian 
day  number 

Computes  the  initial  location  and  final 
properties  of  a  Lagrangian  cell 

Determines  energy  histogram  for  ion  leak 
particles 
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Table  1 
Name  of 

LINKER 

LOCATE 

LOCLAX 

LOSCON 

MAGI)  IS 

MAGE  IT 

MAGFLD 

MANPTS 

MATMUL 

MI  SCAT 

MODEM! 

MODEM2 

MODEM3 

NFSKMF 

NINT 

NORMAL 

OMEGA 


(Continued) 

Routine  Description 

Used  on  the  Honeywell  machine  for  linking 
(overlay) 

Determines  array  indices  and  weighting  func¬ 
tions  for  use  in  link  performance  models 

Calculates  vector  transformation  matrix 

Determines  energy  in  loss  cone  particles  and 
spatial  distribution  parameter 

Determines  altitude  and  length  along  geo¬ 
magnetic  field  line  of  a  point  on  the  field 
1  ine 

Fits  a  dipole  magnetic  field  to  the  local 
magnetic  field  at  specified  point 

Determines  location  of  the  point  on  a  magnetic 
field  line  that  is  at  a  specified  altitude 

Determines  mandatory  integration  points  re¬ 
lated  to  beta  tubes,  low-altitude  ionization 
sources,  and  heavy  particle  ionization  sources 

Performs  matrix  multiplication  between  two 
real  matrices 

Computes  the  Mie  backscatter  and  extinction 
efficiencies  given  the  normalized  particle 
diameter 


Determines  probability  of  message  error  rate 
for  fading  environment  for  three  modems 


Computes  probability  of  bit  error  rate  for 
non-cohercnt  FSK  for  fading  environment 

Determines  a  normalized  integral  over  the 
antenna  pattern  used  in  computing  receiver 
antenna  noise  temperature 

Calculates  the  probability  that  y  lies  between 
x  and  -x  for  a  standard  normal  distribution 

Computes  solid  angle  subtended  from  any  aspect 
for  spheroid,  skewed  spheroid,  torus,  cone, 
and  skewed  cone 
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Table  1  (Continued) 
Name  of  Routine 


Description 

Calculates  magnetic  field  at  specified  point 


0N1MG5 

ORBIT  Determines  satellite  orbit  parameters 

OZONE-  Determines  ozone  concent  rat  ior  s  below  55  km 


PCIIF.M 

PC LOUD 

n  m:r 

PGRODP 

PIlIiAV  n 

PI  I  P  NOM 

PIIOTOD 

PI  101  OR 

P I  ON  !•' 

PL  I  NT P 
PLYVAL 

pmass 

PROPFN 


Determines  change  in  species  densities  at  alti¬ 
tudes  above  100  km  due  to  prompt  radiation 

Computes  the  location,  shape,  and  velocity  of 
a  dust  particle  size  group 

Computes  the  energy  deposited  due  to  prompt 
rad i at  ion 

Computes  the  average  backscatter  and  extinc¬ 
tion  cross  sections  for  a  specified  dust 
particle  group 

Determines  the  heave  disturbed  atmosphere 
prof i lc 

Driver  routine  for  fireball  and  debris  phenom¬ 
enology  module 

Computes  the  fraction  of  0,  remaining  after 
thermal  radiation 

Calculates  the  negative  ion  photodetachment 
and  photodissociation  rates  due  to  thermal 
rad iat ion 

Obtains  the  initial  ion  concentrations  at  a 
point  above  100  km 

Obtains  the  power  law  interpolation,  y  =  ax^ 

Evaluates  a  polynomial  function  given  the 
coefficients  and  argument 

Computes  the  mass  penetrated  between  two  points 
assuming  that  the  p>r  density  varies  exponen¬ 
tially  with  alt i tude 

Driver  routine  for  calculating  propagation 
effects  along  a  given  path 
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Table  1  (Continued) 


Name  of 

PSPM 

PTH  I  NT 

QCOSM 

QFNTN 

radout 

RATE 

RICATT 

R INTER 

RRMSTP 

SANCAN 

SATC 

SATCIR 

SATLNK 

SEPA 

SUNK 

SLOBE 

son, 


Routine  Description 

Priver  routine  for  posr-stabil icat ion  debris 
module  for  specified  debris  region 

Computes  the  integral  of  the  electron  density, 
absorption,  refraction,  scintillation  effects, 
and  Faraday  rotation  along  a  propagation  path 

Computes  cosmic  ion-pair  production  rate 

Computes  Marcum's  Q  function 

Computes  the  radiated  power  for  a  burst  of  a 
given  yield  at  a  specified  altitude 

Calculates  the  reaction  rate  coefficient  for 
a  specified  reaction 

Computes  the  transient  electron  density  based 
on  an  approximate  solution  to  the  Ricatti 
cquat ion 

Computes  intersections  between  ray  path  and 
right  circular  and  skewed  spheroids,  cones, 
magnetically  defined  tubes,  and  torus 

Calculates  ratio  of  rms  to  peak  frequency 
deviation  of  an  FM  system 

Computes  two  integrals  used  in  analytic  signal 
processing  models 

Priver  routine  for  satellite  communication 
calculations 

lietcrmincs  performance  quantities  for  satellite 
communication  circuit 

Determines  link  propagation  quantities 

Calculates  the  angle  between  two  vectors 

Determines  whether  propagation  or  noise  calcu¬ 
lations  have  been  made  for  specified  link 

Determines  maximum  side  lobe  gain 

Provides  the  mass  density  and  dielectric 
constant  for  a  given  soil 
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Table  1 
Name  of 

SOLCYC 

SOLORB 

SOLZEN 

SORB  IT 

SPCTRM 

SPECDP 

SI’ECDQ 

START 

STR 1 1 

STRIP 

sdr\  i.c 
TATSRT 
TEMPZ.I 

TEST 

ri  xk 

THING 

TETEI.D 


{Continued) 

Routine  Description 


Computes  the  10. 7 -cm  solar  flux 

Computes  the  latitude  and  longitude  of  the 
subsolar  point 

Calculates  the  cosine  of  the  solar  zenith 
angle 

Determines  satellite  location  and  velocity  at 
specified  time 

Defines  device  data 

Computes  the  neutral  species  concentrations 
at  the  end  of  Phase  3 

Computes  the  modification  of  the  neutral 
species  due  to  delayed  radiation 

Computes  a  list  of  mandatory  integration  point 
ranked  in  increasing  distance  from  the  trans¬ 
mitter 

Determines  scintillation  quantities  for 
specified  path  through  fireball 

Defines  fireball  strintion  parameters 

Calculates  the  difference  between  two  vectors 

Computes  true  anomaly  of  satellite  orbit 

Determines  atmospheric  temperature  below 
120  km 

Test  input  quantities  to  see  if  they  are 
within  given  bounds 

Computes  the  excitation  temperature  for  a 
given  energy  density  equilibrated  between  t  he 
\ 2 ( V I B ) ,  0('b),  the  electron  kinetic  energy 
and  N(-!>) 

Computes  fireball  temperature  function  for 
high-altitude  bursts 

Determines  location  of  the  point  on  a  tilted 
axis  that  is  at  a  specified  altitude 
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Table  1  (Continued) 


Name  of 
TORBIT 

TOROID 

TUBE 

UN  I  TV 

USERMD 

UVWAV 

UVWD 

VECLIN 

VECM 

VOLUME 

WATER 

WOBD 

WOGD 

WOGP 

WOG1 

WOND 


Routine  Description 

Determines  ground  terminal  location  and 
velocity  at  specified  time 

Computes  the  intersection  points  of  a  ray  path 
with  a  toroidal  region 

Computes  the  intersection  point  of  a  ray  path 
with  a  magnetically  contained  region 

Computes  a  unit  vector  from  a  given  input 
vector 

User  routine  to  compute  circuit  performance 

Determines  average  wind  velocity  over  a  verti¬ 
cal  mixing  length 

Determines  horizontal  wind  field  above  100  km 
altitude 

Computes  the  linear  combination  of  two  vectors 
each  multiplied  by  an  arbitrary  constant 

Computes  the  product  of  a  vector  with  a 
scalar 

Computes  the  volume  of  a  given  geometrical 
region 

Determines  water  vapor  concentration  below 
45  km  altitude 

Computes  the  beta  energy  deposition  coeffi¬ 
cient,  particle  release  rate,  and  average 
energy 

Computes  the  delayed  gamffla-ray  energy  depo¬ 
sition  parameters 

Provides  prompt  gamma  energy  deposition 
parameters 

Initializes  the  gamma-ray  energy  deposition 
parameters 

Provides  the  time -dependent  neutron  energy 
deposition  parameters 
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Table  1  (Continued) 

Name  of  Routine 

Description 

WONP 

Provides  prompt  neutron  energy  deposition 
parameters 

IVONl 

Initializes  the  neutron  energy  deposition 
parameters 

WOXC 

Computes  the  X-ray  energy  containment 

KOXP 

Computes  the  prompt  X-ray  energy  deposition 
integral 

IVOX1 

Initializes  the  neutron  energy  deposition 
parameters 

XMAG 

Computes  the  absolute  magnitude  (length)  of 
a  vector 

ZTTOUT 

Converts  a  local  time  and  Gregorian  calendar 
date  to  the  local  time  and  Gregorian  calendar 
date  at  Greenwich 
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SECTION  2 

ROUTINE  INTERCONNECTIONS  AND  COMMON  BLOCKS 


Table  2  lists  HTSCOM  routines  and  library  routines  called  directly 
and  the  common  blocks  used  by  each  routine.  Table  3  lists  routines  calling  a 
specific  routine  and  Table  4  lists  routines  using  a  specific  common  block. 
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Table  2.  WJSCOM  routines,  library  routines, 
and  common  blocks  used. 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

1  COMMON  BLOCK 

ROUTINES 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

CONWES 

ATMOS I 

ALTODN 

PHEN 

DATE 

ATMOSN 

PROP 

FOOTPC 

ATMOUP 

REGINT 

INPUT 

BREG 

SPECW 

IONOSI 

BURST 

SSYS 

MAGFIT 

CHEMR 

TIME 

PUENOrl 

C  INPUT 

TSYS 

SATC 

CIRDES 

WEDEPO 

SORB IT 

CONST 

WOG 

ZTTOUT 

DBREG 

WON 

DCREG 

DEPDAT 

DEVICE 

FBREG 

HEACOM 

HEAVE1 

JAMMER 

MAGLNK 

1  ADULT 
OPTION 
ORGIN 
PATH 

WOX 

ACCUM 

ALOG 

AMBABS 

ATMOS U 

ALOG 

ATMOSN 

GOX 

GWAT 

HTOS 

UNI  TV 

VECLIN 

VECM 

XMAG 

CONST 

ANLYT2 

EXP 

SORT 

ANTND 

ATMOSU 

ALOG 

ANTMP 

COS 

BETAB 

SIN 

BREG 

FBTINT 

SQRT 

C INPUT 

GEOCOR 

CLINK 

LOCLAX 

CONST 

Table  2  (Continued) 


ROUTINE 

WESCOM  ROUTINES 
CALLED  DIRECTLY 

LIBRARY  ROUTINES 
CALLED  DIRECTLY 

COMMON  BLOCKS  j 

USED 

ANTND 

MAGFLD 

FBREG 

(Continued) 

MATMUL 

MAGLNK 

NINT 

OPTION 

PROPEN 

PROP 

SUBVEC 

SSYS 

UNITV 

TSYS 

VECLIN 

VECM 

XMAG 

ARRLIM 

COORDV 

ALOG 

C INPUT 

VOLUME 

COS 

CONST 

EXP 

DBREG 

RAY 

ATMOS F 

ATMOSU 

COS 

BREG 

he; 

EXP 

CONST 

HELP 

SIN 

HEACOM 

LAGRAN 

SQRT 

HEAVE1 

PEDEP 

FBREG 

SUBVEC 

UN  I  TV 

VECM 

XMAG 

ATMOS I 

BKRIAD 

ALOG 

ALTODN 

EXPINT 

COS 

ATMOUP 

JULIAN 

EXP 

CONST 

OZONE 

SIN 

TIME 

SOLCYC 

SQRT 

ZHTEMP 

SOLORB 

SOLZEN 

TEMPZH 

WATER 

ATMOSU 

BKRIAD 

EXP 

ALTODN 

EXPINT 

ATMOSN 

ATMOUP 

CONST 

TIME 

ZHTEMP 

AZELR 

LOCLAX 

AT  AN  2 

CONST 

MATMUL 

SEPA 

SUBVEC 

XMAG 

BALFIT 

MAGDIS 

SQRT 

CONST 

LENGTH 

BB 

ALOGIO 

BDPSK 

NORMAL 

EXP 

CONST 

SQRT 

DIGMOD 
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Table  2 


Continued 


ROUTINE 


BETAS 


WESCOM  ROUTINES 
CALLED  DIRECTLY 


LIBRARY  ROUTINES 
CALLED  DIRECTLY 


BFIELD 


BKRIAD 


BLKDTA 


BTRACE 


CARDIN 

CHEMD 


CHEMDQ 


CHEMEF 


ATMOSU 

CMDEDT 

DEDEP 

DRATE 

DTNEP 

DTNEQ 

EN2VIB 

INITAL 

IONOSU 

PHOTOR 

SPECDP 

SPECDQ 

VECM 

ATMOS U 

DRATE 

DTNEQ 

IONOSU 

SPECDO 

CHMION 

GNNO 


COMMON  BLOCKS 
USED 


ANTMP 
BETAB 
BREG 
BURST 
CINPUT 
C 


ATMOSN 

ATMOST 

BEDIT 

BREG 

CHMOVL 

CINPUT 

FDSRAT 

OPTION 

OUT 


ATMOSN 

ATMOST 

CINPUT 

FDSRAT 

CHEMAN 

SPECEF 
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ROUTINE 

WESCOM  ROUTINES 
CALLEO  DIRECTLY 

I  LIBRARY  ROUTINES 
CALLED  DIRECTLY 

COMMON  BLOCKS 
USED 

CHEMEF 

'Continued) 

RATE 

TEXK 

SQRT 

CHEMHR 

ATMOS U 

EQLAIR 

EQLMTL 

PHOTOR 

RATE 

ALOG 

EXP 

SQRT 

ATMOSN 

CHEMB 

FDSRAT 

OPTION 

SPEC 

CHEMQ 

RATE 

SQRT 

SPECQ 

CHMION 

ANALYT2 

RICATT 

EXP 

SQRT 

CHEMAN 

SPECEF 

CHXSPC 

ALOG 

EXP 

WEDEPO 

CLOSE 

DOT 

VECLIN 

CMDEDT 

PEDEP 

PHOTOD 

EXP 

SQRT 

ATMOST 

BEDIT 

BREG 

OPTION 

OUT 

WRATE 

CMFEDT 

ATMOS F 

PIONF 

BREG 

EDTOVL 

OPTION 

SOURCE 

SPECEF 

COMB 

ALOG 

EXP 

SQRT 

CONST 

CONSPC 

SQRT 

WEDEPO 

COORDV 

COS 

SIN 

CONST 

COORDX 

ASIN 

COS 

SIN 

SQRT 

CONST 

CORDTF 

AZELR 

COORDV 

COORDX 

GEOCOR 

LOCLAX 

MATMUL 

VECLIN 

XMAG 

CONST 

ERFC 

SQRT 

DIGMOD 
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Table  2 


Continued 


ROUTINE 

WESCOM  ROUTINES 
CALLED  DIRECTLY 

LIBRARY  ROUTINES 
CALLED  DIRECTLY 

COMMON  BLOCKS 

USED 

DATE 

ORIGIN 

DBMHT 

ATMOSU 

EXP 

BREG 

HELP 

SQRT 

BURST 

MAGDIS 

CONST 

MAGFLD 

DBREG 

VECM 

DEVICE 

FBREG 

LENGTH 

OPTION 

PHEN 

DBMLT 

ATMOSU 

EXP 

BURST 

HELP 

SQRT 

CONST 

UNITV 

DEVICE 

VECM 

DBREG 

VOLUME 

FBREG 

WOBD 

PHEN 

DBSTAM 

DATE 

BREG 

GEOCOR 

C INPUT 

MAGFLD 

CONST 

PSDM 

DBREG 

VECM 

OPTION 

XMAG 

ORIGIN 

ZTTOUT 

PHEN 

TIME 

DEBRIS 

AZELR 

ALOG 

BREG 

CHXSPC 

COS 

CONST 

CONSPC 

SIN 

DTUBE 

FZET 

SQRT 

MAGLNK 

GEOCOR 

OPTION 

LEKSPC 

WEDEPO 

MAGFLD 

SEPA 

VECM 

DECODE 

DEDEP 

EDEPB 

COS 

BREG 

EDEPG 

SIN 

BTUBE 

EDEPND 

CINPUT 

GEOCOR 

DLREG 

VECM 

FBREG 

XMAG 

MAGLNK 

OPTION 

DELPSK 

ERFC 

ATAN 

CONST 

SQRT 

DIGMOD 

DIFINT 

ALOG 

EXP 

DISTF 

EXP 

DOT 
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Table  2 


(Continued) 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

DPSK 

ERFC 

EXP 

CONST 

SQRT 

DIGMOD 

DQPSK 

ERFC 

QFNTN 

SQRT 

DIGMOD 

DRATE 

RATE 

EXP 

WRATE 

SQRT 

WRATEH 

DSTINT 

PCLOUD 

C INPUT 

PGROUP 

CONST 

RINTER 

DCREG 

UN  I  TV 

FBREG 

VECM 

OPTION 

PROP 

DTNEP 

ALOG 

ATMOST 

EXP 

FDSRAT 

SQRT 

WRATE 

DTNEQ 

HELP 

ALOG 

ATMOST 

SQkV 

FDSRAT 

OUT 

WRATE 

DUSTMI 

HELP 

ALOG 

BREG 

SOIL 

BURST 

DCREG 

PHEN 

EANOM 

DCOS 

DMOD 

DSIN 

EDEPB 

ATMOSU 

ALOG 

BREG 

AZELR 

COS 

BURST 

BFIELD 

EXP 

CONST 

COORDX 

SIN 

DBREG 

CROSS 

SQRT 

FBREG 

EIF 

E2 

GEOCOR 

LOCLAX 

MAGFLD 

MATMUL 

PMASS 

RINTER 

SUBVEC 

UN  I  TV 

VECLIN 

VECM 

VOLUME 

WOBD 

XMAG 

MAGLNK 
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ROUTINE 

WESCOM  ROUTINES 
CALLED  DIRECTLY 

LIBRARY  ROUTINES 
CALLED  DIRECTLY 

COMMON  BLOCKS 

USED 

EDEPG 

ATMOSU 

ACOS 

BURST 

AZELR 

ALOG 

CONST 

COORDX 

ASIN 

DBREG 

DISTF 

COS 

FBREG 

HELP 

EXP 

MAGFLD 

SIN 

PMASS 

SQRT 

RINTER 

SUBVEC 

UNITV 

WOGD 

XMAG 

EDEPND 

PMASS 

EXP 

BREG 

SUBVEC 

SQRT 

BURST 

WOND 

DEVICE 

XMAG 

EIF 

ALOG 

EXP 

ELDEN 

VECM 

BREG 

CHEMD 

BURST 

ENEF 

CINPUT 

CHMOVL 

DEVICE 

OPTION 

SOURCE 

ENEF 

ATMOS F 

ALOG 

ATMOSN 

ATMOSU 

EXP 

BREG 

CHEMEF 

SQRT 

CHMOVL 

CHEMQ 

CINPUT 

DEPEP 

CONST 

IONOSU 

HEACOM 

PIONF 

OPTION 

TEXK 

SOURCE 

UNITV 

SPECEF 

VECM 

SPECQ 

EN2VIB 

EXP 

ATMOST 

EQCROS 

HELP 

SQRT 

CONST 

LENGTH 

EQLAIR 

SQRT 

SPEC 

EQLMTL 

SQRT 

SPEC 

ERF 

ERFC 

SQRT 

CONST 

ERFC 

ERF 

EXP 

CONST 

SQRT 

EUXFIT 

WOXC 

ALOG 

BREG 

BURST 

EXINTP 
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Table  2  (Continued) 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

EXPINT 

ALOG 

EXP 

E2 

EIF 

EXP 

SQRT 

FAC 

ALOG 

CONST 

EXP 

SQRT 

FBEN 

ATMOS F 

ALOG 

BREG 

ATMOSU 

COS 

CHEMB 

BFIELD 

EXP 

CHEMFB 

SIN 

C INPUT 

CHEMHR 

SQRT 

CONST 

FBREG 

DIFINT 

LENGTH 

GEOCOR 

MAGLNK 

HDINT 

LOCLAX 

MAGDIS 

MAGFLD 

MATMUL 

RATE 

SEPA 

SUBVEC 

TFUNC 

UN  I  TV 

VECLIN 

VECM 

XMAG 

OPTION 

FBI  NT 

1 

ALOG 

BREG 

COS 

BTUBE 

DOT 

SIN 

BURST 

FBEN 

SQRT 

CHEMB 

GEOCOR 

CHEMFB 

LOCLAX 

C  INPUT 

MAGFLD 

CONST 

MATMUL 

DBREG 

RINTER 

DEVICE 

SEPA 

FBREG 

STRII 

MAGLNK 

SUBVEC 

OPTION 

UN  I  TV 

PROP 

VECLIN 

REGINT 

VECM 

STRIIN 

XMAG 

FBITMD 

EXP 

BREG 

SIN 

ALOG 

COS 


EXP 


FBMHI 


ATMOS F 
ATMOSU 
ERF 


BREG 

BURST 

CONST 


Table  2  (Continued) 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

FBMHI 

EUXFIT 

SIN 

PHEN 

(Continued) 

FBITMD 

FBMHT 

IONLEK 

LAGRAN 

LOSCON 

MAGFLD 

WOXC 

SQRT 

FBMHT 

BTRACE 

ACOS 

BREG 

COORDX 

ALOG 

CONST 

EQCROS 

COS 

FBREG 

FBITMD 

EXP 

LENGTH 

LOCLAX 

SIN 

PHEN 

MAGDIS 

MATMUL 

SEPA 

STRIP 

VECLIN 

VECM 

XMAG 

SQRT 

TRACE 

FBMLI 

ATMOSF 

ALOG 

BREG 

ATMOSU 

EXP 

BURST 

WOXC 

SQRT 

CONST 

PHEN 

FBMLT 

ATMOSU 

EXP 

BREG 

COORDX 

SQRT 

CONST 

LOCLAX 

FBREG 

MATMUL 

VECLIN 

XMAG 

PHEN 

FBTINT 

ACCUM 

ALOG 

ANTilP 

DOT 

COS 

BREG 

FBEN 

SIN 

BURST 

GEOCOR 

SQRT 

BTUBE 

MAGFLD 

CHEMB 

SUBVEC 

DEVICE 

UN  I  TV 

FBREG 

VECLIN 

MAGLNK 

VECM 

OPTION 

XMAG 

FNNO 

EXP 

FOOTPC 

CROSS 

ALOG 

BREG 

FBEN 

COS 

CHEMB 

FBINT 

SIN 

CHEMFB 

FPINT 

CINPUT 

FPLNK 

CONST 

FPPLOT 

FBREG 

GEOCOR 

FPINTC 

HTOS 

LENGTH 
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Table  2  (Continued) 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

FOOTPC 

LOCLAX 

MAGLNK 

(Continued) 

MAGDIS 

OPTION 

MAGFLD 

PATH 

RINTER 

PROP 

SEPA 

REG I NT 

SUBVEC 

SSYS 

TLTFLD 

UN  I  TV 

VECLIN 

VECM 

XMAG 

TSYS 

FPCIR 

BDPSK 

ALOGIO 

CONST 

CPSK 

ATAN2 

DIGMOD  i 

DELPSK 

SIN 

FPCIRC 

DPSK 

DQPSK 

FSKWB 

MODEM1 

MODEM2 

MODEM3 

NFSKMF 

USERMD 

SQRT 

MODULT 

FPINT 

GEOCOR 

ALOG 

CONST 

■ 

HELP 

ALOGIO 

FPINTC 

HTOS 

SUBVEC 

UN  I  TV 

VECLIN 

OPTION 

FPLNK 

AZELR 

ALOG 

ANTMP 

FBTINT 

ALOGIO 

BREG 

FPCIR 

SIN 

C INPUT 

GAIN 

SQRT 

CLINK 

LOCLAX 

CONST 

MAGFLD 

FBREG 

MATMUL 

FPCIRC 

NINT 

MAGLNK 

SLOBE 

OPTION 

SUBVEC 

PROP 

UN  I  TV 

SSYS 

VECLIN 

VECM 

XMAG 

TSYS 

FPPLOT 

COS 

CONST 

FPINIC 

OPTION 

FQAMB 

ATMOSU 

SQRT 

ATMOS N 

DRATE 

ATMOST 

DTNEQ 

FDSRAT 

OUT 

Table  2  (Continued) 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

FSKWB 

FSKWBY 

SQRT 

CONST 

DIGMOD 

FSKWBQ 

EXP 

FSKWBY 

COMB 

ALOG 

DIGMOD 

FSKWBO 

EXP 

SQRT 

FZET 

ALOG 

SQRT 

GAIN 

DOT 

SIN 

SEPA 

SQRT 

1 

SEPA 

ATAN2 

CONST 

FNNO 

ALOG 

EXP 

GOX 

EXP 

GWAT 

GWPLME 

BF I  ELD 

ACOS 

BREG 

COORDV 

ASIN 

CONST 

DATE 

COS 

FBREG 

GEOCOR 

SIN 

ORIGIN 

MAGFLD 

SQRT 

PHEN 

SEPA 

RAY 

SUBVEC 

TIME 

UVWD 

VECM 

XMAG 

ZTTOUT 

WSTOR 

HDINT 

MAGDIS 

ALOG 

CONST 

EXP 

LENGTH 

SQRT 

HEAVE 

ATMOSU 

EXP 

HEAVE1 

LAGRAN 

PHEAVE 

SQRT 

HELP 

EXIT 

HERCAN 

EXP 

HPCIIEM 

SPECEF 

TSTEF 

HTOS 

SEPA 

COS 

CONST 

XMAG 

SQRT 

IDEFLT 

COORDV 

BURST 

C INPUT 

CIRDES 

CONST 

JAMMER 

MODULI _ 
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ROUTINE 


WESCOM  ROUTINES 
CALLED  DIRECTLY 


LIBRARY  ROUTINES 
CALLED  DIRECTLY 


COMMON  BLOCKS 
USED 


IDEFLT 

OPTION 

(Continued) 

ORIGIN 

PATH 

SSYS 

TSYS 

rN ITAL 

DRATE 

EXP 

1  Cat'll 

EIF 

u  jy  ■ 

INPUT 

BLKDTA 

ALOGIO 

BURST 

CARDIN 

EXIT 

C INPUT 

DECODE 

CIRDES 

IDEFLT 

CONST 

INPUTP 

DEVICE 

ORBIT 

INPTCM 

SPCTRM 

JAMMER 

TEST 

MODULT 

WOG1 

OPTION 

WONT 

ORIGIN 

WOX1 

PATH 

SPECW 

SSYS 

TSYS 

INPUTP 

BB 

BURST 

COORDV 

C INPUT 

CORDIF 

CIRDES 

GEOCOR 

CONST 

RRMSTP 

DEVICE 

XMAG 

INPTCM 

JAMMER 

MODULT 

OPTION 

ORIGIN 

PATH 

SSYS 

TSYS 

IONLEK 

IONOSI 

COORDV 

ALOG 

CONST 

FQAMB 

COS 

IONOUP 

JULIAN 

EXP 

MAGLNK 

QCOSM 

SEPA 

SOLCYC 

SOLORB 

ZTTOUT 

SQRT 

TIME 

IONOSU 

BFZ 

EXP 

CONST 

COORDV 

SQRT 

IONOUP 

FQAMB 

MAGLNK 

QCOSM 

TIME 

Table  2  (Continued) 


ROUTINE 

WESCOM  ROUTINES 
CALLED  DIRECTLY 

LIBRARY  ROUTINES 
CALLED  DIRECTLY 

COMMON  BLOCKS 
USED 

IONOSU 

RATE 

(Continued) 

SEPA 

JAMND 

AZELR 

SQRT 

CLINK 

GAIN 

CONS! 

PROPEN 

JAMMER 

SLOBE 

OPTION 

SUBVEC 

PROP 

UN  I  TV 

SSYS 

VECM 

TSYS 

XMAG 

JULIAN 

TIME 

LAGRAN 

ALOG 

HEAVE  1 

EXP 

SQRT 

LEKSPC 

SQRT 

WEDEPO 

LOCATE 

LOCLAX 

CROSS 

UN  I  TV 

VECM 

LOSCON 

ALOG 

SQRT 

MAGDIS 

ALOG 

CONST 

EXP 

LENGTH 

SQRT 

MAGFIT 

COORDV 

ASIN 

CONST 

*r : 

ATAN 

MAGLNK 

■  Ex  1 1  1 

COS 

VECM 

SIN 

SQRT 

MAGFLD 

DOT 

SQRT 

MAGLNK 

UN  I  TV 

VECLIN 

XMAG 

MANPTS 

ATMOSU 

SIN 

BREG 

DOT 

SQRT 

BTUBE 

GEOCOR 

BURST 

LOCLAX 

C  INPUT 

MAGFLD 

CONST 

MATMUL 

DBREG 

RINTER 

DEVICE 

SUBVEC 

DTU3E 

UNITV 

MAGLNK 

VECLIN 

OPTION 

VECM 

PROP 

XMAG 

REGINT 

MATMUL 

Table  2  (Continued) 


ROUTINE 


MISCAT 


MODEM! 


MODEM2 


MODEM3 


NFSKMF 


NORMAL 

OMEGA 


OZONE 


PCHEM 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

CABS 

CCOS 

COS 

CSIN 

SIN 

LOCATE 

ALOG 

EXP 

LOCATE 

ALOG 

EXP 

LOCATE 

ALOG 

EXP 

COMB 

SANCAN 

EXP 

AZELR 

ACOS 

DOT 

ASIN 

GAIN 

ATAN2 

LOCLAX 

COS 

MATMUL 

OMEGA 

RINTER 

SEPA 

SUBVEC 

SIN 

UNITY 

VECLIN 

VECM 

VOLUME 

XMAG 

EXP 

DOT 

COS 

HELP 

SIN 

LOCLAX 

MATMUL 

SQRT 

SUBVEC 

XMAG 

SQRT 

EANOM 

ACOS 

TATSRT 

ASIN 

ATAN 

COS 

SIN 

SQRT 

TAN 

EXP 

SIN 

TEXK 

EXP 

SORT 

COMMON  BLOCKS 
USED 


CONST 

MAGLNK 


CONST 

ORIGIN 

PATH 
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Table  2 


(Continued) 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

PCHEM 

SPECEF 

(Continued) 

TSTEF 

PCLOUD 

LOCLAX 

BREG 

flATMUL 

CONST 

DCREG 

FBREG 

PEDEP 

EXP 

BREG 

SQRT 

BURST 

f  %  1  "  I ' T  ^ ;'Ki'  *  f  ! 

CONST 

WOGP 

DEVICE 

WONP 

DTUBE 

WOXP 

XMAG 

PGROUP 

MI SCAT 

ALOG 

CONST 

EXP 

DCREG 

PHEAVE 

EXP 

HEAVE  1 

SQRT 

PHENOM 

ATMOS F 

COS 

BREG 

ATMOS  I 

SIN 

BURST 

1  fna 

CINPUT 

1  !aHl  *> 

CONST 

COORDX 

DBREG 

DATE 

DCREG 

DBMHT 

DEVICE 

DBMLT 

FBREG 

LENGTH 

MAGLNK 

FBMHI 

OPTION 

FBMHT 

ORIGIN 

FBI1LI 

PHEN 

FBMLT 

GEOCOR 

GWPLME 

HELP 

MAGDIS 

MAGFLD 

MAGFIT 

RADOUT 

XMAG 

VECM 

ZTTOUT 

TIME 

PHOTOD 

AZELR 

ALOG 

BREG 

EXP 

SIN 

SQRT 

PHOTOR 

SUBVEC 

BREG 

XMAG 

FBREG 

FDSRAT 
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Table  2 


(Continued) 


ROUTINE 

WESCOM  ROUTINES 
CALLED  DIRECTLY 

LIBRARY  ROUTINES 
CALLED  DIRECTLY 

COMMON  BLOCKS 

USED 

PIONF 

ATMOS F 

ALOG 

BREG 

ATMOSU 

EXP 

BURST 

AZELR 

SIN 

CONST 

DEBRIS 

SQRT 

DEPDAT 

EUXFIT 

DEVICE 

HPCHEM 

OPTION 

PCHEM 

SPECEF 

PMASS 

TSTEF 

SUBVEC 

VECLIN 

VECM 

WOXP 

XMAG 

PLINTP 

ALOG 

PLYVAL 

PMA5S 

ATMOSU 

EXP 

CONST 

SQRT 

PROPEN 

COS 

CINPUT 

SIN 

CONST 

CMFEDT 

SORT 

EDTOVL 

DSTINT 

MAGLNK 

FBINT 

OPTION 

HTOS 

PROP 

MANPTS 

REGINT 

PTH I NT 

START 

SUBVEC 

UNITV 

VECLIN 

VECM 

XMAG 

PSDM 

ARRLIM 

COS 

CONST 

UVWAV 

SPEEDS 

UVWD 

RAY 

WSTOR 

PTH  I  NT 

ACCUM 

ALOG 

CINPUT 

BFIELD 

SQRT 

CONST 

DOT 

MAGLNK 

ELDEN 

OPTION 

GEOCOR 

PROP 

MAGFLD 

REGINT 

SUBVEC 

UMTV 

Vt  CLIN 

XMAG 

QCOSM 

COS 

CONST 

EXP 

SIN 
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Table  2  (Continued) 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  OIRECTLY 

CALLED  DIRECTLY 

USED 

QFNTN 

ALOGIO 
!  EXP 

RADOUT 

HELP 

ALOG 

BREG 

EXP 

FBREG 

PHEN 

RATE 

EXP 

CHEMR 

RICATT 

EIF 

ALOG 

EXP 

RINTER 

CCS 

CONST 

■TV*  S  ^fl 

SIN 

MAGLNK 

DOT 

GEOCOR 

LOCLAX 

MAGFLD 

MATMUL 

SUBVEC 

TOROID 

TUBE 

UNITV 

VECLIN 

VECM 

XMAG 

SORT 

RGEOM 

RRMSTP 

PLYVAL 

ALOGIO 

SANCAN 

ERF 

EXP 

CONST 

FAC 

HERCAN 

SQRT 

SATC 

AN7ND 

1 

CHEMB 

ATMOS I 

CIRDES 

DECODE 

,  j 

CLINK 

GEOCOR 

JAMMER 

HTOS 

OPTION 

i 

IONOSI 

PATH 

, 

JAMND 

PROP 

MAGFIT 

SSYS 

PROPEN 

TIME 

SATCIR 

SATLNK 

SLINK 

SUBVEC 

TORBIT 

UNI  TV 

VECLIN 

VECM 

XMAG 

STSYS 

SATCIR 

BDPSK 

ALOGIO 

C INPUT 

CPSK 

ATAN2 

CONST 

DELPSK 

SIN 

DIGMO  ■ 

DPSK 

SQRT 

■KMH 

I 
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Table  2 


ROUTINE 


SATCIR 

(Continued) 


SATLNK 


SLINK 


SOLCYC 


SOLORB 


SOLZEN 


SORB IT 


WESCOM  ROUTINES 
CALLED  DIRECTLY 


DQPSK 

FSKWB 

HODEM1 

MODEM2 

M0DEM3 

NFSKMF 

USERMD 


LIBRARY  ROUTINES 
CALLED  DIRECTLY 


DECODE 


COORDV 

EANOM 

TATSRT 


ACOS 

COS 

SIN 

SQRT 

ASIN 

ATAN2 

COS 

EXIT 

SIN 

SQRT 


COMMON  BLOCKS 
USED 


OPTION 

SIGCOM 


SPECDQ 


C INPUT 

CLINK 

CONST 

OPTION 

PROP 

SIGCOM 

SSYS 

TSYS 

CONST 

CIRDES 

CLINK 

SSYS 

TSYS 

CONST 

DCREG 

ORIGIN 

TIME 

CONST 

CONST 

TIME 

CONST 

TIME 


CONST 

OPTION 

PATH 


SPECW 

ATMOS N 
ATMOST 
WRATE 
WRATEH 

ATMOST 

WRATE 


Table  2  (Continued) 


ROUTINE 


TEMPZH 


TFUNC 

TLTFLD 

TORBIT 


TOROID 


WESCOM  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

FBREG 

PROP 

REG I NT 

ACCUM 

ALOG 

CONST 

CROSS 

SQRT 

PROP 

DOT 

HELP 

SUBVEC 

STRIIN 

UNITV 

VECLIN 

SQRT 

BREG 

FBREG 

ATAN 

CONST 

SQRT 

TAN 

ALTODN 

CONST 

TIME 

ZHTEMP 

ALOG 

EXP 

SQRT 

MATMUL 

VECLIN 

XMAG 

CONST 

COORDX 

C INPUT 

LOCLAX 

CONST 

MATMUL 

VECLIN 

VECM 

XMAG 

PATH 

DOT 

LOCLAX 

MATMUL 

VECLIN 

VECM 

SQRT 

RGEOM 

DOT 

MAGFLD 

SUBVEC 

RGEOM 

Table  2  (Continued) 


WESCOK  ROUTINES 

LIBRARY  ROUTINES 

COMMON  BLOCKS 

ROUTINE 

CALLED  DIRECTLY 

CALLED  DIRECTLY 

USED 

TUBE 

(Continued) 

HU  I 

UNITV 

SQRT 

USERMD 

DIGMOD 

UVWAV 

RAY 

SPEEDS 

UVWD 

COS 

CONST 

EXP 

RAY 

WSTOR 

VECLIN 

VECM 

VOLUME 

HELP 

COS 

CONST 

SIN 

SQRT 

WATER 

ALOG 

CONST 

EXP 

TIME 

SIN 

WOBD 

DEVICE 

WOGD 

EXINTP 

DEVICE 

WOG 

WOGP 

EXINTP 

WOG 

WOG1 

DEVICE 

SPECW 

WOG 

WOND 

EXINTP 

PLINTP 

WON 

WONP 

EXINTP 

PLINTP 

WON 

WON1 

ALOG 

DEVICE 

EXP 

SPECW 

SQRT 

WON 

WOXC 

EXINTP 

ALOG 

WOX 

EXP 

WOXP 

EXINTP 

WOX 

WOX1 

ALOG 

DEVICE 

EXP 

SPECW 

WOX 

XMAG 

SQRT 

ZTTOUT 

EXIT 

CONST 

TIME 
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Table  3.  WESCOM  routine  called  by  a 
specific  routine. 


SUBROUTINE 

CALLED  BY 

AMBABS 

PROPFN 

ANLYT? 

CHNION 

ANTNO 

SATC 

ARRLIM 

P5DM 

ATMOSF 

CMFFDT 

ENTF 

FBf  N 

FBfHI 

FOMLI 

PHENOM 

PIONF 

ATMOSI  CUN WES 

PHFMOM 
SATC 

ATKOSU  AMLAB5 

ANTE  D 
ATMOSF 
CHFMD 
CHF  MBS 
CHEMHR 
DBMMT 
dbmlt 
EDFPM 
EDEPG 
ENEF 
FBFN 
FBMHl 
FRNLl 
FBF'lT 
FQAMB 
HEAVE 
MANPT5 
PI  CNF 
.  PM A  55 


SUBROUTINE 

CALLED  BY 

A2ELP 

ANTND 

COPDTF 

DEPRIS 

EDLPB 

EDEPG 

FPLNF 

JAMND 

NINT 

PHENOM 

PHCTOD 

PIONF 

PROPEN 

SATLflK 

BALFIT 

6TPACE 

OB 

INPIJTP 

BDPSK 

F  PC  IP 
SATCIR 

BETAS 

Af  o 

BFIELO 

EDEPB 

FBFN 

r  13 1  NT 

GwPLMt 

PHENOM 

FTHINT 

PINTER 

BFZ 

BKPIAD 

IONOSU 

BKPJAD 

AT NO 5 I 
ATMOSU 

BLKDTA 

INPUT 

BTRACE 

FbPHT 

SUBROUTINE 

CALLED  BY 

CARDIN 

INPUT 

CHEMD 

EL  or  K 

CHfMDO 

BETAS 

FBFN 

CHEMEF 

ENFF 

CHEMHP 

FBFN 

CHEMQ 

FNf  F 

CHMON 

CHFMEF 

CHX5PC 

DEBRIS 

CLOSC 

RINTFR 

CHDEDT 

CHEMD 

cmfedt 

PROPFN 

COMB 

FSK'WPY 

NFSKMF 

CONSPC 

DEF'P  I  S 

COOROV 

ARRl.IM 
COR DTP 
GWF'LME 
IDFFIT 
INPIJTP 
IOMOSI 
IONOSU 
MAGFIT 
SOFPIT 

SUBROUTINE 

1  CALLED  BY 

SUBROUTINE 

CALLED  BY 

SUBROUTINE 

CALLED  BY 

COORDX 

CORDTF 

EOF  PR 

DOT 

EIF 

EDFPP 

E2 

EDfPG 

fbtint 

INI TAL 

FPMHT 

GAIN 

PICATT 

frmlt 

MAGFLD 

PHLNOM 

MANPTS 

ELDEN 

PTHINT 

Tore  IT 

MINT 

OMEGA 

CORDTF 

INPUTP 

PTHINT 

RIOTER 

ENEF 

ft  pf  r 

CPSK 

FPCIR 

SEf  A 

STPII 

EN2Vlb 

CHEMD 

SATCIR 

TOPOID 

TUt'E 

EOCROS 

FBHHT 

CROSS 

EDf  PR 

FOOTPC 

LOCLAX 

5TR  1 1 

DPSK 

FPCIR 

SATCIR 

EQLAIR 

CHEMHR 

DATE 

CONrfES 

DOPSK 

FPCIR 

SATCIR 

eqlmtl 

CHEMHR 

DBSTAM 

GWPLME 

PHLNOM 

dratl 

CHfMD 

CHEMDO 

ERF 

FRFC 

FQMHI 

SANCAN 

dbmht 

PHLNOM 

FOAMP 

INI TAL 

ERFC 

CPSK 

DELPSK 

dbmlt 

PMRNOM 

DSTINT 

PR OPEN 

DPSK 

DGP5K 

ERf 

DB5TAM 

PHLNOM 

dtnep 

CHEMD 

EUXFIT 

FBMHI 

DEBRIS 

PEDEP 

PIONT 

FIONF 

DTNEO 

CHI  MO 

CHEMDO 

FQAMB 

EXINTP 

WOGD 

DECODE 

INPUT 

WOGP 

SAT  C 

WOND 

SLINK 

DUSTMI 

PHFNOM 

WONP 

woxc 

DEDEP 

CHEMD 

FNEF 

EANOM 

ORPIT 

WOXP 

FBEN 

SORB  IT 

EXPINT 

ATMOS  I 
ATMOSU 

DELPSK 

FPCIR 

SATCIR 

EDEPR 

DEDEP 

E2 

EDEPR 

EDEPG 

FEOFF 

FAC 

SANCAN 

DIFINT 

FBEN 

EDEPND 

DEDEP 

FBEN 

FBI  NT 

DISTF 

EDfPG 

FBTINT 

FOOTPC 

SUBROUTINE 


CALLED  BY 


SUBROUTINE 


CALLED  BY 


SUBROUTINE 

CALLED  BY 

SUBROUTINE 

CALLED  BY 

SUBROUTINE 

CALLED 

FBINT 

FOOTPC 

GAIN 

FPLNK 

hercan 

SANCAN 

pporpn 

JAMND 

NINT 

HPCHFM 

PIONr 

FBITMD 

FBMHI 

FBMHT 

GEOCOR 

satlnk 

ANTND 

CORDTF 

HTOS 

AMBABS 

FOOTPC 

fbmhi 

PHE  NOM 

DBSTAM 

FPINT 

PROPEN 

DEBRIS 

SATC 

FBMHT 

FBMHI 

PHFNOM 

DEDEP 

EDFPR 

FBEN 

FBINT 

IDEFLT 

INPUT 

FBMLI 

PHENOM 

FBTINT 

FOOTPC 

inital 

CHE  MO 

fbmlt 

PHE  MOM 

FPINT 

GWPLME 

INPUTP 

INPUT 

coni'  f  r. 

FBTINT 

ANTND 

MAf  PTS 

INPUTP 

INPUT 

FPLNK 

PHENOM 

PTH1NT 

ionlek 

FBMHI 

FNNO 

GNNO 

RINTER 

SAT  C 
satlnk 

IONOSI 

CONWLS 

FOOTPC 

CONWE.S 

GNNO 

CHFMEF 

SATC 

FPCIR 

FPl.f'K 

GOX 

AMPABS 

IONOSU 

CHFMD 
CHE MOO 

FPINT 

FOOTPC 

GWAT 

AMPABS 

ENEF 

FPLNK 

FOOTPC 

gwplme 

PHENOM 

JAMND 

SATC 

FPPLOT 

FOOTPC 

HDINT 

FBEN 

JULIAN 

ATMOSI 

IONOSI 

FQAMB 

IONOSI 

IUNOSU 

HEAVE 

ATMOSF 

lagran 

ATMOSF 

FBMHI 

FSKWB 

FPCIR 

HELP 

ATMOSK 

HEAVE 

SATC1R 

DBFHT 

DBMLT 

DTNEO 

LEKSPC 

DEBRIS 

F5KWBQ 

FSKWPY 

DUSTMI 

FDFPG 

LOCATE 

MOrEMl 

F5KWBY 

FSKWB 

FOCRCS 

MODEM2 

FPINT 

OMEGA 

MODEM 3 

F2ET 

DEBRIS 

PrtENOM 

RADOUT 

SOIL 

STRI  I 
VOLUME 
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Table  3  (Continued) 


SUBROUTINE 

CALLED  BY 

SUBROUTINE 

CALLED  BY 

SUBROUTINE 

LOCLAX 

ANTND 

AZLLR 

MANPTS 

PPOPEN 

PCLOUD 

CORDTF 

MATMUL 

ANTND 

PEDEP 

EDEPB 

AZELR 

FBEM 

COPDTF 

FBI  NT 

EDEPB 

PGROUP 

fbmht 

FBLN 

FB^  LT 

FBI  NT 

FOOTPC 

F  BMHT 

PHEAVE 

FPL  MX 

FBMLT 

MANPTS 

FPLNK 

PHENOM 

N I  NT 

MANPTS 

OMEGA 

NINT 

PHOTOD 

PCLOUD 

OMf  GA 

RINTER 

PCLOUD 

TOPRIT 

RINTER 

PHOTOR 

TOPOID 

TLTFLD 

TOPBIT 

LOSCON 

FUMHI 

TOROID 

PIONF 

MISCAT 

PGRCUP 

MAGDIS 

balfit 

BTRACE 

MODEM  1 

FPCIR 

PLINTP 

DBMHT 

FBEN 

SATC  I R 

j 

FBMHT 

FOOTPC 

HDINT 

M0DEM2 

FPCIR 

SATCIR 

PLYVAL 

PHENOM 

MODEM3 

FPCIP 

PMASS 

MAGFI T 

CONWES 

PHFNOM 

SATCIR 

SATC 

NFSKMF 

FPCIP 

5ATCIR 

PROPEN 

MAGFLD 

ANTNP 

DBMHT 

NINT 

A*  T'  r 

DBSTAM 

FPL*  * 

DEBRIS 

FDEPH 

normal 

BDPSK 

PSDM 

EDEPG 

FBFN 

PTHINT 

FBI  NT 

FBMHI 
f  BTINT 

OMEGA 

ONEMG5 

NINT 

MAGFIT 

QCOSM 

FOOTPC 

FPLNK 

OFNTN 

GWPLME 

MANPTS 

PHENCM 

PTHINT 

PINTFR 

ORBIT 

INPUT 

ATMOSI 

RADOUT 

OZONE 

TUFF 

PCHEM 

PIONF 

CALLED  BY 
DSTINT 


ATMOSF 

CMCEDT 


DSTINT 

HEAVE 


CONWES 


CMDEDT 


CHEMD 

CHEMHR 


CMFEDT 

ENEF 


WOND 

WONP 


PRH5TP 


EDLPB 
FDEPG 
EOF  PND 
PEDEP 
PIONF 

ANTND 
JAMND 
SAT  C 


DBSTAM 

PROPEN 


IONOSI 

IONOSU 

DQPSK 


PHFNOM 


Table  3  (Continued) 


SUBROUTINE 


RATE 


RICATT 


RINTER 


RRMSTP 


SANCAN 

SATC 


SATCIR 


SATLNK 


SEPA 


SLINK 

SLOBE 


SOIL 


CALLED  BY 


CHFMEF 

CHLMHR 

CHFMO 

DRATF 

FBEN 

IOMOSU 

CHMION 

DSTINT 

EDEPB 

FDFPG 

FBI  NT 

FOOTPC 

MANPTS 

NINT 

INPUTP 

NFSKMF 

CONWES 

SATC 

SATC 

azelr 

DEBRIS 

FBEN 

FBINT 

FBMHT 

FOOTPC 

GAIN 

GEOCOR 

GWPLME 

HTOS 

IONOSI 

IONOSU 

NINT 

SATC 

FPLNK 

JAMND 

SATLNK 

DUSTMI 


SUBROUTINE 

CALLED  BY 

SUBROUTINE 

SOLCYC 

ATMOSI 

IONOSI 

TATSRT 

SOLORB 

ATMOSI 

IONOSI 

TEMPZH 

solzfn 

ATMOSI 

TEST 

SORBIT 

CONWES 

TEXK 

SPCTRM 

INPUT 

TFUNC 

SPECDP 

CMEMD 

TLTFLD 

SPECDO 

CHEMD 

CHEMDQ 

TORBIT 

START 

PROPEN 

TOROID 

STRII 

FBINT 

TUBE 

STRIP 

FBMHT 

UN  I  TV 

SUBVEC 

ANTND 

ATMOSF 

AZELR 

EDEPB 

FDEPG 

FOFPND 

FBFN 

FBINT 

FBTINT 

FOOTPC 

FP1NT 

FPLNK 

GWPLME 

JAMND 

MANPTS 

NINT 

OMFGA 

PErEP 

PHOTOR 

PIONF 

PROPEN 

PTHINT 

RINTER 

SATC 

SATLNK 

STRII 

TUBE 


CALLED  BY 


ORBIT 

SORBIT 

ATMOSI 

INPUT 

CHEMEF 

ENFF 

PCHEM 

FBEN 

FOOTPC 

SATC 

RINTER 

PINTER 

AME.ABS 

ANTND 

ATMOSF 

DBMLT 

DSTINT 

EDEPB 

EDFPG 

FNEF 

FBFN 

FBINT 

FBTINT 

FOOTPC 

FPINT 

FPLNK 

JAMND 

LOCLAX 

MAGFIT 

MAGFLD 

MANPTS 

NINT 

PROPEN 

PTHINT 

RINTER 

SATC 

SATLNK 

STRII 

TUBF 
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Table  J  i continued; 


SUBROUTINE 

CALLED  BY 

SUBROUTINE  1 

1 

USERMD 

FPCIR 

SAT  C I R 

UVWAV 

PSIM 

UVWD 

GWPLME 

PSDM 

VECLIN 

AMBABS 

ANTND 

CLOSF 

CORDTF 

FDFPB 

FBFN 

FBINT 

FBMHT 

FBMLT 

FBTINT 

FOOTPC 

FPINT 

F  PLNK 

VOLUME 

MAGFLD 

MANPTS 

N  I  NT 

WATER 

PI  ONF 
PROPLN 

PTH  I  NT 
RINTLR 

WOBD 

SATC 

STR  I  I 
TLTFLD 

WOGD 

TORBIT 

TOROID 

WOGP 

WOG1 

VFCM 

AMBABS 

ANTND 

ATMOSF 

WOND 

CHL  MD 
DBMHT 

WONP 

DBMLT 
DBSTAM 
DEBP  I  S 

WON  1 

DEPLP 

DSTINT 

FDEPB 

FLOFN 

woxc 

E  Nf  F 

FBI  N 

FBINT 

WOXP 

FBMHT 

FBTINT 

FOOTPC 

WOXI 
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CALLED  BY 


SUBROUTINE 


CALLED  BY 


FPL  NK 

GWPLME 

JAMND 

LOCLAX 

MAGFIT 

MANPTS 

MINT 

PHFNOM 

P I ONF 

PROPEN 

RINTER 

SATC 

SATl NX 

TORBIT 

TOROID 

ARPLIM 
DBMLT 
LDF  PB 
MINT 

ATMOSI 

BETAS 
DBMLT 
f  DF  PB 

EDEPG 

PEDFP 

INPUT 

EOF PND 

Pt'DEP 

INPUT 

EUXFIT 
FBNHJ 
FBML  I 

PEIEP 
PI  ONF 

INPUT 


XMAG 


ZTTOUT 


AMHABS 

ANTND 

ATNOSF 

AZFLR 

CORDTF 

DBSTAM 

DtDFP 

FDEPB 

EDEPG 

FDFPND 

THEM 

FRINT 

fejmht 

FBMLT 

FBTINT 

FOOTPC 

f  PLNK 

GWPLME 

HTOS 

INPUTP 

JAMND 

MAGFLD 

MANPTS 

NINT 

OMF  GA 

PEDEP 

PHFNOM 

PHOTOR 

P I  ONF 

PROPEN 

PTH I Nr 

RINTLR 

SATC 

5ATLNX 

SEPA 

TLTFLD 

TORBIT 

CONWES 

DBSTAM 

GWPLME 

IONOSI 

PHLNOM 


Table  4.  Routines  usinq  a  specific 
common  block. 


COMMON  BLOCK 

USED  IN 

COMMON  BLOCK 

USED  IN 

COMMON  BLOCK 

Al  TOON 

CONNER 

RE  TAB 

ANTNO 

RTUBE 

ATHOSI 

BETAS 

ATNPS'i 

RLKOTA 

UHEG 

CON^F  S 

TFHPZH 

ANTNO 

ATmh$f 

BURST 

AM  HP 

AN  T  NO 

BETAS 

RF  T  AS 

CHEPf) 

FRTINT 

CNf)F.oT 

F  Pi  Nk 

C*F  EOT 

ninT 

DRHHT 

DBSTah 

ATMUSN 

CONNF  R 

DERMIS 

amhah.r 

of  oEp 

ATPHSU 

uhsthi 

CHFMn 

t  Of  PH 

CHEHOfJ 

EOfPNO 

CHEMHR 

UPEn 

ENf  F 

FNfF 

FOAMH 

E.uxfjt 

SPECOP 

FREN 

ERIN! 

atmost 

CHE 

FRITmI) 
f  i 

CHEMOO 

CMDEnT 

F  RHMT 

OTNEP 

FHBl  I 

OTNECJ 

frmlt 

E  N?V  I B 
FQAMM 

ehitnt 

footpc 

CHE*»AN 

I  N  I  T  AL 

FPLNk 

SPECOP 

CiNPLHF 

SPECOQ 

"ANPTS 

PCI  our: 

CHE“P 

pfoEp 

atmoup 

CONNF  S 

PHENHM 

ATMOSI 

phuTiio 

ATMOStl 

phiiTiim 

pthnf 

R  A  DC  KIT 

bfoit 

CHE  Hfi 

STRIP 

CHDFOT 

C  HE  MF  R 

I 

USED  IN 


DFOE  P 
FBINT 
FRTINT 
MANPT  S 


COnwfs 

BETAS 

OHHHT 

OBMl  T 

OUSTmT 

FDfcPH 

fcO£PG 

f  Of  PND 

FU>EN 

F IIXF  I  T 

faint 

FRMMI 

Fmmi  i 

FRTINT 

I DF  Fl  T 

INPUT 

INP'ITP 

HANPTS 

PFOEP 

PHENOM 

PTllNF 


fHF.Mt  F 
ChnImn 


CMf.HMP 

FREN 

F  R 1 N  T 

FRTINT 

FlKlTPC 

RATC 


FRFN 

FBINT 

FOtiTPC 
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Table  4  (Continued) 


COMMON  BLOCK 

USED  IN 

COMMON  BLOCK 

USED  IN 

COMMON  BLOCK 

USED  IN 

CHE*M 

CONMES 

CONST 

CONNES 

IONdSU 

8LK0TA 

a*hab$ 

JAMND 

RATE 

ANTNn 

MAGDIS 

ARRLIH 

magft  t 

ATMQSF 

MANPTS 

CMMOVL 

CMfcMD 

ATMOS! 

NFS*mF 

flofn 

ATMUSU 

nInT 

ENE  F 

AZELR 

omega 

balfit 

ORBIT 

CINPUT 

CflNWES 

bhpsk 

OZONE 

ANT  NO 

BETAS 

PfLOuO 

ARRL  I M 

BFIFLf) 

PEDFP 

BETAS 

RlKf>TA 

PGMOtjP 

CHEMO 

BTRACE 

PHENOM 

ChE“OQ 

COMB 

PIUNF 

obstam 

coordv 

PMASS 

OEDEP 

COOPDX 

PROPEN 

OSTINT 

COROTF 

PSDM 

ELOEN 

DBMHT 

PTHINT 

enef 

obmlt 

dcosm 

FBEN 

ubstam 

rinter 

fbint 

debris 

sanc an 

FOUTPC 

DELPSK 

satcir 

FPLNK 

P>PSK 

SATLNK 

IOEFLT 

DS1 INT 

SF  PA 

INPUT 

EOEPH 

SLOBF 

Inputp 

EOEPG 

SOI  CYC 

manpts 

ENF.F 

SOL  ORB 

PhENOM 

EOCRUS 

SOlZEN 

PROPEN 

ERF 

SORB IT 

pthint 

F  RFC 

STRIJ 

satcir 

F  AC 

TEMPZH 

satlnk 

FBFN 

tltfld 

tqrbit 

FBINT 

TORBIT 

FBMHI 

UVWD 

CIROES 

CONHES 

FBMMT 

VOLUME 

IDEFL! 

FBMLI 

MATER 

input 

FBMLT 

Z  T  TOUT 

Inputp 

FOUTPC 

satc 

FPCIR 

DBRFG 

CONWFS 

SLINK 

fpint 

ARRl IM 

FPLNK 

BFTAS 

Cl  INK 

antnd 

FPPLOT 

DflMHT 

FPLNK 

FS«M 

ORMLT 

J  AMNO 

GEUCIJR 

DBS! AM 

satc 

GKPLME 

DEDEP 

satlnk 

HD  I  NT 

EDFPB 

slink 

hToS 

FDFPG 

ideflt 

FBINT 

INPUT 

MANPTS 

INPUTP 

PHENOM 

I ONUS I 

47 


Table  4  (Continued) 


COMMON  BLOCK  USED  IN 


COMMON  BLOCK  USED  IN 


COMMON  BLOCK  USED  IN 


DCREG 


DEPDAT 


DEVICE 


DIGMUO 


DTUBE 


CON*ES 
OSTInT 
DUSTmI 
PCI  OUO 
PGROUP 
PMfcNOH 
SOIL 

CON«E8 

0LAOTA 

PCMEM 

PIONF 

conmes 

BETAS 

OBHMT 

DBMLT 

EHEPNO 

ELOEN 

fbint 

fbtint 

INPUT 

INPUTP 

MANPTS 

PEOEP 

PMfcNOM 

PIONF 

MlftD 

rogd 

*OG1 

«ONl 

mOXI 


BOPS* 

CPSK 

HELPS* 

dps* 

OOPSK 

FPCIW 

FSKWH 

fsmhy 
nfskmf 
SATCIR 
OSf PMO 


OFRRlS 

MANPTS 

PEOEP 


EDTOVL 


FRRFG 


FOSRAT 


FPCIRC 


fpjntc 


HFaCOM 


CNFEOT 

PROPEN 

CONNE8 

ANTNO 

AT*n3F 

BETAS 

DBNMI 

ORNLT 

OEOEP 

03TINT 

EOFPH 

EOEPG 

FBEN 

FBINT 

FB*HT 

FBMLT 

FBTlNT 

FOOT PC 

fplnk 

GWPLME 

NINT 

pclhud 

PMENOM 

PHOTyR 

RADOtl 

START 

STRIP 

CHFMO 

CME^OO 

CHE  mHR 

DTNfP 

OTnFR 

FOAMfl 

PHOTOR 

FPCIR 

FPUNk 

FOOTPC 

FPINT 

FPPLflT 

CONRES 
ATMOSF 
f  NEF 


HE  AVF  1 


INPTC« 


JONOUP 


JAMMER 


LENGTH 


MAGLNlf 


CONMES 

ATMOSF 

HFAVE 

LAGRAN 

pmeave 

INPUT 

INPIJTP 

IONOSI 

IONOSU 

CONMf.S 

IOFFLT 

INPUT 

INPUTP 

JAHNO 

SATC 

balfit 

bTRAcF 

DBMHT 

EQCR(»S 

F0FN 

FBMMT 

FOOTPC 

HOINT 

MAGDIS 

PHENOM 

CONMES 

ANTNO 

hfielo 

DEBRIS 

OFOFP 

EOFPB 

FBFN 

FBINT 

FBTlNT 

FOOTPC 

FPLNk 

IONUSI 

IONOSU 

M  A(fF  |  T 

HAGFID 

MANPTS 

OMEGA 

PHENOM 

PROPEN 

pthjnt 

RINTEH 


48 


Table  4  (Continued) 
COMMON  BLOCK  USED  IN 


MODULI  CONHEfl 

FPCIH 
ioeflt 

INPUT 

INPUTP 

S4TCIR 

OPTION  CONNES 

ANTNO 
8LKOTA 
C  ARO I N 
CHE"D 
CHEMhR 
CNOEPT 
CHFEDT 
OPMHT 
06STAM 
OFRRIS 
DFOEP 
DSTlNT 
FLOEn 
fnff 

FRFN 

frint 

FBTInT 

FOOTPC 

FPI*  T 

fplnk 

FPPLOT 

lOFFtT 

INPUT 

INPUTP 

JAMNO 

RANPT3 

NlNT 

PCHER 

PHENOM 

PIONF 

PROPER 

PTMINT 

8ATC 

3ATCIR 

Satlnk 

SORB  IT 


COMMON  BLOCK 

j  USED  IN 

ORIGIN 

CONNES 

DATE 

orstar 

GNPL*F 

IOEFlT 

INPUT 

INPUTP 

ORHIT 

PHFNOR 

SOLCVC 

OUT 

CHE  HD 

CHDFDT 

DTNFQ 

F  QAMB 

PATH 

CONNES 

ideflt 

INPUT 

INPUTP 

UPPIT 

SATC 

SOPBIT 

TOPRIT 

PMEN 

CONNER 

DRHHT 

DBMLT 

DBSTa" 

DU3THI 

F0MH I 

FRHHT 

FRMLI 

FBMLT 

GnPLHF 

PHENflH 

padout 

PROP 

CONMES 

ANTNO 

DSTlNT 

F8INT 

fplnk 

JARNO 

RANPT3 

PRfjPEN 

pthint 

SATC 

SATLNK 

START 

strii 

COMMON  BLOCK 

USED  IN 

RAY 

APPL  t  H 
GNPLHF 
PSDH 
KVNAV 

uvwo 

RFGINT 

CONNES 

fbint 

RANPTS 

PROPER 

PTMINT 

start 

RGFOR 

RINTEP 

TOROID 

TORE 

SIGCOM 

SATC  IR 
SATLNK 

SOURCF 

CHFEDT 

flden 

FNEF 

SPEC 

CHFHHR 

FOLAIR 

FOLHTl 

SPECEF 

chehef 

CRH ION 

CHFEDT 

enef 

hpcheh 

pchfr 

PIONF 

SPECO 

CHFHO 

ENEF 

SPECN 

CONNES 

INPUT 

SPCTRM 

NOG1 

NON! 

NOXI 

SPEEDS 

PSDH 

UVNAV 

Table  4  (Continued) 


COMMON  BLOCK 

USED  IN 

COMMON  BLOCK 

USED  IN 

ssvs 

CONMES 

NEDEPO 

CONNf s 

ANTND 

blkhta 

FPL** 

CMXSPC 

ioeflt 

CONSPC 

* 

INPUT 

DERMIS 

INPUTP 

LFKSPC 

JAMNf) 

SATC 

SATLNK 

NOG 

CONTES 

WOOD 

SUNK 

wOGP 

wnr.i 

STRIIN 

FPINT 

STWII 

WON 

con wts 
wOwD 

TIME 

CONNfS 

mONP 

ATHUSI 

ATmDSL 

WON  t 

DBSTAH 

10N0ST 

IONOSU 

JULIAN 

uzone 

wOX 

CUNWf  S 

wOXC 

wnxP 

WflXl 

PHENfjM 

SATC 

WRATE 

CMDEDT 

SOLCVC 

DRATE 

SOLOBB 

OTNEP 

SOlZEN 

OTnEO 

TEMPZH 

INITal 

WATER 

SPfcCDP 

ZTTHijT 

SPECPC 

trace 

HTPACE 

WPATEH 

orate 

FBPHT 

SPECOP 

TSTEP 

HPCM£M 

NSTflR 

GWPLME 

PCHFM 

PSO* 

piunf 

uvwn 

T3VS 

CPNNfS 

ZHTEMP 

ATMOSI 

ANTNO 

ATwnsi) 

fplnk 

TFWPZH 

IDEFtT 

INPUT 

INPUTP 

JAHNl) 

SATC 

SATLNK 

SUNK 
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SECTION  3 

ROUTINE  DESCRIPTIONS 


A  description  of  inputs  and  outputs  transferred  through  the  call 
statement  and  a  list  cf  inputs  and  outputs  transferred  through  labeled  common 
blocks  are  given  in  consent  statements  at  the  beginning  of  each  routine.  A 
description  of  quantities  in  labeled  common  blocks  is  given  in  Table  5. 


Most  of  the 
of  routines  developed 
descriptions  routines 
been  modified  for  use 


routines  used  in  the  KTSCOM  code  are  modified  versions 
for  the  WHIT!  code.  However,  as  indicated  in  the  following 
developed  for  the  SATL,  ROSCOli,  and  KOI:  codes  have  also 
in  the  KF.SCOM  code. 


51 


I 

Table  5.  Quantities  stored  in  labeled 
common  blocks. 


COW-ION  BLOCK 

QUANTITY 

DESCRIPTION 

ALTODN  (Atmosphere 
Coninon) 

ALTKM 

Altitudes  used  for  atmospheric  calculations 
in  routine  ATMOSF 

ZHT 

Altitudes  used  for  temperature  and  ozone 
calculations  in  routine  ATMOSI 

ANTMP  (Noise  Common) 

BEAN 

Antenna  beam  efficiency 

DEPTH 

Effective  distance  into  the  fireball  pro¬ 
ducing  thermal  radiation  (km) 

FREQN 

Frequency  (MHz) 

KEYN 

Code  describing  location  of  fireball  in 
antenna  beam 

KPLME 

Code  describing  fireball  z-axis  curvature 

N 

Fireball  index 

NCFB 

Index  of  fireball  closest  to  receiver  in 
beam  center 

RBO 

Antenna  beam  width  (radians) 

RC 

Side  lobe  gain  relative  to  RGO 

RGO 

Main  beam  gain 

RXKO 

Side  lobe  specification 

SRFBE 

SRF8X 

Slant  ranges  to  entrance  and  exit  inter¬ 
sections  with  fireball 

TANT 

Normalized  antenna  temperature 

TASUM 

Antenna  temperature  summation  for  fireball 
loop  (K) 

TIME 

Calculation  time  (s) 

UV 

Transformation  matrix  from  tangent  plane  to 
earth  centered  coordinates 

VBEAM 

Antenna  pointing  direction  vector  (km) 

VR 

Receiver  location  vector  (km) 

VSPA 

Sight  path  vector  to  fireball  (km) 

VZZ 

Fireball  reference  vector  (km) 

ATMOSN  (Atmosphere 
Coninon) 

AC02  \ 

; 

Concentration  of  CO2*  H,  ...  0^  (cm  ) 

A03  > 

Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

ATMOST  (Atmosphere 

1 

Common) 

nil 

HI 

Concentration  of  CO2,  H,  ...  O^  (cm  ) 

Density  scale  height  (km) 

P 

Pressure  (dynes  cm  ) 

RHO 

Density  (gm  cm  ) 

TEM 

Temperature  (K) 

ATMOUP  (Atmosphere 

AC02T  \ 

Common) 

AH02T  I 

: 

Concentration  of  C02»  HO2.  ...  03  (cm"J) 

A03T  / 

FDAY 

Day  profile  of  mean  molecular  weight 

FHR120 

Density  scale  height  term 

GAMT 

Coefficient  of  M.  in  equation  for  Y^ 

Density  scale  height 

HRHO 

HR0110 

Density  scale  height  at  110  km 

_„PP 

Pressure  (dynes  cm"2) 

RHO 

Density  (gm  cm"'*) 

SA 

(TIF-TZ)/TIF 

SB 

Approximate  derivative  of  density  scale 
height  at  110  km 

SNIZ 

Species  densities  at  120  km  (cm’**) 

TAU 

Temperature  gradient  at  120  km 

TIF 

Exospheric  temperature  (K) 

TT 

Temperature  (K) 

TZ 

Temperature  at  120  km  (K) 

BEDIT  (D-region  Edit 

EDP 

Energy  deposited  by  prompt  radiation 
(ergs  gm-1) 

Common) 

EXPD 

Fraction  of  O3  remaining  after  thermal 
radiation 

INDEX 

Indices  of  bursts  affecting  D-region  neutral 
species 

NPBC 

Index  of  last  burst  in  INDEX 
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COMMON  BLOCK 


BETAB  (Beta  Sheath 
Common) 


BREG  (Burst  Common) 


QUANTITY 

UtiLKIM  iUrt 

DBBETA 

Absorption  through  beta  sheath  (dB) 

ETAS 

Beta  sheath  reflection  coefficient 

FSCA 

Scaled  frequency  used  for  reflection  coef¬ 
ficient  calculation 

HPBETA 

Altitude  of  field  point  (km) 

MB ETA  j 

0  Beta  sheath  absorption 

1  Absorption  and  reflection  coefficient 

NBETA 

Fireball  index 

RHOBET 

Fireball  air  density  (gm  cm  ) 

AENE 

.3 

Apogee  Ng  (cm  ) 

AHF 

Apogee  reference  altitude  (km) 

AHMAXF 

Apogee  maximum  altitude  (km) 

AHMINF 

Apogee  minimum  altitude  (km) 

AHSNE 

Apogee  Ne  scale  height  along  longitudinal 
axis  (cm"3  km-1) 

AHSNK 

Apogee  density  scale  height  along  longitudinal 
axis  (km) 

AHTF 

Apogee  temperature  scale  height  along  longi¬ 
tudinal  axis  (km) 

ALPHUV 

UV  parameter 

ari:of 

Apogee  mass  density  at  fireball  bottom 
(gm  cm-3) 

ARLF 

Apogee  longitudinal  dimension  (km) 

ARTF 

Apogee  transverse  dimension  (km) 

ATF 

Apogee  temperature  (k) 

ATILTF 

Apogee  tilt  angle  (radians) 

IBI 

Apogee  earth  centered  vector  location  (km) 

AXL 

Apogee  longitudinal  positions  (km) 

DBV3 

Debris  volume  when  fireball  is  3000  K  (knr) 

DECANB 

Magnetic  declination  angle  at  burst  ooint 
(radians) 

DIPANB 

Magnetic  dip  angle  at  burst  point  (radia.is) 

ENERGY 

Number  of  photons  in  4.135  to  7.255  eV  band 
(1 . E-18  photons) 

EQRHOF 

Mass  density  at  bottom  of  fireball  at  t 
(gm  cm-3) 

■  > 

Burst  point  density  scaling  parameter 

IO'  i 

Fraction  of  hydroenergy  converted  to  UV 

HB 

J _ 

Burst  altitude  (km) 

1 

j 

I 


Table  b  (Continued; 


COMMON  BLOCK  QUANTITY 


BREG  (Continued) 


HB0T2 


HSTOP 


HS2000 


IBURST 

IDAYB 

IFBRM 


MRGID 


RDZERO 

REFT 

RHOB 

RHOFO 

RHUV 

RHZERO 

RH2000 

RH3000 

RIONL 

RLOSC 

RM 

RUZERO 

TAPOGE 


DESCRIPTION 


Altitude  of  fireball  bottom  when  fireball 
temperature  Is  2000  K  (km) 

Density  scale  height  at  burst  point  (km) 

Average  scale  height  between  bottom  and  center 
of  Initial  fireball  (km) 

Stopping  altitude  for  ballistic  fall  of  Ioniza¬ 
tion  (km) 

Density  scale  height  at  bottom  of  fireball 
when  fireball  temperature  is  2000  K  (km) 

Burst  index 

Time  of  day  at  burst  time  (day  or  night) 
Fireball  merge  flag 

0  if  fireball  bottom  goes  above  HSTOP,  1  if 
below 

Apogee  fireball  shape 

1  if  burst  causes  heave 

1  If  low  altitude  burst,  0  if  high  altitude 
burst 

Merge  index^"'"'>>^ 

Pressure  at  burst  potM^Mynes  cm"*) 

Charge  exchange  characterlstH^Cadius  (km) 
Downward  UV ‘radius  (km) 

Initial  downward  radius  (km) 

Initial  vortex  surface  temperture  (k) 

i  -3 

Burst  point  mass  density  (gm  cm  ) 

Initial  bottom  ma^s  density  (gm  cm"'5) 

Horizontal  UV  radios  (km) 

Initial  horizontal  Vadius  (km) 

Fireball  density  when  fireball  temperature  is 
2000  K  (gm  cm-3)  \ 

Fireball  density  when  fireball  temperature  Is 
3000  K  (gm  cm"3)  \ 

Ion  leak  characteristic  Vidius  (km) 

Loss  cone  characteristic  radius  (km) 

Magnetic  containment  radius^  (km) 

Initial  upward  radius  (km) 

Apogee  time  (s)  \ 

Burst  time  (s) 


Table  5  (Continued) 


COMMON  BLOCK  QUANTITY 


BREG  (Continued) 


BTUBE  (Beta  Tube 
Common) 


BURST  (Burst  Input 
Common) 


VMAXD 

VMFBO 


VMIND 

W 

WCHEX 

WILD 

WLCD 

WUV 

XEMAX 


HBETA1 

HBETA2 

KBS 

ICONJ 


DESCRIPTION 


Temperature  at  burst  point  (K) 

Magnetic  containment  time  (s) 

Initial  fireball  temperature  (K) 

Magnetic  freezing  time  (s) 

Time  for  all  parcels  to  fall  (s) 

Merge  time  (s) 

Time  of  completion  of  torus  formation  (s) 

Time  of  end  of  first  expansion  phase  (s) 

Time  at  end  of  second  expansion  phase  (s) 

Time  when  fireball  temperature  is  3000  K  (s) 

Earth  centered  vector  to  burst  location  (km) 

Heavy  particle  velocity  (km  s’1) 

Earth  centered  vector  location  of  dipole  for 
burst  (km) 

Earth  centered  vector  location  of  heavy 
particle  reference  location 

Heavy  particle  velocity  (km  s’1) 

Burst  yield  (MT) 

Charge  exchange  yield  (MT) 

Ion  leak  yield  (MT) 

Loss  cone  downward  yield  (MT) 

UV  yield  (MT) 

Position  of  parcel  in  apogee  fireball  that 
stops  at  equator  (km) 

Position  of  parcel  in  apogee  fireball  that 
goes  above  HSTOP  (km) 

Wake  parameter 

Wake  parameter 

Altitude  of  lower  beta  tube  intersection  (km) 
Altitude  of  upper  beta  tube  intersection  (km) 
Beta  sheath  flag 

Burst  location  flag 
Device  index 
Surface  index 
Number  of  bursts 
Burst  time  (s) 
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Table  5 


COMMON  BLOCK 


BURST  (Continued) 

CHEMAN  (F-region  Chem¬ 
istry  Common) 


CHEMB  (Fireball  Chem¬ 
istry  Common) 


QUANTITY 


ALPHA 1 
ALPHA2 
AMOLP 
BETA1 


EXPAND 

Q 

RH021 
XK( I )M 


FEQN 

FUVTC 


Q 

RAL 

RHO 

RHOALQ 

RHOEQ 

RHOEQM 

RHOEQN 

RHOHB 

RHO  I 

RHOUQ 

RU 


DESCRIPTION 


Earth  centered  vector  to  burst  location  (km) 

O 

Atomic  Ion  recombination  rate  (cm  /s) 

3 

Molecular  Ion  recombination  rate  (cm  /s) 

.3 

Molecular  Ion  concentration  (cm  ) 

Effective  loss  rate  due  to  volume  expansion 

<*-’>  3 
Energy  returned  from  chemistry  (eV/cm  ) 

Time  Interval  (s) 

Time  average  volume  expansion  ratio 

-3  -1 

Total  Ion  pair  production  rate  (cm  s  ) 

Ratio  of  density  at  time  t2  to  density  at 
time  t-j 

Mean  rate  coefficient  for  reaction  I 

Entrainment  fraction 
Entrainment  fraction  at  2000  K 
Fraction  of  hydroenergy  converted  to  UV 
Conjugate  region  flag 
Footprint  calculation  flag 

-3  -1 

Ion  pair  production  rate  (cm  s  ) 

.3 

Aluminium  density  (gm  cm  ) 

.3 

Air  density  (gm  cm  ) 

_3 

Aluminium  density  at  time  TEQM  (gm  cm  ) 

_3 

Air  density  at  TEQ  (gm  cm  ) 

Air  density  at  TEQM  (gm  cm"3) 

Air  density  at  TEQN  (gm  cm"3) 

Burst  point  air  density  (gm  cm  ) 

.3 

Initial  air  density  (gm  cm  ) 

-3 

Uranium  density  at  time  TEQM  (gm  cm  ) 

Uranium  density 

Magnetic  equilibrium  time  (s) 

Time  after  burst  when  temperature  Is  TFEM  (s) 

Fireball  temperature  (K) 

Non-equilibrium  metal  ion  temperature  (K) 

Non-equilibrium  neutral  air  species  tempera¬ 
ture  (K) 

Initial  temperature  (K) 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

CHEMB  (Continued) 

TOB 

Burst  time  (s) 

XI 

Scaled  distance  from  fireball  bottom 

CHEMFB  (Fireball  Chem- 

EfIEFP 

Electron  density  (cm  ) 

is try  Common) 

FPXL 

Distance  from  earth's  surface  measured  along 
geomagnetic  field  line  (km) 

HSNEFP 

Electron  density  scale  height  (km) 

HSTFFP 

Temperature  scale  height  (km) 

JNFBFP 

Fireball  index 

TFBFP 

Temperature  (K) 

CHEMR 

AR,  BR, 

CR 

Reaction  rate  coefficient  parameters 

CHMOVL  (F- region  Chem- 

ENEC 

Electron  density  (cm  ) 

is try  Common) 

ENPC 

Positive  ion  density  (cm  ) 

HPC 

Altitude  (km) 

IALT 

Altitude  index 

IONC 

Ionization  flag 

NPBX 

Number  of  prior  bursts 

VEIC 

electron- ion  collision  frequency  (s'1) 

VEMC 

Electron-neutral  collision  frequency  (s-^) 

VHEAVE 

Heave  velocity  (km/s) 

Ion-neutral  collision  frequency  (s-1) 

VIMC 

VPOC 

Earth  centered  vector  to  field  point  (kn) 

CINPUT  (Input  Common) 

IOAY 

Time  of  day  (day  or  night) 

JNDB 

Debris  region  index 

JNDC 

Dust  cloud  index 

JNFB 

Fireball  index 

JTIME 

Calculation  time  index 

LNCNT 

Output  line  count 

NPROB 

Problem  number 

NT 

Number  of  calculation  times 

PROB 

Problem  heading 

T 

Calculation  times  (s) 

i 

TESTIP 

Input  test  flag 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

CIRDES  (Circuit  Common) 

■■ 

System  indexes  for  downlink 

System  indexes  for  satellite-to-satellite  link 

mm 

Ground  terminal  index 

ESI 

Satellite  terminal  index 

ITI3 

Satellite  terminal  index 

ITI4 

Ground  terminal  index 

IUL 

System  indexes  for  uplink 

NC 

Number  of  circuits 

CLINK  (Link  Common) 

INDEX 

Circuit  index 

Link  index 

1 

EKi  ■■ 

System  indexes  for  specified  link 

LFLAG 

Propagation  calculation  flag 

LINK 

Circuit  index 

LI,  L2, 

L3,  L4 

Terminal  indexes 

NFLAG 

Noise  calculation  flag 

CONST  (Constant  Common) 

G 

Acceleration  of  gravity  (km  s  ) 

PI 

TT 

RADIAN 

tt/180 

RE 

Earth's  radius 

RHOG 

Air  density  at  earth's  surface  (gm  cm  ) 

XNHAT 

Earth  centered  unit  vector  to  geographic 
pole  (km) 

DBREG  (Debris  Region 

DB 

Beta  sheath  thickness  (cm) 

Common) 

DISTD 

Debris  distribution  parameter 

FEI 

Fraction  of  beta  energy  that  escapes  debris 
region 

FEO 

Fraction  of  escaping  energy  deposited  in 
beta  sheath 

HBETA 

Beta  source  altitude  (km) 

HD 

Debris  altitude  (km) 

HMAXD 

Maximum  debris  truncation  altitude  (km) 

HMIND 

Minimum  debris  truncation  altitude  (km) 

HSND 

Debris  scale  height  (km) 

IXDBB(L) 

Burst  index  associated  with  L th  dust  cloud 

Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

DBREG  (Continued) 

IXDBFB 

Fireball  index  associated  with  L th  debris 
region 

KBET 

Beta  sheath  flag 

Debris  region  geometry  flag 

NDB 

Number  of  debris  regions 

RLBETA 

North-south  radius  of  beta  source  plane  (km) 

RLD 

Longitudinal  debris  radius  (km) 

RTBETA 

East-west  radius  of  beta  source  plane  (km) 

RTD 

Transverse  debris  radius  (km) 

TILTD 

Debris  tilt  angle  (radians) 

VDBO 

Earth  centered  vector  to  debris  reference 
point  (km) 

VOLD 

Debris  volume  (km3) 

Beta  source  fission  yield  (ITT) 

Gamma  source  fission  yield  (MT) 

Larmor  radius  outside  debris  region  (cm) 

DCREG  (Dust  Cloud 

ABAR 

Particle  size  group  boundary  values  (cm) 

Common) 

DIECR 

Real  part  of  index  of  refraction 

DIECI 

Imaginary  part  of  index  of  refraction  | 

FNP 

Fraction  of  particles  In  each  group 

IMAXP 

Index  of  maximum  particle  group  lofted 

IXDCFB(L) 

Fireball  index  associated  with  L th  dust  cloud 

JGROUP 

Group  Index 

NDC 

Number  of  dust  clouds 

NPG 

Number  of  particle  groups 

RHODP 

Particle  mass  density  (gm  cm'3) 

SIGB 

Average  backscatter  cross  section  for  each 
group  (m2) 

SIGE 

Average  extinction  cross  section  for  each 
group  (m2) 

ViP 

Total  number  of  particles  in  all  groups 

\MLOD 

Dust  loading  factor  (MT/MT) 

DEPDAT  (Prompt  Energy 

ERGU(L) 

Photon  energy  for  UV  group  L  (erg) 

Common) 

PREFF(L.I) 

Pressure  efficiency  for  energy  group  L  and 
species  I 

SIGU(L,I ) 

UV  absorption  cross  section  for  energy  group 

L  and  species  I  (cm2) 
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Table  5  (Continued) 


COMMON  SLOCK 

QUANTITY 

DESCRIPTION 

DEPDAT  (Continued) 

THRESH(I) 

Effective  threshold  for  absorption  of  UV 
photons  by  species  I  (erg) 

DEVICE  (Device  Common) 

FEDG 

Fraction  of  weapon  energy  radiated  as  prompt 
gammas 

Fraction  of  weapon  energy  in  hydro 

E  .  ■ 

Fraction  of  weapon  energy  in  neutrons 

FEDX 

Fraction  of  weapon  energy  radiated  as  X-rays 

FF 

Fission  fraction 

FMAL 

Fraction  of  weapon  mass  that  is  aluminum 

FMUR 

Fraction  of  weapon  mass  that  is  uranium 

FTHERM 

Fraction  of  fissionable  material  that  is 
fissionable  by  thermal  neutrons 

WI 

Weapon  yield  (MT) 

WMASS 

Weapon  mass  (g) 

DIGMOD  (Signal  Proces- 

A 

Ratio  of  bit  period  to  decorrelation  time 

sing  Common) 

mm 

Scintillation  Index  function 

1  -  1/MARY 

D I SPAR 

Dispersion  parameter 

Bit  energy-to-noise  ratio 

Chip  energy-to-noise  rat’o 

GAMMA 

Normalized  energy-to-noise  ratio 

GPR 

Function  of  number  of  modulation  states 

MARY 

Number  of  modulation  states 

NINFOB 

SUMMARY) 

S4 

Scintillation  Index 

TAU 

Decorrelation  time  (s) 

WT 

Product  of  chip  period  and  frequency  separation 

DTUBE  (Debris  Kinetic 
Energy  Common) 

LDKE 

Radiation  source  flag 

EDTOVL  (F-reglon  Burst 

HFTEST 

Test  altitude  (km) 

Edit  Common) 

JNFEDT 

Last  burst  index 

JPBEDT 

Index  of  last  burst  producing  F-reglon 

Ionization 

VPO 

Earth  centered  vector  to  test  location  (km) 
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Table  5  (Continued) 


COMMON  iLOCK 

QUANTITY 

DESCRIPTION 

FBREG 

AFB,  AFC, 
BFB,  BFC 

Coefficients  for  parcel  location 

BVAL 

Oval  of  Cassini  parameter 

CFB,  CFC 

Coefficients  for  parcel  location 

ENE1DF 

Diffusion  ionization  parameter  (cm"3) 

HF 

Fireball  reference  altitude  (km) 

HLDF 

Hydrostatic  scale  length  (km) 

HMAXF 

Maximum  altitude  (km) 

HMINF 

Minimum  altitude  in  burst  region 

HMINFC 

Minimum  altitude  in  conjugate  region  (km) 

HSNF 

Density  scale  height  along  longitudinal  axis 
(km) 

KINDF 

Shape  indicator 

KPB 

Burst  region  parcel  indicator  (0  =  none) 

KPC 

Conjugate  region  parcel  indicator 

NFB 

Number  of  fireballs 

POWER 

Thermal  power  radiated  in  specified  bands 
(10"'°  photons  s"') 

Radial  expansion  velocity  (km  s-1) 

ROOT 

RHOF 

Bottom  density  (gm  citi  ) 

RLF 

Longitudinal  radius  (km) 

RTF 

Transverse  radius  (km) 

RLFC 

Conjugate  longitudinal  radius  (km) 

RGEOM 

Radial  geometry  factor  used  in  FBEN 

RMAXDF 

Radius  from  earth  center  to  maximum  altitude 
on  field  line  (km) 

RMAXF 

Distance  from  earth  center  to  field  line  at 
equator  (km) 

RPDF 

Radius  from  earth's  center  to  diffusion  ioni¬ 
zation  peak  (km) 

RTREF 

Radius  of  reference  Oval  of  Cassini  (km) 

RVL 

LAFBP  -  vortex  longitudinal  radius  (km) 

RVT 

LAFBP  -  vortex  transverse  radius  (km) 

SPLEX 

Power  law  exponent 

SSCAI 

Inner  scale  length  (km) 

SSCAO 

Outer  scale  length  (km) 

TOEB 

Debris  temperature  (!<) 

TF 

Temperature  (K) 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

FBREG  (Continued) 

TILTF 

Tilt  angle  (radians) 

TVORT 

VARNEN 

Vortex  surface  temperature  (K) 

Electron  density  variance  (cm"^ 

VFBO 

t 

Earth  centered  vector  defining  central  field 
line 

VFBVEL 

Earth  centered  vector  defining  velocity  of 
fireball  (km  s-*) 

VOV 

Vortex  volume  (km  ) 

XLLDF 

Hydrostatic  scale  length  along  field 
line  in  diffusion  fireball  (km) 

XLMDF 

XLPDF 

Netural  density  vertical  scale  height 
in  diffusion  fireball  (km) 

Distance  from  earth's  surface  to  peak  of 
diffusion  ionization  (km) 

FDSRAT  (Fireball  Photo¬ 
dissociation  Common) 

FDIS3  ) 

: 

Photodissociation  rate  (s) 

FDis9  / 

FDT1  ) 

: 

Photodetachment  rate  (s) 

FDT8  f 

FPCIRC  (Footprint  Cir¬ 
cuit  Common) 

— 

Noise  bandwidth  (MHz) 

Ambient  noise  density  (w  Hz’1) 

Disturbed  noise  density  (W  Hz’1) 

ENTOT 

Electron  density  integral  (cm  ) 

FREQGH 

Frequency  (GHz) 

Frequency  selective  bandwidth  (Hz) 

I  MOD 

Modulation  index 

KALSN 

Signal-to-noise  option 

KFMFP 

FM  modulation  flag 

PERFP 

Probability  of  bit  or  message  error  (s’1) 

SIGFS 

Ambient  signal  strength  (w) 

SIGND 

Disturbed  signal  strength  (W) 

SPHIL 

rms  phase  fluctuation  (radians) 

SPSIL 

rms  angle  of  arrival  fluctuation  (radians) 

S4L 

Scintillation  strength  index 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

FPCIRC  (Continued) 

TAUZL 

Decorrelation  time  (s) 

FPINTC  (Footprint  In¬ 
terpolation  Conmon) 

HR 

KFMFP 

Receiver  altitude  (km) 

FM  modulation  flag 

i 

Nil 

Scan  number  ending  first  set  of  contours 

NI2 

Scan  number  starting  second  set  of  contours 

NN 

Total  number  of  scans 

NPK 

Number  of  paths  for  specified  scan 

NPKM 

Index  of  path  with  maximum  propagation  effects 

QT 

Propagation  quantities 

QTI 

Interpolated  propagation  quantities 

SKD 

Distance  from  origin  to  paths 

TLAMAX  1 
TLAMIN  j 

Minimum  and  maximum  latitudes  for  scan 
(deg) 

TLATFL 

Latitude  of  central  field  line  (deg) 

TLATRD 

Latitude  of  calculation  paths  for  scan  (deg) 

TLATRI 

Latitude  of  paths  for  interpolated  propagation 
quantities  (deg) 

TLOMAX  ) 
TLOMIN  j 

TLONFL 

Minimum  and  maximum  longitudes  for  scan 
(deg) 

Longitude  of  central  field  line  (deg) 

TLONRD 

Longitude  of  calculation  paths  for  scan  (deg) 

TLONRI 

Longitude  of  paths  for  interpolated  propaga¬ 
tion  quantities  (deg) 

VPO 

Earth  centered  vector  to  field  point  (km) 

VR 

Earth  centered  vector  to  receiver  terminal  (km) 

VXP 

Earth  centered  path  direction  vector  (km) 

HEACOM  (Heave  Interpo¬ 
lation  Common) 

HPINT 

NINT 

Altitude  (km) 

Number  of  interpolated  points 

RHOINT 

Density  (gm  cm  ) 

TINT 

Time  after  burst  (s) 

HEAVE1  (Heave  Common) 

A  ) 

( 

Calculation  quantities  used  in  heave  model 

_ 1 

ZRSS  ' 

] 

l 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

INPTCM  (Input  Common) 

BCORDK 

1 

Transfer  of  internal  input  quantities  between 

:  ( 

routines  INPUT  and  INPUTP 

VROO  * 

IONOUP  (Ionosphere 
Common ) 

BA  1 

BB  ) 

Booker  interpolation  parameters 

IONOUP  (Continued) 

BNER 

_3 

Electron  density  at  reference  point  (cm  ) 

BZ 

Transition  altitudes  (km) 

UR 

Reference  altitude  (km) 

HI 

Upper  altitude  of  production  rate  fits  (km) 

XI  ON 

Ionosphere  fit  quantities 

JAMMER  (Jarnner  Common) 

DJAM 

Jammer  noise  density  (W  Hz”1) 

IJS 

Index  of  satellite  terminal 

IJT 

Index  of  ground  terminal 

NJ 

Number  of  jammers 

VJO 

Earth  centered  vector  to  jammer  (km) 

LENGTH  (Field  Line 

COSLO 

Colatitude  of  earth  Intercept  (cosine) 

Geometry  Common) 

RMX 

Radius  from  earth  center  to  maximum  altitude 
(km) 

XLMAX 

Length  along  field  line  from  earth's  surface 
to  maximum  altitude  (km) 

MAGLNK  (Magnetic  Dipole 

COSLTO 

Colatitude  of  south  pole  (cosine) 

Common) 

PH  10 

Latitude  of  south  pole  (radians) 

SINLTO 

Colatitude  of  south  pole  (sine) 

VMFO 

Earth  centered  vector  to  north  pole 

XMVO 

Magnetic  moment  (gauss  km  ) 

MODULT  (Modulation 

AMDAT 

Modulation  parameters 

Common) 

LPMIDX 

Modulation  Index 

LPMAOD 

Modulation  type 

NM 

Number  of  modems 

OPTION  (Option  Common) 

IBI ,  ICXI, 
IGI,  IGDI * 
IILI, 

ILCI, 

INDI, 

IUVI,  IXI 

Radiation  options 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

OPTION  (Continued) 

ID1,  ID2, 
KALCD1 , 
KALCD2, 

M01 ,  MD2 

Dummy  variables 

IFPE 

IFPLOT 

IFPP 

IFPS 

I  OUT 

Fcocprint  output  options 

■  Mlil 

1 

Output  options 

KALCFP 

Footprint  calculation  option 

KALCS 

Satellite  circuit  calculation  option 

LFDO 

LFIDO 

LFIN 

LFINT 

LFNO 

LFPLOT 

File  numbers 

MAMB 

MFREGI 

MNATN 

MWIND 

Model  options 

ORIGIN  (Origin  Common) 

IDAYSO 

Day  of  the  month 

IMONSO 

Number  of  the  month 

IYRSO 

Number  of  the  year  in  the  1900s  (1974  is  74) 

PLATO 

North  latitude  (radians) 

PLONO 

East  longitude  (radians) 

SBARO 

Sunspot  number 

ZTO 

Local  zone  time  (decimal  hours) 

OUT  (Output  Common) 

A 

Attachment  rate  (s) 

ALPHAD 

Electron-ion  recombination  coefficient 
(cm3  S-1 ) 

3  -1 

Ion-ion  recombination  coefficient  (cm  s  ) 

ALPHAI 

0 

Detachment  rate  (s) 

PATH  (Propagation  Path 

LPSGOM 

Satellite  geometry  option 

Common ) 

LPTGOM 

Ground  terminal  geometry  option 

NGT 

Number  of  ground  terminals 

NST 

Number  of  satellites 

SORBP 

Satellite  orbit  parameters 
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COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

PATH  (Continued) 

TORBP 

Ground  terminal  motion  parameters 

Earth  centered  vector  to  satellite  (km) 

VSOV 

Earth  centered  satellite  velocity  vector 
(km  s-') 

VTIO 

Earth  centered  vector  to  ground  terminal  (km) 

VTO 

Earth  centered  vector  to  ground  terminal  (km) 

VTOV 

Earth  centered  ground  terminal  velocity  vector 
(km  s-') 

PHEN  (Phenomenology 

PHI 

Common ) 

• 

Temporary  storage  of  quantities  used  in 

■ 

phenomenology  models 

PH6 

PROP  (Propagation 

OBSTOP 

Maximum  circuit  signal  loss  (dB) 

Common ) 

DBSUM 

Absorption  due  to  ionization  (dB) 

DBWVAO 

Absorption  due  to  water  vapor  and  molecular 
oxygen  (dB) 

Jammer  noise  density  (W  Hz”*) 

DJAMS 

DNOISA 

Ambient  noise  density  (W  Hz-1) 

Disturbed  noise  density  (W  Hz"1) 

ONOISE 

DREGDB 

D-region  absorption  (dB) 

ENFAR 

Faraday  rotation  (radians) 

ENSUM 

Electron  density  integral  (cm-2) 

FREQ 

Frequency  (MHz) 

FZERO 

Frequency  selective  bandwidth  (Hz) 

FZEROR 

Frequency  selective  bandwidth  for  reverse  prop¬ 
agation  direction  (Hz) 

IXFBNM 

Fireball  index 

KDIR 

Propagation  direction  flag 

KODE 

Propagation  effects  code 

LJLINK 

Link  index 

LL 

Circuit  index 

LLINK 

Propagation  index 

LNINK 

Noise  index 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

PROP  (Continued) 

Calculation  option 

1  JL 

Sl'.nt  range  between  terminals  (km) 

rms  phase  fluctuation  (radians) 

rms  angle  of  arrival  fluctuations  (radians) 

SPSIR 

rms  angle  of  arrival  fluctuation  for  reverse 
propagation  direction  (radians) 

S4S 

Scintillation  strength  index 

TAUZ 

Decorrelation  time  (s) 

TAUZR 

Decorrelation  time  for  reverse  propagation 
direction  (s) 

VRO 

Earth  centered  vector  to  receiver  terminal  (km) 

VROV 

Earth  centered  receiver  terminal  velocity 
vector  (km/s) 

VSP 

Earth  centered  propagation  direction  vector 
(km) 

VXO 

Earth  centered  vector  to  transmitter  terminal 
(km) 

VXOV 

Earth  centered  transmitter  terminal  velocity 
vector  (km  s"^) 

XVBAR 

Relative  velocity  between  striations  and 
propagation  path  (km  s"1} 

RAY  (Atmospheric 

Wind  Common) 

ARRAY 

Air  parcel  location 

REGINT  (Region  Propa- 

DBCODE 

Absorption  maximum 

gation  Common) 

DBDC 

Dust  cloud  attenuation  (dB) 

DLO 

Minimum  integration  step  size  (km) 

DL2 

Maximum  integration  step  size  (km) 

DSNOM 

Nominal  integration  step  size  (km) 

ENCODE 

Ionization  maximum 

INDMP 

Mandatory  points  flag 

KODEA 

geometrical  region  code 

KODEG 

Incremental  absorption  code 

KODEP 

Ionization  code 

NGAMT1A 

Number  of  low-altitude  ionization  sources 

NHOUT 

Output  flag 

NMP 

Number  of  mandatory  points 

_ L 

NPBFRG 

Number  of  F-region  ionization  bursts 

Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

REGINT  (Continued) 

SGAMMA 

Slant  range  from  transmitter  to  point  of  maxi¬ 
mum  ionization  due  to  low-altitude  ionization 
source  (km) 

SMP 

Slant  range  from  transmitter  to  mandatory 
points  (km) 

SMPMAX 

Maximum  slant  range  to  mandatory  point  (km) 

SMPMIN 

Minimum  slant  range  to  mandatory  point  (km) 

SRBT 

Slant  range  from  transmitter  to  beta  tube 
intersections  (km) 

SRFB 

Slant  range  from  transmitter  to  fireball 
intersections  (km) 

RGEOM  (Geometry  Common) 

Position  vector  (km) 

SPHAT 

Path  direction  vector  (km) 

XMHAT 

Magnetic  pole  direction  vector  (km) 

SIGCOM  (Circuit  Coirmon) 

BWL 

Noise  bandwidth  (MHz) 

DJAMSL 

Jammer  noise  density  (W  Hz”1) 

DNOSA 

Ambient  noise  density  (W  Hz"1) 

DNOSD 

Disturbed  noise  density  (W  Hz"1) 

ENTOT 

_  O 

Electron  density  integral  (cm”  ) 

FREQGH 

Frequency  (GHz) 

FZEROL 

Frequency  selective  bandwidth  (Hz) 

I  MOD 

Modulation  index 

ISIGPR 

Signal  processing  flag 

KALSN 

Signal -to-noise  ratio  flag 

NCIR 

Circuit  index 

SIGFS 

Ambient  signal  strength  (W) 

SIGND 

Disturbed  signal  strength  (W) 

SPHIL 

rms  phase  fluctuations  (radians) 

SPSIL 

rms  angle  of  arrival  fluctuation  (radians) 

S4L 

Scintillation  strength  index 

TAUZL 

Decorrelation  time  (s) 

SOURCE  (Heavy  Particle 
Ionization  Command) 

LPRT 

Ionization  flag 

SPEC  (Fireball  Ioniza¬ 

BAL 

Number  of  aluminum  atoms 

tion  Common) 

BN 

Number  of  nitrogen  atoms 

BO 

Number  of  oxygen  atoms 

Table  5  (Continued) 


COMMON  BLOCK  QUANTITY 


SPEC  (Continued)  BU 


DESCRIPTION 


Number  of  uranium  atoms 


Species  concentration  (cm"J) 


Equilibrium  constants 


SPECEF  (F-region  Ioni¬ 
zation  Common) 


Normalized  values  for  BAL  to  BU 


Species  concentrations  (cm"J) 


SPECQ  (F-region  Ioniza¬ 
tion  Common) 


SPECW  (Device  Spectrum 
Common) 


TVIBN2 
T1F,  T2F 


FEDGD 

FEDHD 

FEDND 

FEDXD 


O 

Pressure  (dynes  cm  ) 

-3  -1 

Ion-pair  production  rate  (cm  s  ) 

Air  density  (gm  cm  ) 

Electron  temperature  (K) 

Molecular  nitrogen  vibrational  temperature  (K) 

Time  at  beginning  and  end  of  deionization 
time  step  (s) 


O 

Species  concentrations  (cm  ) 

Electron  temperature  (K) 

Gas  temperature  (K) 

Molecular  nitrogen  vibrational  temperature  (K) 

Prompt  gamma  yield  fraction 
Hydro  yield  fraction 
Neutron  yield  fraction 
X-ray  yield  fraction 
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Table  5  (Continued 

COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

SPECW  (Continued) 

FFD 

Fission  fraction 

Fraction  of  weapon  mass  in  aluminum 

Fraction  of  weapon  mass  to  uranium 

SPECGD 

Gamma  ray  spectrum  specification 

SPECND 

Neutron  spectrum  specification 

X-ray  spectrum  specification 

WMASSD 

Weapon  mass  (gm) 

SPEEDS  (Atmospheric 
Wind  Common) 

UVW 

Air  velocity  (m  s~^) 

SSYS  (Satellite  Sys¬ 

I  SAD 

Antenra  direction  option 

tems  Common) 

ISASL 

Side  lobe  option 

ISATYP 

Antenna  specification  option 

ISMOD 

Not  used 

ISRN 

Input  noise  specification  option 

ISTMOD 

Modulation  index 

NSANT 

Number  of  antennas 

NSRCV 

Number  of  receivers 

NSXMT 

Number  of  transmitters 

SADAT 

Antenna  quantities 

SASLG 

Side  lobe  gain  relation  to  main  beam  (dB) 

SRBW 

Receiver  noise  bandwidth  (MHz) 

SRNTEM 

Noise  temperature  («) 

STFREQ 

Transmitter  frequency  (MHz) 

STPWR 

Transmitter  power  (W) 

VSAD 

Earth  centered  antenna  direction  vector  (km) 

STRIIN  (Scintillation 

CTHETZ 

Colatitude  of  point  (cos) 

Common) 

DNESB 

Electron  density  variance  (cm"*’) 

DSMX 

Maximum  slant  range  of  ionization  from 
transmitter  (km) 

EN 

Power  law  exponent 

FREQH 

Frequency  (Hz) 

I  STOP 

Striatlon  velocity  flag 

LINK 

Propagation  path  index 

RP 

Radius  from  earth  center  to  field  point  (km) 

RTH 

Equivalent  radius  (km) 

Table  5  (Continued 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

STRIIN  (Continued) 

ST 

Slant  range  between  transmitter  and  re¬ 
ceiver  (km) 

SX 

Slant  range  from  transmitter  to  field  point 
(km) 

VBHAT 

Unit  vector  (earth  centered)  normal  to 
magnetic  field  (km) 

VSPH/ r 

Earth  centered  propagation  direction  vector 
(km) 

VV1 

Earth  centered  velocity  vector  for  trans¬ 
mitter  (km  s“T ) 

VV2 

Earth  centered  velocity  vector  for  fireball 
(km  s- ■ ) 

VV3 

Earth  Centered  velocity  vector  for  receiver 
(km  s-1) 

XU 

Inner  scale  length  (km) 

XLO 

Outer  scale  length  (km) 

TIME  (Origin  Common) 

CHI 

Zenith  angle  of  the  sun  (radians) 

FST 

Fractional  summer  -  one  on  1  July  and  zero 
on  1  January  in  northern  hemisphere  and  re¬ 
versed  in  southern  hemisphere 

FYR 

Fractional  season  year  -  zero  on  1  January  in 
northern  hemisphere  and  zero  on  1  July  in 
southern  hemisphere 

GAT 

Greenwich  apparent  time  (decimal  hours) 

HL 

Local  apparent  time  (decimal  hours) 

Day  of  the  month  in  Greenwich  time  zone 

Time  of  day  (+  1  day,  -  1  night) 

IMONS 

Number  of  the  month  in  the  Greenwich  time  zone 

IYRS 

Number  of  the  year  in  the  1900s  (1974  is  74) 

PLAT 

North  latitude  of  origin  (radians) 

PLON 

East  longitude  of  origin  (radians) 

SBAR 

Average  10.7  -  cm  solar  flux  (l.OxlO-22 

W/m2  Hz)) 

UT 

Universal  time  corresponding  to  the  zone  time 

ZT  (decimal  hours) 

ZT 

Zone  time  for  the  15-degree  longitude  inter¬ 
val  containing  PLON  (decimal  hours) 

TRACE  (Geometry 

AHMIN 

Altitude  of  lowest  air  parcel  in  apogee 

Common) 

_ L 

fireball  (km) 
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Table  5  (Continued) 


COMMON  BLOCK 


TRACE  (Continued) 


TSTEF  (Heavy  Particle 
Ionization  Common) 


TSYS  (Ground  Terminal 
Systems  Common) 


QUANTITY  DESCRIPTION 


Altitude  of  air  parcel  in  apogee  fireball  (km) 
Conjugate  flag 

Altitude  separation  between  air  parcels  in 
apogee  fireball  (km) 

DELV  Velocity  change  between  air  parcels  in  apogee 
fireball  (kms-1) 

HALT  Altitude  of  air  parcel  at  time  T 

RAPOG  Radius  from  earth's  center  to  air  parcel 

apogee  altitude  (km) 

T  Calculation  time  (s) 

Fireball  apogee  time  (s) 

Time  when  air  parcel  reaches  apogee  (s) 

Time  when  air  parcel  reaches  stopping  altitude 

(s) 

VHBOT  Velocity  of  lowest  air  parcel  in  apogee  fire¬ 
ball  at  apogee  (km  s"1) 

XJP  Air  parcel  index 

XLT  Distance  from  earth's  surface  to  air  parcel 

at  time  T  measured  along  magnetic  field  line 
(km) 

.3 

BUF2(1 )  Change  in  electron  concentration  (cm  ) 

BUF2(2)  Change  in  N(*S)  concentration  (cm-3) 

BUF2(3)  Change  In  N(2D)  concentration  (cm-3) 

_3 

BUF2(4)  Change  In  0  concentration  (cm  ) 

3 

BUF2(5)  Change  in  N2  concentration  (cm-J) 

BUF2(6)  Change  In  0?  concentration  (cm-3) 

BUF2(7)  Change  in  N+  concentration  (cm  ) 

BUF2(8)  Change  in  0+  concentration  (cm  ) 

BUF2(9)  Change  in  0  excitation  energy  and  electron 
thermal  energy  per  electron  (ev/electron) 

BUF2( 11)  Change  in  n£  concentration  (cm-3) 

BUF2( 12)  Change  in  ot  concentration  (cm-3) 

^  1 

DAI XX  X-ray  energy  deposited  (?rg  gm  ) 

Fluence  of  the  L th  UV  group  (erg  gm-1) 

FREQ02  Frequency  (MHz) 

I02M0D  Modulation  index 

ITASL  Side  lobe  option 

ITATYP  Antenna  specification  flux 
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Table  5  (Continued) 


COMMON  BLOCK 

QUANTITY 

DESCRIPTION 

TSYS  (Continued) 

ITRN 

Input  noise  specification  option 

ITTMOD 

Modulation  index 

KALCSN 

Signal-to-noise  calculation  option 

NTANT 

Number  of  antennas 

NTRCV 

Number  of  receivers 

NTXMT 

Number  of  transmitters 

STN02 

Signal  to  noise  ratio  (dB) 

TADAT 

Antenna  quantities 

TASLG 

Side  lobe  gain  relative  to  main  beam  (dB) 

TRBW 

Receiver  noise  bandwidth  (MHz) 

TRNTEM 

Noise  temperature  (K) 

TTFREQ 

Transmitter  frequency  (MHz) 

TTPWR 

Transmitter  power  (W) 

WEDEPO  (Heavy  Particle 
Ionization  Common) 

ALPHA1 

Elastic  portion  of  the  heavy-particle  stopping 
cross  section  for  0  on  N  (eV  cm2  atom-1) 

ALPHA2 

Elastic  portion  of  the  heavy-particle  stopping 
cross  section  for  A1  on  N  (eV  cm2  atom-1) 

AM 

Mass  of  mass-15  particle 

BETA1 

Coefficient  of  the  square  root  of  the  energy 
in  the  inelastic  portion  of  the  heavy-particle 
stO£ping  cross-section  for  0  on  N 
(/ 5v  cm2  atom-l ) 

BETA2 

Coefficient  of  the  square  root  of  the  energy 
in  the  inelastic  portion  of  the  heavy-particle 
stojjpinq  cross-section  for  A1  on  N 
(/eV  cm2  atom’l ) 

cm  to  km  (km  cm'1) 

CMTKM 

EMTERG 

MT  to  ergs  (ergs  MT"1) 

MT  to  eV  (ev  MT'1) 
ev  to  ergs  (ergs  ev’1 ) 

EMTTEV 

EPSLN 

IWK 

Number  of  energy  groups  used  in  approximating 
the  heavy  particle  energy  deposition 

SM 

Mass  of  mass-15  particle 

WOG  (Gamma  Deposition 
Common) 

PMG 

Penetrated  air  mass  array  for  which  energy 
deposition  parameters  precomputed  (gm  cm-2) 

UBARGP 

Prompt  gamma  ray  energy  deposition  coefficient 
(cm2  gm-1 ) 

UBARGD 

Delayed  gamma  ray  energy  deposition  coeffi¬ 
cient  (cm2  gm-1) 

Table  5  (Continued) 


COMMON  BLOCK  QUANTITY 


WON  (Neutron  Deposi¬ 
tion  Common) 


lilt 


EXPQNE 


UBARNP 


DESCRIPTION 


Burst  point  nonuniform  air  correction  factor 

Deposition  point  nonuniform  air  correction 
factor 

Mean  free  path  array 

Neutron  elastic  scatter  energy  deposition 
coefficient  (cm2  gm-1) 

Penetrated  air  mass  array  for  which  energy 
deposition  parameters  precomputed  (gm  cm~2) 

Mean  total  cross  section  (cm2  gm"^) 

Neutron  capture  energy  deposition  coefficient 
(cm2  gm-1) 

Prompt  neutron  energy  deposition  coefficient 
(cm2  qm-1) 


WOX  (X-ray  Deposition 
Common) 


WRATE  (D-region  Chem¬ 
istry  Common) 


FCONTX 


UBARX 


XK243 


X-ray  energy  containment  within  mass  depth 
PMXBAR 

O 

Normalizing  air  mass  penetrated  (gm  cm  ) 

Normalized  penetrated  air  mass  array  for  which 
energy  deposition  parameters  precomputed 
(gm  cm-2) 

2  -1 

X-ray  energy  deposition  coefficient  (cm  gm  ) 


Reaction  rate  coefficients 


Initial  ionization  and  dissociation  parameters 


WRATEH  (D-region 
Chemistry  Common) 


Reaction  rate  coefficients 


WSTOR  (Atmospheric 
Wind  Common) 


Day  of  month 
Number  of  month 

Number  of  the  year  in  the  1900s  (1974  is  74) 

Altitude  index 

Local  time  (decimal  hours) 


Table  5  (Continued) 

COMMON  BLOCK  QUANTITY 


DESCRIPTION 


WTOR  (Continued) 

U 

East-west  wind  velocity  (ms-1) 

North-south  wind  velocity  (ms’1) 

V 

WW 

Vertical  wind  velocity  (ms-1) 

ZHTEMP  (Atmospheric 

TZH 

Temperature  array  (K) 

Common) 

CONTROL  ROUTINE 


This  routine  is  the  main  control  routine  for  the  WESCOM  code.  A 
simplified  flow  chart  for  the  routine  is  shown  in  Figure  1. 

The  problem  loop  refers  to  a  given  input  specification  and  the 
problem  number  is  changed  each  time  control  is  transferred  to  the  input 
module  and  new  input  data  obtained. 

After  input  data  are  obtained,  routines  ATMOSI,  MAGFIT,  and  IONOSI 
are  called  to  initialize  the  ambient  atmosphere  and  magnetic  field  models.  Next, 
a  loop  over  calculation  times  is  started  and  routine  PHENOM  is  called  to  obtain 
the  fireball  and  debris  region  geometry  and  intrinsic  properties.  The  location 
of  the  satellites(s)  at  the  calculation  time  is  found  by  calling  routine  SORBIT. 

Next,  a  test  is  made  to  see  if  circuit  or  footprint  calculations 
have  been  requested  and  either  routine  SATC  or  FOOTPC  is  called.  When  the 
calculation  time  loop  is  completed  control  is  returned  to  routine  INPUT  to  obtain 
the  next  problem  case. 
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Figure  1.  Flowchart  for  WESCOM  control  routine. 
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SUBROUTINE  AMBABS 

This  routine  determines  resonance  absorption  due  to  molecular 
oxygen  and  water  vapor  along  a  specified  propagation  path  for  a  specified 
frequency. 

The  absorption  is  determined  by  integrating  the  incremental  absorp¬ 
tion  computed  at  points  along  the  path.  The  integration  is  started  at  the 
lowest  point  on  the  path  and  integration  step  sizes  in  path  length  corresponding 
to  3-km  altitude  steps  are  used  (but  not  larger  than  20-km  path  length  step).  If 
the  lowest  point  on  the  path  is  between  the  transmitter  and  receiver  terminals, 
the  integration  between  the  lowest  point  and  receiver  is  computed  first  and  then 
the  integration  between  the  lowest  point  and  the  transmitter  determined.  The 
integration  in  a  given  direction  is  stopped  when  either  the  incremental  absorption 
goes  to  zero  or  the  product  of  the  incremental  absorption  and  the  remaining  path 
length  becomes  less  than  0.1  dB. 

Routine  GWAT  is  called  to  obtain  the  incremental  absorption  due  to 
water  vapor  and  routine  GOX  is  called  to  obtain  the  incremental  absorption  due 
to  molecular  oxygen.  Inputs  to  these  routines  include  the  atmospheric  pressure, 
temperature,  and  water  vapor  concentration.  These  quantities  are  obtained  at 
each  calculation  point  from  routiie  ATMOSU. 
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SUBROUTINE  ANTND 

This  routine  is  the  driver  for  the  antenna  noise  temperature  module. 
Given  a  receiver  location,  it  obtains  the  ambient  noise  density  (W  Hz"*) 
and  the  noise  density  due  to  thermal  radiation  from  each  fireball  visble  to  the 
receiver.  A  simplified  flowchart  for  the  routine  is  shown  in  Figure  2. 

After  determining  geometry  quantities  associated  with  the  receiver 
location,  the  antenna  and  receiver  characteristics  are  determined.  Then  the 
ambient  noise  density  is  determined  if  the  receiver  noise  characteristics  are 
given  as  input.  If  the  system  noise  characteristics  are  given,  the  ambient  noise 
is  assumed  to  be  included  in  the  input  specification. 

Next,  a  fireball  loop  is  started  and  a  normalized  integral  over  the 
antenna  pattern  is  determined  by  calling  routine  NINT.  Then,  an  index  is  set 
up  which  ranks  the  fireballs  according  to  the  normalized  integral,  in  descending 
order.  This  index  is  used  in  the  next  fireball  loop  to  minimize  the  number  of 
thermal  emission  calculations. 

A  new  fireball  loop  is  begun  to  sum  the  contributions  from  each  fire¬ 
ball  to  t'-e  total  antenna  noise  temperature,  TASUM.  Each  fireball  is  checked 
to  see  if  the  unattenuated  contribution  assuming  a  maximum  fireball  temperature 
is  at  least  five  percent  of  the  current  value  of  the  antenna  noise  temperature. 

Tf  not,  it  is  skipped.  If  it  is,  th6  actual  thernrl  emission  from  the  fireball 
is  determined  by  calling  routine  FBTINT.  After  railing  routine  FBTINT  the  fire¬ 
ball  contribution  is  checked  again  and  further  calculations  skipped  if  its  con¬ 
tribution  is  less  than  five  percent.  If  the  fireball  thermal  radiation  is  still 
significant  the  beta  sheath  absorption  is  determined  for  low-altitude  fireballs 
by  calling  routine  BETAS  and  the  path  absorption  for  both  low-  and  high-altitude 
fireballs  estimated  from  the  D-region  absorption  determined  for  the  path  between 
the  transmitter  and  receiver.  This  absorption  is  generally  a  reasonable  estimate 
of  the  absorption  outside  the  fireball  since  only  fireballs  in  the  main  beam  will 
produce  significant  thermal  noise.  However,  this  approximation  can  be  tested  by 
exercising  the  noise  test  option  as  part  of  the  input  specification.  When  this 
is  done  the  absorption  between  the  fireball  edge  and  the  receiver  antenna  is 
found  by  calling  routine  PROPEN. 
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DETERMINE  RECEIVER  LOCATION 
AND  ANTENNA  AND  RECEIVER 
CHARACTERISTICS 


Figure  2. 


"t  for  subroutine  ANTND. 
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The  final  calculation  in  routine  ANTND  is  the  noise  contribution 
produced  by  emissions  along  the  main  path.  This  is  computed  from  the 
absorption  along  the  path  between  transmitter  and  receiver. 


SUBROUTINE  ARRLIM 


This  routine  determines  the  debris  region  location  and  size  from 
debris  marker  particles  that  have  been  displaced  from  the  debris  stabilization 
location  by  atmospheric  winds. 

The  location  of  the  debris  center  is  found  by  averaging  the  locations 
of  the  marker  particles.  The  minimum  and  maximum  latitude  and  east-longitude 
positions  are  saved  for  use  in  determining  the  debris  radius. 

The  debris  radius  is  determined  as  the  average  of  the  north-south 
and  east-west  dimensions  plus  the  radius  at  debris  stabilization  time.  The 
debris  characteristic  radius  and  distribution  parameter  are  found  for  use  in 
beta  particle  ionization  calculations. 


SUBROUTINE  ATMOS F 

This  routine  computes  properties  of  the  natural  or  heaved  atmosphere 
at  points  above  100  km.  The  routine  is  a  modified  version  of  routine  ATMOSF 
developed  for  the  WOE  code.  A  simplified  flow  chart  for  the  routine  is  shown 
in  Figure  3* 

If  there  are  no  bursts  prior  to  the  calculation  time,  the  properties 
of  the  air  at  the  specified  altitude  are  the  same  as  in  the  natural  atmosphere 
and  are  obtained  from  routine  ATMOSU. 

If  one  or  more  bursts  prior  to  the  calculation  time  can  cause  heave 
(bursts  above  120  km)  and  if  the  operational  mode  discussed  below  is  equal 
to  1,  wake  heave  parameters  are  determined  for  each  fireball.  The  wake  heave 
causes  the  air  motion  below  the  fireball  to  be  consistent  with  the  air  motion 
in  the  fireball.  A  loop  over  bursts  is  made  and  the  energy  deposited  at  110  km 
beneath  the  specified  point  is  calculated  for  each  burst. 

If  any  of  the  bursts  produce  heave  the  Lagrangian  coordinates  of 
the  bottom  and  top  of  each  fireball  producing  wake  heave  are  found  and  then 
routine  HEAVE  is  called.  Routine  HEAVE  determines  the  location  of  the  air  that 
is  at  the  specified  altitude  at  the  calculation  time  and  at  the  time  of  each 
prior  burst  (whether  it  is  a  heave  burst  or  not).  The  heaved  air  is  assumed 
to  move  along  a  vertical  path.  The  HEAVE  routine  also  provides  the  ratio 
of  the  density  at  each  altitude  to  the  density  at  the  initial  altitude  (altitude 
at  time  of  first  burst)  and  the  vertical  mass  scale  height.  Routine  ATMOSU  is 
then  called  to  determine  the  density,  pressure,  temperature,  neutral  particle 
concentrations,  and  mean  molecular  weight  at  the  initial  altitude. 

There  are  four  operational  modes  for  routine  ATMOSF  determined  by  the 
input  parameter  MODE.  For  MODE  =  1,  an  output  array  is  prepared  giving  the 
altitude,  density  or  density  ratio,  and  density  scale  height  at  the  time  of  each 
burst  and  at  the  calculation  time.  The  mass  density  (gm/cm'*)  is  stored  in  the 
array  for  the  initial  altitude  and  the  density  ratio  returned  by  the  HEAVE 
routine  is  stored  for  other  altitudes.  When  the  time  interval  between  bursts 
or  between  the  last  burst  and  the  calculation  time  is  large,  the  altitude  and 
density  ratio  at  intermediate  times  are  found  for  use  in  interpolation  cal¬ 
culations. 


84 


Figure  3.  Flowchart  for  subroutine  ATMOSF. 
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For  MODE  =  0  output  quantities  are  given  for  the  final  altitude 
only.  The  gas  temperature  and  pressure  are  computed  from 


where  the  subscripts  i  and  f  refer  to  the  initial  and  final  altitude 
respectively  and  y  is  the  gas  constant.  The  third  operational  mode  (MODE  =  -1) 
is  used  when  calls  to  ATMOSF  are  made  for  the  same  vertical  column  of  air  and  the 
same  calculation  time  as  a  previous  call,  but  at  a  different  final  altitude.  For 
MODE  =  -1,  the  energy  deposition  matrix  needed  as  input  to  the  HEAVE  routine  is 
not  recalculated.  The  output  for  MODE  =  -1  is  the  same  as  for  MODE  =  0. 

The  final  operational  mode  (MODE  =  2)  is  used  when  calls  to  ATMOSF  are 
made  for  the  same  air  parcel  as  a  previous  call  but  at  an  earlier  time.  The  out¬ 
put  is  the  same  as  for  MODE  =  0  and  also  includes  the  altitude  of  the  air  parcel. 


1 
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SUBROUTINE  ATMOS I 


This  routine  determines  properties  of  the  atmosphere  as  a  function 

of  altitude  at  a  specified  location  and  time  for  use  by  routine  ATMOSU.  The 

routine  is  a  modified  version  of  routines  ATMOSU  and  SPECMIN  developed  for 

ROSCOE.  The  major  modifications  include  combining  the  initialization  computations 

into  a  single  separate  routine  and  the  use  of  an  expontential  formulation  described 
★ 

by  Booker  for  fitting  atmospheric  properties  as  a  function  of  altitude.  The 
formulation  provides  continuous  derivatives  and  replaces  the  high  order  poly¬ 
nomials  used  in  <he  models  developed  for  ROSCOE. 

The  first  step  in  the  atmosphere  initialization  process  is  the 
determination  of  solar  quantities  needed  to  evaluate  the  solar-flux  dependent 
Fourier  coefficients  used  in  computing  the  temperature  gradient  at  the  lower 
boundary  (120  km)  of  the  high-altitude  model  and  the  exospheric  temperature.  The 
required  solar  quantities  are  obtained  by  calls  to  routine  JULIAN,  SOLCYC,  SOLORB, 
and  SOLZEN  (a  call  to  routine  ZTTOUT  used  in  the  ROSCOE  model  to  determine  the 
date  and  time  at  Greenwich  is  made  prior  to  calling  routine  ATMOSI) . 

Next,  routine  TEMPZH  is  called  to  obtain  the  atmospheric  temperature 
profile  at  the  location  and  season  specified  as  input.  This  temperature  profile 
is  used  to  determine  the  molecular  scale  temperature  used  in  evaluating  atmospheric 
intrinsic  properties.  Atmospheric  properties  at  the  boundary  of  the  low-  and  high- 
altitude  models  are  determined  and  the  solar  quantities  are  used  to  evaluate  the 
high-altitude  temperature  model. 

The  remaining  computations  in  routine  ATMOSI  are  calculations  of  minor 
neutral  species  concentrations  at  the  altitude  matrix  given  by  ALTKM  in  the 
labeled  common  block  ALTODN.  This  altitude  matrix  is  different  than  the  one 
(AHT)  used  to  evaluate  the  temperature  and  water  vapor  profiles.  The  calculations 
are  similar  to  those  given  in  routine  SPECMIN  developed  for  ROSCOE  except  that  the 


* 

Booker,  H.G.,  "Fitting  of  Multi-Region  Ionospheric  Profiles  of  Electron 
Density  by  a  Single  Analytic  Function  of  Height,"  J.  Atmos,  Phys . ,  vol  39, 
pp  619-123,  1977. 
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exponential  formulation  is  used  to  obtain  integrals  and  derivatives  and  inter¬ 
polation  of  altitude  dependent  quantities  and  the  species  concentrations  at  the 
matrix  altitude  values  are  the  output  that  will  be  used  by  routine  ATMQSU. 
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SUBROUTINE  ATMOSU 

This  routine  determines  properties  of  the  atmosphere  at  a  specified 
time  and  point  in  space.  The  routine  is  a  modified  version  of  routine  ATMOSU 
developed  for  ROSCOE. 

The  routine  uses  data  prepared  by  routine  ATMOSI  (during  an  initial¬ 
ization  call)  to  determine  atmospheric  quantities  at  a  specified  altitude.  For 
altitudes  below  120  km  a  low-altitude  model  is  used  to  determine  the  intrinsic 
properties  and  major  species  concentrations.  The  basic  parameter  in  the  mode 
is  g/TM,  the  ratio  of  the  acceleration  of  gravity  to  the  molecular  scale 
temperature  determined  in  routine  ATMOSI. 

For  altitudes  above  120  km  a  high-altitude  model  is  used  that 
determines  species  concentrations  and  intrinsic  properties  in  terms  of  the 
temperature  gradient  at  120  km  (TAU) ,  the  exospheric  temperature  (TIF)  and 
the  species  concentrations  at  the  lower  boundary  determined  in  routine  ATMOSI. 

After  the  intrinsic  properties  and  major  species  concentrations  are 
determined,  the  minor  species  concentrations  are  evaluated  if  the  calculation 
flag  JJ  is  greater  than  2.  The  minor  species  concentrations  are  found  by  calling 
routine  BKRIAD  to  interpolate  the  values  determined  in  routine  ATMOSI. 
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SUBROUTINE  BALFIT 

This  routine  fits  a  fourth  order  curve  to  the  ballistic  equation 
solution  for  use  in  determining  the  position  of  an  ionized  air  parcel  on 
a  geomagnetic  field  line. 

The  fit  is  used  for  locations  between  two  points  on  the  field  line 
(2q  and  2^)  given  as  input.  First,  the  velocity  and  acceleration  of  the  air 
parcel  at  2  and  the  time  for  the  parcel  to  reach  2^  are  determined.  Then 
coefficients  C  and  D  are  found  for  use  in  the  following  equation  for  parcel 
position  versus  time  between  2,  and  2^. 

2  =  2q  +  vt-v  +  ct-v2 

♦  C(t-to)3  ♦  D(t-to)4 

where  2  is  the  distance  along  the  geomagnetic  field  line  from  the  earth's 
surface  to  the  air  parcel  and  the  time  t  must  be  between  tQ  and  tj  (time  for 
the  parcel  to  reach  2j). 


90 


SUBROUTINE  BETAS 

This  routine  computes  the  one-way  path  absorption  and  the  back- 
scatter  coefficient  for  normal  incidence  on  the  beta  sheath  ionization  region. 
The  routine  is  a  modified  version  of  routine  BETAS  developed  for  the  WEPH  code. 

The  calculations  are  for  a  sharply  bounded  fireball/debris  geometry 
where  it  is  assumed  that  the  debris  is  within  the  fireball  and  that  the  air 
temperature  and  density  outside  the  fireball  have  ambient  values.  Only  the 
path  absorption  calculation  is  used  in  the  WESCOM  code. 

Beta  sheath  ionization  relations  are  used  to  compute  the  ion-pair 
production  rate,  q,  at  two  points  within  the  beta  sheath:  a  point  at  the 

fireball  surface  and  a  point  outside  the  surface  where  the  ion-pair  production 

_2 

rate  is  reduced  by  about  a  factor  of  e  .  For  each  of  these  points  (assumed 
to  be  at  the  same  altitude)  routine  CHEMDQ  is  called  to  obtain  the  electron 
density  and  electron-neutral  collision  frequency.  Then  the  electron  density 
within  the  sheath  is  approximated  with  an  exponential  function. 


where 


N  .  =  electron  density  at  fireball  surface 
ei 


X  =  distance  from  fireball  surface 
H  =  ionization  scale  height  determined  from 

electron  density  calculations  at  two  points. 
The  one-way  path  absorption  is  found  from 


Substituting  for  N 


Adc  “  4,6  x  10  N  'i 
BS  -  r  e 


7 


2  2 
0)  +  V 


dX 
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dB 


where 


4.6  X  104  NelHu 

*BS '  Jv 

d)  =  wave  frequency  (radians  s  ) 
v  =  collision  frequency  (s  ) 

H  =  scale  height  (km). 
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SUBROUTINE  BTRACE 


This  routine  determines  the  motion  of  an  ionized  air  parcel  along 
a  geomagnetic  field  line  given  the  location  and  velocity  of  the  air  parcel  at 
an  initial  time. 

The  procedure  used  is  to  divide  the  field  line  into  segments  and 
to  calculate  the  time  when  the  parcel  reaches  the  end  of  each  segment  and  the 
parcel  velocity  and  acceleration  at  these  times.  The  parcel  motion  is 
determined  by  calling  routine  BALFIT  which  uses  the  parcel  properties  at  one 
point  to  determine  the  parcel  location  and  properties  at  the  end  of  a  time  step. 

After  the  parcel  properties  are  found  at  equal  intervals  along  the 
field  line  (including  the  portion  between  the  starting  point  and  the  earth)  a 
4th  order  interpolation  procedure  is  used  to  determine  the  location  of  the 
parcel  at  a  specified  input  time.  If  the  parcel  crosses  the  equator  a  flag 

(CONJ)  is  set  and  the  time  for  the  parcel  to  reach  the  equator  is  given  as 

output  in  addition  to  the  time  for  the  parcel  to  reach  the  earth's  surface. 


93 


SUBROUTINE  CHEMD 

This  routine  is  the  driver  routine  for  the  D-region  chemistry  module. 

Th  .outine  is  a  modified  version  of  routine  CHEMD  developed  for  ROSCOE.  Figure 
4  shows  a  simplified  flow  diagram  for  the  routine. 

Inputs  to  the  routine  are  the  location  of  the  D-region  point  at  which 
ionization  is  to  be  calculated  and  the  calculation  time.  Outputs  from  the 
routine  are  electron  density,  total  positive  ion  density,  and  collision 
frequencies. 

Subroutines  ATMOSU  and  IONOSU  are  called  to  find  the  ambient  air 
properties  at  hc«  Subroutine  DRATE  is  called  to  determine  the  values  of  the 
reaction  rate  coefficients.  The  initial  values  of  the  neutral  species  are  set 
equal  to  the  ambient  values  at  hc. 

Next,  routine  DEDEP  is  called  to  determine  the  energy  deposited 

per  gram  at  h  and  t  by  delayed  radiation  and  the  ion-pair  production  rate 
C  c  10 

is  computed  by  multiplying  by  1.79  x  10  times  the  density  and  adding  the 
ambient  ion-pair  production  rate.  Routine  CMDEDT  is  called  to  determine  which 
of  the  previous  bursts  must  be  considered  as  prompt  energy  sources.  For  each 
burst  to  be  considered  the  prompt  energy  deposited  per  gram  and  the  photo¬ 
dissociation  rate  for  ozone  are  stored. 

Next,  a  time  loop  over  each  burst  that  must  be  considered  plus  the 
calculation  time  is  established  and  calculations  are  started  for  the  first  time 
interval.  The  time  interval  is  between  t^  and  t 2  where  tj  is  the  time  of  the 
first  burst  that  must  be  considered  and  t^  is  either  the  time  of  the  next  burst 
that  must  be  considered  or  the  calculation  time  (t£) .  Routine  INITAL  determines 
the  change  in  neutral  species  caused  by  the  burst  at  the  end  of  Phase  2.  Routine 
PH0T0R  determines  the  photodissociation  and  photodetachment  rates  caused  by 
thermal  radiation  from  previous  bursts  at  the  end  of  the  time  interval  . 
Routine  SPECDP  determines  the  neutral  species  concentrations  at  and  Routine 
DTNEP  determines  the  electron  and  total  positive  ion  densities  at  t0. 

After  the  last  time  interval  =  t£)  routines  SPECDQ  and  DTNEQ  are 
called  to  determine  the  effect  of  delayed  radiation  on  the  neutral  and  ion  species 
concentrations.  The  electron  and  positive  ion  densities  returned  by  CHEMD  are 
the  larger  of  those  due  to  prompt  and  delayed  radiation. 
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ROUTINE  CHEMDQ 

This  routine  is  a  driver  routine  to  determine  ionization  at  a 
specified  D-region  altitude.  The  ionization  determined  is  the  same  as  would  be 
obtained  with  routine  CHEMD  for  a  steady-state  ionization  source  if  the  effect 
of  prompt  energy  deposition  on  the  neutral  species  is  neglected.  If  the 
ionization  source  is  not  specified,  an  ambient  value  is  determined  and  used  to 
calculate  the  ionization. 

First,  routine  ATMOSU  is  called  to  obtain  the  ambient  neutral  species 
at  the  a.titude  of  interest  and  routine  DRATE  is  called  to  determine  the  values 
of  the  reaction  rate  coefficients.  Then,  if  the  input  value  of  the  ion-pair 
production  rate  (Q)  is  zero,  routine  IONOSU  is  called  to  obtain  the  ambient 
ion-pair  production  rate. 

Next,  routine  SPECDQ  is  called  to  determine  the  change  in  neutral 
species  and  routine  DTNEQ  is  called  to  determine  the  electron  and  ion  densities 
due  to  the  ion-pair  production  rate.  The  final  calculations  are  for  the  electron 
and  ion  collision  frequencies. 
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SUBROUTINE  CHEMEF 

This  routine  computes  E-  and  F-region  chemistry  for  times  after 
the  prompt  energy  deposition  phase  has  been  completed.  The  routine  is  a 
modified  version  of  routine  CHEMEF  developed  for  ROSCOE. 

Input  and  output  is  provided  through  the  labeled  common  block 
SPECEF.  Common  block  SPECEF  contains  the  initial  and  final  times,  t^  and 
t,,,  for  the  specified  time  interval,  and  a  subscripted  variable  (dimension  2) 
for  each  chemical  parameter.  The  first  array  location  for  each  variable 
contains  the  value  at  time  t^,  and  the  second  location  will  contain  the  up¬ 
dated  value  at  time 

The  operations  performed  by  CHEMEF  are  sequential,  with  very 
little  branching  and  no  looping  or  iteration.  The  operations  are  as  follows: 

1.  Compute  the  heavy-particle  temperature  at  times  tj  and  t^. 

2.  If  the  electron  temperature  TEL(l)  differs  from  the  excita¬ 
tion  temperature  TV1BN2(1),  compute  the  total  excitation 
energy  density  and  equilibrate  it  to  define  a  new  effective 
excitation  energy.  Routine  TEXK  computes  the  new  temperature. 

3.  Compute  the  average  volume  expansion  parameter  and  mean  reaction 
rate  coefficients.  The  rate  coefficients  at  each  temperature 
are  obtained  from  routine  RATE. 

4.  Call  routine  CHMION  to  compute  the  change  in  species  due  to 
the  ion  reactions. 

5.  Compute  the  mean  value  of  each  species  due  to  the  ion  chem¬ 
istry. 

6.  Compute  change  in  species  due  neutral  reactions.  Routine 
GNNO  is  used  to  determine  changes  due  to  reactions  of  atomic 
nitrogen. 

7.  Compute  the  change  in  pressure  and  heavy  particle  temperature 
due  to  the  chemical  energy  release. 

8.  Compute  the  change  in  excitation  energy  density  and  the  new 
equilibration  temperature.  Routine  TEXK  computes  the  new 
temperature.  Set  the  electron  temperature  equal  to  the 
excitation  temperature. 
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ROUTINE  CHEMHR 

This  subroutine  computes  the  electron  and  total  positive  ion  den¬ 
sities  at  a  specified  point  within  a  heated  region  overlay  for  a  specified 
evaluation  time.  This  routine  is  a  modified  version  of  routine  CHEMIR  developed 
for  ROSCOE.  Figure  5  shows  a  simplified  flow  diagram  for  the  routine. 

First,  the  fraction  of  the  ions  that  remain  atomic  (f  )  ,  the 

A 

molecular  ion  production  rate  (q^)  ,  and  the  atomic  ion  production  rate  (q^) 
are  computed.  If  f  >  0.99  and  the  point  of  interest  is  within  the  fireball 
the  ionization  sources  considered  are  nonequilibrium  air  ionization  and  the 
atomic  ion  production  rate.  If  f^  >  0.99  and  the  point  of  interest  is  outside 
the  fireball,  only  the  ionization  due  to  the  molecular  ion  production  rate  is 
considered  and  a  simple  recombination  loss  equation  is  used  to  determine 
the  electron  and  positive  ion  density. 

If  f^  <.  0.99  the  thermal  equilibrium  properties  of  air  and  metal 
species  are  determined  by  calling  routine  EQLAIR  and  EQLMTL  (see  Figure  6). 

An  iteration  loop  is  used  to  insure  species  conservation.  After  each  call  to 
routine  EQLAIR  the  0,  density  is  modified  to  account  for  reactions  with 
nitrogen. 

Next,  effective  reaction  rate  coefficients  are  computed  (see  Figure 
7).  In  computing  the  reaction  rate  coefficients,  the  equilibrium  values  of 
the  species  are  used  if  the  temperature  is  greater  than  T^  (approximately 
2000  K) .  If  the  temperature  is  less  than  T£N  ,  and  Tpj  >  T£N  ,  non- 
equilibrium  values  for  NO  ,  NC^  ,  0  ,  and  0^  are  determined.  If 

Tpi  <.  T^n  ,  ambient  values  for  NO  and  N09  are  used  and  the  0  and  0^ 
arc  assumed  to  be  equal  to  the  thermal  equilibrium  values.  A  partial  correction 
to  the  neutral  species  is  made  to  account  for  effects  of  delayed  radiation.  For 
points  outside  the  fireball  the  effect  of  photodetachment  due  to  fireball 
thermal  radiation  is  determined  by  calling  routine  PH0T0R. 

After  the  effective  reaction  rate  coefficients  are  determined  the 
nonequilibrium  metal  ionization  and  the  equilibrium  ionization  due  to  the 
molecular  ion-pair  production  rate  arc  evaluated.  Then  the  nonequilibrium  air 
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Figure  5.  Flowchart  for  subroutine  CHEMHR. 


Figure  6.  Detail  for  equilibrium  thermal  properties  of  air  and 
metal  species  block  In  Figure  5. 
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Figure  7  Detail  for  effective  rea 
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rate  block  in  Figure  5 


ionization  and  the  ionization  due  to  the  atomic  ion  production  rate  are 
determined. 

The  electron  density  returned  by  routine  CHEMiR  is  the  sum  of 
that  due  to  thermal  ionization  and  that  due  to  delayed  radiation. 
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SUBROUTINE  CHMION 

This  routine  computes  the  change  in  species  concentrations  due  to 
F-  and  F-region  ion  chemistry.  The  routine  is  essentially  the  same  as  routine 
CHMION  developed  for  ROSCOE. 

Initial  values  for  the  species  concentrations  are  obtained  from 
labeled  common  block  SPECEF  and  mean  reaction  rate  coefficients  and  the  ion 
production  rate  are  obtained  from  labeled  common  block  CHEMAN.  The  new  species 
concentrations  at  the  time  t^  will  be  returned  through  SPECEF  common  (under 
variable  subscript  2).  The  energy  released,  DELEV,  is  returned  to  CHP'TF  through 
CHEMAN  common. 

Mean  ion  loss  rates  due  to  each  reaction  are  computed  accounting  only 
for  volume  expansion.  A  check  is  made  to  determine  if  the  N2  or  02  could  be 
significantly  depleted  during  the  time  interval.  If  so,  improved  mean  values 
for  N2  and  02  are  determined.  Next,  if  the  loss  rate  due  to  reaction  with  NO 
is  within  10  percent  of  the  total  N*  loss  rate,  an  improved  mean  value  of  NO  is 
computed. 

In  final  preparation  for  the  ion  solutions  the  total  ion  production  rate 
is  proportioned  into  N+,  0+,  and  M+  production  rates.  The  N+  and  0+  solutions 
are  computed  by  routine  ANLYT2.  The  transient  solution  for  M+  is  computed  by 
routine  RICATT.  This  is  then  combined  with  the  equilibrium  solution  for  M+. 
Finally,  the  change  in  neutral  species  concentrations  is  computed  and  the  energy 
release  determined. 
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SUBROUTINE  CMDEDT 

This  routine  edits  the  bursts  to  form  a  list  for  use  in  D-region 
chemistry  calculations.  The  routine  is  essentially  the  same  as  routine 
CMDEDT  developed  for  ROSCOE.  The  bursts  chosen  are  those  which  affect  neutral 
species  or  ionization  due  to  prompt  energy  deposition  or  photodissociation 
of  0^  .  For  each  burst  the  prompt  energy  deposited  and  the  photodissociation 
rate  for  0^  are  computed  and  stored.  A  flow  diagram  for  routine  CMDEDT  is 
shown  in  Figure  8. 

A  burst  is  included  in  the  burst  list  if  the  prompt  energy  depo¬ 
sition  or  photodissociation  of  0^  changes  the  neutral  species  sufficiently 
to  affect  the  reaction  rate  coefficients  used  for  equilibrium  ionization 
calculations.  The  equilibrium  electron  density  can  be  approximated  by 


where  D  ,  A  ,  and  a  are  effective  reaction  rates.  Tests  are  made 
to  see  if  the  reduction  of  0  ,  N  ,  and  O^A)  results  in  a  detachment 
rate  larger  than 

0 . 1  ( D  ♦  /o!) 

and  whether  the  production  of  0^  and  0  result  in  an  attachment  rate  larger 
than  D  .  In  performing  the  tests  ambient  values  of  the  neutral  species  are 
used  to  compute  A  and  1  .  Since  the  detachment  rate,  D  ,  depends  on  the 
composition  of  negative  ions,  a  fit  to  the  ratio  of  O2  to  total  negative 
ion  density  is  used  and  only  detachment  from  0~  considered.  Changes  in  the 
neutral  species  can  also  affect  the  effective  recombination  coefficient  but 
by  using  an  effective  value  which  is  small  bursts  affecting  a  will  be 
included  in  the  bursts  affecting  the  detachment  rate. 

Routines  EDEPD  and  PHOTOD  are  called  to  determine  the  prompt  energy 
deposition  and  the  photodissociation  rate  for  0,  due  to  thermal  radiation. 
Approximate  relations  are  used  to  estimate  tne  amount  of  0  ,  07 (  A)  ,  and  N 

produced  by  the  burst  and  the  decay  of  these  species  between  the  burst  and 
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calculation  time. 


If  the  burst  does  not  affect  the  neutral  species  sufficiently  to  be 
included  in  the  burst  list,  a  test  on  the  prompt  ionization  remaining  at  the 
calculation  time  is  made  to  see  if  it  is  significant  in  comparison  to  that 
product  by  delayed  radiation.  First  a  test  is  made  to  see  if  a  saturation 
impulse  is  important  by  comparing  the  positive  ion  density  due  to  an  ionization 
impulse  with  the  positive  ion  density  due  to  the  ion-pair  production  rate.  The 
ionization  due  to  prompt  energy  deposition  is  not  important  if 


<  °.l  \ 

at  y  a 

or  rearranged 

t^efa  t  >  10  . 


If  the  burst  could  be  important  a  further  test  is  made  with  the 
actual  ionization  produced  by  the  burst.  For  this  test  the  electron  density 
at  the  calculation  time  due  to  the  impulse  is  estimated  from 


impulse): 


D  +  A  e 


-(A+D)t 


1  +  N  at 
e 


D  +  A 


and  compared  with  the  electron  density  due  to  delayed  radiation. 

The  detachment  rate  used  in  computing  the  electron  density  due  to 
prompt  radiation  is  computed  using  the  neutral  species  concentrations  at  the 
time  of  burst.  This  maximizes  the  electron  density  due  to  prompt  radiation. 
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SUBROUTINE  CMFEDT 

This  routine  edits  the  burst  to  form  a  test  for  use  in  0-  and 

F-region  chemistry  calculations.  The  bursts  chosen  are  those  that  produce 

7-3 

a  specified  level  of  initial  ionization  (set  at  3  x  10  cm  for  the  WESCOM 
code)  at  a  specified  altitude  (currently  chosen  as  200  km  unless  the  lowest 
altitude  on  the  path  is  above  200  km). 

The  air  properties  and  neutral  species  densities  at  the  test  point 
are  found  for  each  burst  from  routine  ATM0SF.  Then  routine  PIONF  is  called  to 
determine  the  initial  ionization  caused  by  the  burst  and  the  result  compared 
to  the  minimum  values  used  to  select  bursts  and  E-  and  F-region  ionization 
sources.  In  computing  prompt  ionization  in  routine  PIONF  air  motion  due  to 
heave  is  neglected.  While  this  can  result  in  overestimating  the  initial  ion¬ 
ization,  it  reduces  the  calculation  time  in  routine  PIONF. 
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SUBROUTINE  DBMHT 


This  routine  determines  time-dependent  quantities  for  debris  regions 
produced  by  high-altitude  bursts. 

For  high-altitude  bursts,  geometrical  debris  regions  are  only 
used  for  gamma  radiation  effects.  For  beta  radiation  the  altitude  and 
dimensions  of  an  ellipse  in  the  horizontal  plane  are  defined. 

Up  to  four  debris  regions  can  be  defined  to  model  the  location 
of  the  fission  debris  for  a  given  burst.  One  debris  region  is  used  to  model 
the  ionized  debris  distributed  within  the  fireball.  Two  debris  regions  are 
used  to  model  debris  settling  down  the  geomagnetic  field  lines  in  the  burst 
and  opposite  hemispheres.  For  bursts  detonated  below  120  km  a  debris  region 
is  used  to  model  neutral  debris  falling  back  over  the  burst  point. 

First  the  ambient  atmospheric  scale  heights  at  the  bottom  of  the 
fireball  and  at  the  debris  settling  altitude  (taken  to  be  150  km)  are  found 
by  calling  routine  ATMOSU.  Then  parameters  related  to  the  fireball  are 
defined  and  for  times  greater  than  the  fireball  apogee  times  the  amount  of 
debris  ch.it  is  in  the  fireball  (FDI),  and  the  amount  of  debris  that  has 
reached  the  settling  altitude  in  each  region  (FDNB  and  FDNC)  are  calculated. 
These  values  are  estimated  from  the  air  parcel  data  calculated  for  the  fire¬ 
ball  region  assuming  that  the  debris  density  distribution  is  similar  to  the 
fireball  ionized  air  mass  distribution.  For  bursts  below  about  120  km 
(FNFUT  >  0.1)  and  times  after  the  magnetic  equilibrium  time  (TEQ),  the 
fraction  of  debris  that  will  fall  back  over  the  burst  point  is  accounted  for. 

If  the  ionized  debris  fraction  is  greater  than  five  percent,  an 
ionized  debris  region  is  defined.  A  skewed  spheroid  debris  geometry  is  used 
and  the  current  fireball  scale  height  is  used  for  the  debris  distribution. 

A  beta  source  plane  is  defined  at  the  altitude  HBETA  determined  in  routine 
FBMHT  that  includes  effects  of  magnetic  freezing  of  the  beta  tube. 

After  the  ionized  debris  region  is  defined,  tests  are  made  to  see  if 
more  than  five  percent  of  the  debris  has  fallen  down  the  field  lines  to  150-km 
altitude.  If  the  amount  is  more  than  five  percent,  a  debris  region  is  defined. 
A  spheroid  geometry  is  used,  but  by  the  use  of  truncation  altitudes  a  pancake- 
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like  region  is  defined.  Beta  source  planes  at  150  km  are  defined  for  use  in 
determining  beta  ray  ionization.  If  the  wind  option  (given  as  input)  is 
exercised,  separate  beta  source  planes  arc  defined  for  each  hemisphere.  Other¬ 
wise,  a  single  beta  sou.ce  plane  in  the  burst  hemisphere  is  defined. 

For  bursts  below  120  km  and  times  after  TEQ  a  debris  region  that 
falls  back  over  the  burst  point  until  reaching  150  km  is  modeled. 
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SUBROUTINE  DBMLT 


This  routine  determines  time-dependent  quantities  for  debris  regions 
produced  by  low-altitude  bursts. 

For  low-altitude  bursts,  a  geometrical  debris  region  is  defined  for 
use  in  computing  effects  of  delayed  beta  and  gamma  radiation.  The  right 
circular  spheroid  (untilted)  and  toroid  geometries  used  for  the  fireball 
(KI\l)F  =  1  and  KINDF  =  3)  are  also  used  for  the  debris  region.  For  times 
before  toroid  formation  the  debris  region  dimensions  are  found  from 

RTD  =  RTBF 

RLD  =  RTBF 

where  RTD  and  RLD  are  the  debris  axis  dimensions  perpendicular  and 
parallel  to  the  vertical  direction  respectively  and  RLF  and  RTF  are  similar 
dimensions  for  the  fireball.  The  quantity  RTBF  is  the  horizontal  debris 
radius  defined  by  routine  ;BMLT.  The  debris  region  is  truncated  with  minimum 
and  maximum  altitude  slices  if  the  fireball  region  is  truncated. 

The  vector  location  of  the  debris  origin  is  determined  and  routine 
VOLUME  is  used  to  obtain  the  debris  volume.  The  debris  stabilization  time 
and  debris  geometry  at  stabilization  time  arc  stored  for  use  in  modeling  wind 
effects.  If  the  debris  stabilizes  above  200  km,  the  vertical  radius  is  limited 
to  50  km  and  if  the  debris  stabilizes  below  200  km  the  debris  vertical  radius 
at  stabilization  time  is  not  allowed  to  be  less  than  4  km. 

Several  additional  quantities  are  calculated  for  use  in  beta 
ionization  calculations.  First,  the  dimensions  of  a  horizontal  ellipse  that 
can  be  used  to  determine  the  maximum  extent  of  the  beta  tube  arc  calculated 
from  the  quantities  RLBF  and  RTBF  defined  in  routine  FBMLT.  Then,  quantities 
required  in  beta  sheath  calculations  are  determined. 
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SUBROUTINE  DBSTAM 


This  routine  is  a  driver  routine  for  determining  post-stabilization 
debris  geometry  due  to  atmospheric  winds. 

The  atmospheric  wind  model  is  used  to  determine  debris  geometry 
for  debris  regions  above  100-km  altitude  and  for  times  after  the  debris  stabili¬ 
zation  times.  For  earlier  times  or  for  debris  altitudes  below  100  km  the  debris 
geometry  determined  by  the  fireball/debris  models  (routines  DBMLT  and  DBM1T)  are 
used.  These  models  include  a  nominal  debris  expansion  due  to  winds,  but  do  not 
include  debris  lateral  motion  (debris  offset)  due  to  atmospheric  winds. 

For  debris  regions  above  100  km  and  times  after  the  debris  stabili¬ 
zation  time  the  geographic  location  of  the  debris  center  and  the  local  time  at 
the  debiis  stabilization  time  (start  time  for  wind  effects  model)  are  computed. 
The  debris  parameters  at  stabilization  time  are  written  out  on  the  detailed  out¬ 
put  file  (file  LFDO)  if  detailed  phenomenology  data  has  been  requested  (input 
option).  Then  routine  PSDM  is  called  with  the  debris  data  at  stabilization  time 
to  determine  the  debris  geometry  at  the  current  calculation  time  due  to 
atmospheric  winds. 
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SUBROUTINE  DEBRIS 

This  routine  is  a  driver  routine  to  determine  the  energy  deposi¬ 
tion  at  a  specified  point  due  to  charge  exchange,  loss  cone,  and  ion  leak 
particles.  Figure  9  shows  a  simplified  flowchart  for  the  routine. 

First,  a  check  is  made  to  see  if  the  deposition  point  is  in  the 
magnetic  conjugate  region  of  the  burst.  If  it  is,  energy  deposition  due  to 
debris  kinetic  energy  sources  is  neglected. 

If  the  deposition  point  is  in  the  burst  region,  the  reduced 
range  for  loss  cone  and  charge  exchange  particles  is  computed  and  the  parti¬ 
cle  energy  lost  in  transversing  the  mass  penetrated  from  the  burst  to  the 
deposition  points  is  found  by  calling  routine  FZET.  Then,  tests  are  made  to 
see  if  the  energy  due  to  loss  cone  particles  is  to  be  computed.  These  tests 
include  an  input  option,  a  test  on  the  amount  of  debris  kinetic  energy  in 
loss  cone  particles,  and  a  test  on  the  upper  intersection  altitude  for  loss 
cone  particles.  Only  the  upper  intersection  altitude  is  tested  (rather  than 
testing  to  see  if  the  field  point  is  between  the  lower  and  upper  intersection 
altitudes  as  is  done  for  beta  particles)  because  vertical  air  motion  (heave) 
can  move  the  air  in  which  the  particles  initially  deposit  their  energy. 

If  the  loss  cone  particle  energy  deposition  is  to  be  computed, 
routine  CONSPC  is  called  to  obtain  the  energy  spectrum  and  the  energy  deposited 
neglecting  the  spatial  distribution  is  computed.  Next,  the  effect  of  the 
spatial  distribution  is  found  by  computing  the  location  of  the  100-km  point  on 
the  geomagnetic  field  line  passing  through  the  deposition  point  of  interest. 

The  distance  of  this  point  from  where  the  field  line  passing  through  the  burst 
point  reaches  1  10  km  is  used  to  compute  the  distribution  function. 

Following  calculations  for  the  loss  cone  particles,  the  energy 
deposition  neglecting  the  spatial  distribution  for  ion  leak  and  charge  exchange 
particles  is  computed.  To  account  for  the  spatial  distribution  for  charge 
exchange  and  ion  leak  particles  a  plane  normal  to  the  geomagnetic  field  lines 
at  the  burst  point  is  used  as  the  distribution  plane.  The  position  of  .ne 
magnetic  field  line  that  passes  through  the  field  point  in  the  distribution 
plane  is  determined  for  use  in  the  distriFution  function. 
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COMPUTE  REDUCED  RANGE  AND 
ENERGY  CORRESPONDING  TO 
PM  FOR  LOSS  CONE  AND 
CHARGE  EXCHANGE  PARTICLES 


DETERMINE  ENERGY  SPECTRUM 
A',0  DEPOSITION.  MODIFY 
DEPOSITION  TO  ACCOUNT  FOR 
SPATIAL  DISTRIBUTION 


DETERMINE  ENERGY  SPECTRUM. 
COMPUTE  REDUCED  RANGE  AND 
ENERGY  CORRESPONDING  TO  PM. 
DETERMINE  ENERGY  DEPOSITION 


^  CHARGE  ^ 
EXCHANGE 
CALCULATIONS? 


DETERMINE  ENERGY  SPECTRUM 
ANO  ENERGY  DEPOSITION 


MODIFY  ION  LEAK  AND  CHARGE 
EXCHANGE  ENERGY  DEPOSITION 
TO  ACCOUNT  FOR  SPATIAL  DISTRIBUTION 


Figure  9.  Flowchart  for  subroutine  DEBRIS. 


SUBROUTINE  DEDEP 

This  routine  determines  the  delayed  energy  deposition  rate 
(ergs  gm  *  s  *)  due  to  weapon  radiation. 

A  loop  over  debris  regions  is  initiated  and  routines  EDEPG  and 
EDEPB  are  called  to  obtain  the  energy  deposited  per  gram  by  gamma  rays  and 
beta  particles,  respectively.  The  call  to  routine  EDEPG  is  skipped  if  the  air 
density  at  the  field  point  is  less  than  10  ^gm  cm  ^  (about  125  km)  and  the 
call  to  routine  EDEPB  is  skipped  if  the  field  point  is  outside  the  beta  tube. 

Delayed  energy  deposition  due  to  neutrons  is  computed  by  setting 
up  a  burst  loop  and  calling  routine  EDEPND.  Routine  EDEPND  determines  the 
neutron  energy  deposition  due  to  capture  reactions  and  elastic  scattering 
reactions  occurring  after  the  time  used  for  prompt  energy  deposition. 
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SUBROUTINE  DSTINT 


This  routine  determines  the  attenuation  due  to  progagation  through 

dust  clouds. 

A  loop  over  dust  clouds  defined  by  phenonemonology  module  is 
started  and  for  each  dust  cloud  routine  PCLOUD  is  called  to  determine  the 
geometry  for  each  dust  cloud  particle  group.  Then  routine  RINTER  is  called  to 
see  if  the  ray  path  intersects  the  region.  If  it  does,  routine  PGROUP  is  called 
to  determine  the  extinction  coefficient  and  the  attenuation  is  calculated. 

If  detailed  output  has  been  requested  the  dust  cloud  intersections  and 
attenuations  are  written  out  on  file  LFDO, 
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SUBROUTINE  DTNEP 

This  routine  determines  the  electron  and  total  positive  ion  den¬ 
sities  at  a  specified  time  after  burst  due  to  prompt  radiation.  The  routine 
is  similar  to  routine  DTNEP  developed  for  ROSCOE. 

First  a.  and  A  and  initial  values  of  a.  and  D  are  determined. 
1  d 

rhese  values  are  used  to  compute  initial  values  for  Ng  and  N+  . 

The  initial  value  of  Ng  is  tested  against  the  electron  density 

that  would  be  produced  by  delayed  radiation  for  the  same  reaction  rate 

coefficients.  If  the  electron  density  due  to  prompt  radiation  is  less  than  1/10 

the  value  due  to  delayed  radiation  or  if  the  electron  density  due  to  prompt 

3  -3 

radiation  is  less  than  10  electrons  cm  ,  the  electron  and  positive  ion 
densities  due  to  prompt  radiation  are  set  to  zero  and  no  further  calculations 
are  made. 

If  after  the  above  tests  prompt  ionization  can  still  be  important, 
the  positive  and  negative  ion  compositions  and  new  effective  values  for  ad 
and  D  are  computed.  These  revised  values  are  then  used  to  recompute  the 
electron  and  total  positive  ion  densities. 
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SUBROUTINE  DTNEQ 

This  routine  determines  the  steady-state  reaction  rate  coefficients 
and  the  corresponding  electron  and  total  positive  ion  densities  due  to  delayed 
radiation.  The  routine  is  essentially  the  same  as  routine  DTNEQ  developed  for 
ROSCOE. 

The  ion-pair  production  rates  producing  0*  and  NO*  are  determined 

from  the  total  ion-pair  production  rate.  Next  attachment  and  detachment  rates 
are  computed.  The  detachment  rate  includes  the  effects  of  photodetachment 
from  thermal  radiation  from  each  burst  occuring  prior  to  the  calculation  time. 
Then  positive  and  negative  ion  transfer  rates  are  determined  including  effects 
of  photodissociation  from  thermal  radiation. 


Initial  values  for  the  ratios  R  through  R  ^  (positive  ion 
ratios)  and  R  j  through  (negative  ion  ratios)  are  chosen  and  used  to 

compute  effective  reaction  rate  coefficients  and  equilibrium  electron  and  ion 
densities.  The  electron  and  ion  densities  are  then  used  to  recompute  the 
ratios  R^  through  R^0  and  R^  through  R^g  .  The  above  procedure  is 
repeated  until  the  effective  recombination  coefficient  and  effective  detachment 
rate  change  by  less  than  1  percent. 


SUBROUTINE  DUSTMI 


This  routine  computes  time-independent  quantities  used  in  dust 
region  calculations. 

First,  the  scaled  burst  height  (ratio  of  burst  height  to  the 
pressure  equilibrium  radius  computed  in  routine  FBMLI)  is  tested  to  see  if 
the  burst  will  cause  a  dust  cloud.  If  the  scaled  burst  height  is  less  or 
equal  to  3  a  dust  cloud  is  defined.  For  the  first  cloud  defined  the  particle 
sice  group  parameters  are  computed  (same  for  all  dust  clouds).  Eight 
particle  groups  are  used  in  the  current  model  (see  Table  6),  with  the 
boundaries  chosen  to  keep  the  spatial  overlap  of  adjacent  groups  relatively 
uniform.  The  fraction  of  the  particles  in  each  size  group,  FNP(I),  is  computed 
from  the  assumed  size  distribution  with  0.001-  and  10.0-cm  minimum  and  maximum 
particle  sizes. 

For  each  dust  cloud  the  mass  loading  factor,  XMLOD,  is  computed  and 
the  diameter  of  the  largest  particle  lofted,  ATRY,  is  computed  rounded  to  the 
nearest  particle  group  boundary;  ie,  the  nearest  value  of  ABAR.  The  index  of 
the  largest  group  lofted,  IMAXP,  is  then  the  group  for  which  the  maximum  size 
is  the  rounded  value  of  ATRY.  Finally,  the  total  number  of  particles  included 
in  all  the  groups,  designated  TNP,  is  computed.  It  should  be  noted  that  all 
these  particles  are  not  lofted;  only  those  contained  in  particle-size  groups 
up  through  IMAXP  are  lofted. 
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Table  6.  Time-independent  particle  group  parameters. 


Group 

Index 

I 

Minimum 

Di ameter 
(cm) 

ABARD ) 

Maximum 

Diameter 

(cm) 

ABARU+1) 

Fraction  of 
Particles 

FNP(I) 

1 

0.001 

0.09 

» 1.000 

2 

0.09 

0.4 

1.356E-6 

3 

0.4 

1.0 

1.463E-8 

4 

1.0 

2.0 

8.750E-10 

5 

2.0 

3.25 

9.587E-1 1 

6 

3.25 

5.25 

2.222E-1 1 

7 

5.25 

7.5 

4.540E-12 

!  8 

7.5 

10.0 

1.370E-12 
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FUNCTION  EDEPB 

This  routine  computes  the  energy  deposition  rate  (ergs  gm  *  s”1) 
due  to  beta  particles.  The  model  is  based  on  formulations  developed  for  ROSCOE. 

A  simplified  flowchart  for  the  routine  is  shown  in  Figure  10. 

The  average  beta  particle  energy,  beta  particle  release  rate,  and 
energy  deposition  coefficient  are  found  from  routine  WOBD.  Then  routine 
ATMOSU  is  called  to  determine  atmospheric  properties  at  the  beta  source  altitude. 
If  the  debris  is  distributed  in  a  constant  altitude  plane,  a  correction  term 
is  found  for  use  in  computing  the  debris  integral  along  the  geomagnetic  field  line 
through  the  field  point.  The  correction  term  is  a  function  of  the  location  of  the 
field  line  at  the  beta  source  altitude. 

For  beta  source  altitudes  above  80  km  a  debris  pancake  (elliptical 
region  in  a  constant  altitude  plane)  is  used  to  model  the  beta  source.  For 
lower  beta  source  altitudes,  the  actual  debris  shape  (spheroid  or  toroid)  is  used. 
The  beta  source  geometry  is  used  to  determine  the  integral  of  debris  along  the 
geomagnetic  field  line  intersecting  the  field  point,  the  altitude  of  the  debris 
center,  and  the  magnetic  field  strength  at  the  debris  center.  For  .?icld  points 
within  the  debris  region  (source  altitudes  below  80  km)  the  debris  region  is 
divided  into  two  parts;  one  below  and  one  above  the  field  point. 

The  integrals  are  used  to  compute  the  energy  deposition  within  the  beta 
tube.  For  beta  source  altitudes  below  80  km  the  energy  deposited  within  the 
beta  sheath  is  computed  if  requested  (input  option  to  routine).  The  beta  sheath 
energy  deposition  is  computed  and  added  to  that  calculated  for  the  beta  tube. 
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Figure  10.  Flowchart  for  function  EDEPB 


FUNCTION  EDEPG 

This  routine  computes  the  energy  deposition  rate  (ergs  gm-1  s-1J 
at  a  given  point  i '  space  due  to  gamma  radiation  from  the  fission  debris.  A 
simplified  flowchart  for  routine  EDEPG  is  presented  in  Figure  11. 

First,  the  debris  truncation  altitudes  and  the  vertical  dimensions 
are  calculated.  Then,  the  mass  penetrated  from  the  field  point  to  the  debris 
center  is  computed.  If  a  point  source  calculation  has  been  requested  (input 
option),  the  mass  penetrated  and  slant  range  between  the  field  point  and  debris 
center  are  used  to  compute  the  gamma  ray  energy  deposition. 

If  a  point  source  has  not  been  requested,  the  dimensions  and  volume 
of  the  debris  region  are  found  and  the  mass  penetrated  is  computed  to  the  debris 
edges.  Then,  tests  are  made  to  see  if  analytical  solutions  can  be  used  or 
whether  a  numerical  solution  is  required  to  compute  the  energy  deposition.  A 
numerical  solution  is  used  when  the  debris  geometry  is  a  toroid  or  tilted  spheroid 
or  has  been  truncated  so  that  the  top  is  below  the  center,  when  the  debris  vertical 
or  radial  distribution  is  not  uniform,  and  when  the  step  sizes  based  on  allowable 
changes  in  mass  penetrated  are  smaller  than  the  debris  size. 

The  coordinate  system  used  in  computing  the  step  sizes  is  one  where 
the  r  coordinate  is  measured  along  a  constant-altitude  plane  at  the  altitude 

of  the  debris  center.  The  Ar  and  Az  step  sizes  for  which  the  change  in  mass 

_2 

penetrated  from  the  field  point  is  less  than  AM  gm  cm  are  determined  from 
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Figure  11.  Flowchart  for  function  EDEPG. 
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z =  z  coordinate  of  top  of  debris  region 

Tj  ,  r^  =  minimum  and  maximum  horizontal  distances  from  field  point 
to  debris  plane  at  z  =  zq 

Mj  ,  M2  =  mass  penetrated  from  field  point  to  r^  ,  zq  and  r ,  zq 
a  *  radius  of  sphere,  radius  of  toroid  arm,  or  semithickness  of 

pancake. 

If  Ar  and  Az  are  greater  than  the  horizontal  and  vertical  dimensions 
of  the  debris  region,  A^(M)  is  essentially  constant  over  the  debris  region.  If 
Ay(M)  is  constant  and  the  slant  range  to  the  debris  center  is  greater  than 
three  times  the  largest  horizontal  debris  radius,  a  point-source  approximation 
is  used  to  compute  the  ion-pair  production  rate.  If  A^(M)  is  constant  and  the 
slant  range  is  smaller  than  above,  analytical  expressions  are  used  for  spherical 
and  nancake  debris  geometries.  For  all  other  conditions,  a  numerical  integration 
procedure  is  used. 

When  the  numerical  integration  procedure  is  used,  a  second  set  of  Ar 
and  Az  step  sizes  is  determined  to  ensure  that  the  numerical  integration  gives 
good  agreement  with  the  volume  of  the  debris  region.  These  step  sizes  are  found 
from  the  following: 
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2a 

where  b  is  the  largest  horizontal  radius  of  the  debris  region. 

The  Az  step  size  is  taken  as  the  smaller  of  that  determined  on  the 
basis  of  the  variation  in  A^(M)  and  that  determined  for  volume  accuracy,  with 
a  minimum  step  size  of  1  km. 

A  minimum  Ar  step  size  is  then  chosen  from 


124 


where  equals  the  mass  density  at  the  debris  center,  to  restrict  the  step 

size  for  low  debris  altitudes  where  the  mass  penetrated  is  a  strong  function  of 
r  .  As  a  further  restriction,  the  minimum  Ar  step  size  is  not  allowed  to 
exceed  500  km  for  debris  altitudes  less  than  50  km  or  1000  km  for  debris 
altitudes  greater  than  50  km. 

For  each  Az  step  in  the  integration  procedure,  the  mass  penetrated 
to  two  points  in  the  debris  region  at  altitude  z  is  found  with  routine  PMASS. 
A  new  Ar  step  size  is  computed  from 

(s2  -  Sj)AM 
Ar  =  M1<ln(M2/M1) 

where  s^  is  the  slant  range  to  a  point  in  the  debris  region  closest  to  the 
field  point,  and  s2  is  the  slant  range  to  a  point  chosen  to  make  M2  - 
_2 

about  85  gm  cm  .  The  minimum  of  the  above  Ar  step  size  and  the  one 
previously  determined  is  used  for  the  integration  over  r  .  The  mass  pene¬ 
trated  to  each  volume  element  between  r.  and  r»  is  found  by  interpolating 
between  Mj  and  M2  : 


if  r  becomes  larger  than  r2  ,  the  two  points  r^  and  r2  are  changed  so 
that  r2  becomes  Tj  ,  and  r2  is  again  chosen  so  that  M2  -  Mj  is  about 
85  km.  Within  each  integration  region,  the  Ar  step  size  is  uniform  except 
where  it  is  necessary  to  reduce  the  step  size  to  match  a  debris  boundary.  The 
Az  step  size  is  always  uniform. 

To  eliminate  numerical  integration  when  the  ion-pair  production  rate 
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q  is  small,  a  point-source  estimate  is  made  at  each  step  in  z  with  the 
debris  at  r^  ,  2  ,  If  the  point-source  calculation  of  q  is  less  than  a 
minimum  value  (chosen  as  one-tenth  of  the  natural  ion-pair  production  rate), 
integration  over  r  is  not  performed. 

The  vertical  debris  distribution  is  accounted  for  in  computing  the 
debris  density  at  each  step  in  z  .  For  times  after  debris  stabilization 
the  debris  may  have  a  radial  distribution  defined  by  the  debris  distribution 
parameter.  This  is  approximately  accounted  for  in  the  numerical  integration 
by  first  increasing  the  radial  dimension  of  the  debris  to  include  debris  in 
the  tail  of  the  distribution  and  by  modifying  the  contribution  to  the  energy 
deposition  integral  for  each  step  in  r  to  account  for  the  debris  distribution. 

The  last  calculation  in  the  routine  is  to  account  for  the  reduced 
mass  in  the  fireball  if  the  debris  is  within  a  low-altitude  fireball.  Corrections 
to  the  energy  deposited  at  the  field  point  due  to  the  reduced  mass  in  the  fire¬ 
ball  are  made  when  the  fireball  associated  with  the  debris  region  is  below  40  km 
and  the  fireball  temperature  is  greater  than  300  K. 


FUNCTION  EDEPND 

This  routine  computes  the  energy  depositon  rate  (ergs  gm  1  s*1) 
at  a  given  point  in  space  due  to  neutron  elastic  collisions  and  capture 
reactions. 

Before  computing  the  neutron  energy  deposition,  tests  are  m-»de  to  see 
if  the  energy  deposition  would  be  negligible  in  comparison  to  gamma  and  beta 
energy  deposition  rates.  Detailed  neutron  calculations  for  a  particular  burst 
are  not  made  if 
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and  f  is  the  fraction  of  the  weapon  yield  in  neutron  kinetic  and  capture 
energy. 

For  bursts  considered,  the  mass  penetrated  from  the  burst  point  to  the 
deposition  point  and  the  vertical  mass  above  the  burst  point  used  arc  poor  for 
such  cases  and,  in  general,  other  ionization  sources  will  dominate. 

Once  the  shortest  distance  has  been  found,  the  neutron-decay  beta  flux 
and  the  ion-pair  production  rate  are  determined. 
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SUBROUTINE  ENEF 

This  routine  computes  the  electron  density  (cm-3)  and  electron 
collision  frequencies  (s’  )  at  a  specified  location  above  100  km.  A 
simplified  flowchart  for  subroutine  ENEF  is  presented  in  Figure  12. 

First,  the  location  and  properties  of  the  air  that  is  at  the  cal¬ 
culation  point  at  the  calculation  time  are  found  at  the  time  of  all  bursts 
prior  to  the  calculation  time  and  at  the  calculation  time  by  calling 
routines  ATMOSF  and  IONOSU.  Then  the  ion-pair  production  rate  at  the  cal¬ 
culation  time  due  to  gamma  rays  and  beta  particles  is  found  from  routine 
DEDEP. 

The  first  F-region  prompt  source  is  found  by  testing  logical 
parameter  LPRT,  which  was  set  up  in  routine  CMFhDT.  The  prompt  ionization 
due  to  this  burst  is  computed  by  routine  PIONF.  Routine  PIONF  also  modifies 
the  neutral  particle  concentrations  according  to  the  number  of  ions  produced 
and  determines  a  new  excitation  temperature. 

The  end  of  the  deionization  time  interval  is  taken  to  be  the  next 
prompt  source  (or  the  calculation  time  if  there  are  no  more  prompt  sources), 
subject  to  the  constraint  that  the  mass  density  of  the  air  cannot  change  by 
more  than  a  factor  of  two  during  the  interval.  If  the  density  ratio  during 
an  inter-burst  time  interval  changes  by  more  than  a  factor  of  two,  an 
intermediate  time  is  taken  such  that  the  density  change  to  this  time  is 
exactly  a  factor  of  two.  An  exponential  interpolation  is  used  to  compute  the 
density  at  the  intermediate  time  and  linear  interpolation  to  compute  the 
altitude  of  the  air  parcel.  The  remaining  interval  will  be  further  checked 
on  the  next  iteration  and  broken  up  again  if  necessary.  In  addition,  to 
prevent  large  time  steps  immediately  after  burst,  calculations  are  made  at 
1  and  3  seconds  after  each  burst. 

To  avoid  the  possibility  of  overlooking  large  peaks  or  nulls  in 
the  Jensity  time-history  during  a  long  inter-burst  time  interval,  an  array 
of  intermediate  times  and  the  properties  of  the  air  at  those  times  has  been 
computed  during  the  call  to  ATMOSF,  above.  These  times  are  merged  with  the 
burst  time  array  for  use  in  breaking  up  time-intervals  and  interpolating  to 
find  the  properties  of  the  air  at  intermediate  times. 
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Figure  12.  Flowchart  for  subroutine  ENEF. 
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Once  the  time-interval  for  deionization  is  established,  the 
changes  in  ion  and  neutral  species  densities  are  founJ  from  routine  CHEMFF. 

If  the  mass  density  decreases  during  the  time  interval,  the  electron,  ion, 
and  neutral  particle  concentrations  are  modified  due  to  volume  expansion 
and  the  gas  temperature  is  decreased  due  to  adiabatic  expansion  in  routine 
CHEMEF.  If  the  density  increases,  the  particle  concentrations  and  gas 
temperature  are  modified  due  to  cold  air  entrainment  and  volume  contraction. 

A  new  excitation  temperature  is  obtained  from  routine  TEXK. 

If  the  end  of  the  time  interval  just  computed  is  an  intermediate 
time,  the  next  time  interval  is  taken  to  the  time  of  the  next  prompt 
source.  Control  is  returned  to  the  block  which  determines  whether  the 
remaining  interval  need  be  subdivided  again.  If  the  end  of  the  time  interval 
just  completed  is  not  an  intermediate  time,  control  is  returned  to  the 
beginning  of  the  burst  loop.  If  there  is  another  prompt  ionization  source, 
the  the  prompt  ionization  from  that  source  is  computed  and  all  the  sub¬ 
sequent  calculations  repeated. 

If  there  are  no  more  prompt  sources,  the  electron  density  due 
to  natural  and  weapon  ion-pair  production  rates  is  found  from  routine  CHEMQ 
with  neutral  and  ion  species  for  the  most  recent  time  interval.  The  electron 
density  at  the  calculation  time  is  then  computed  as  the  maximum  of  that  due 
to  prompt  radiation  and  the  ion-pair  production  rates.  The  logical  parameter 
ION  is  set  to  equal  to  1,  2,  3,  or  7,  depending  on  whether  the  dominant 
ionization  source  is  natural,  beta  particle,  gamma  ray.,  or  prompt 
ionization,  respectively. 

The  final  step  in  routine  ENEF  is  the  computation  of  electron- 
neutral  and  electron-ion  collision  frequencies. 
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ROUTINE  FBEN 

This  routine  determines  fireball  ionization,  collisi>n  frequencies, 
and  incremental  absorption  at  a  specified  point  within  the  fireball,  A 
simplified  flowchart  for  the  routine  is  shown  in  Figure  13. 

For  fireballs  produced  by  low-altitude  bursts  and  for  fireballs 
produced  by  high-altitude  bursts  for  times  less  than  the  time  to  reach  apogee 
(t  )  the  air  at  the  calculation  point  is  assumed  to  have  come  from  a  position 
in  the  initial  fireball  which  was  at  the  same  number  of  scale  heights  above  the 
bottom  as  the  calculation  point  is  at  the  calculation  time.  The  initial  mass 
density  is  taken  as  the  atmospheric  density  at  the  initial  point  at  the  time  of 
the  burst  and  is  calculated  by  routine  ATMOSU  or  ATMOSF. 

The  mass  density  at  the  calculation  time  is  computed  from  the  initial 
mass  density  and  the  ratio  of  the  initial  and  current  fireball  scale  heights. 

For  fireballs  produced  by  low  altitude  bursts  the  mass  density  is  taken  as  the 
larger  of  the  above  and  the  pressure  equilibrium  value.  If  the  pressure 
enui librium  value  is  larger,  the  additional  mass  density  is  assumed  to  be  due 
to  entrainment. 

For  fireball  reference  altitudes  below  100  km  the  ion-pair  production 
rate  caused  by  beta  particles  and  gamma  rays  is  found  from  routine  DEDEP.  For 
toroid  geometries  it  is  possible  for  the  calculation  point  to  be  within  the  hole 
of  the  toroid.  For  this  case  routine  CHEMDQ  is  used  to  compute  the  steady-state 
ionization  neglecting  effects  of  prompt  ionization  and  assuming  ambient 
temperature. 

For  calculation  points  within  the  fireball  routine  CHFMHR  is 
called  to  obtain  the  ionization  and  collision  frequencies.  If  the  burst 
altitude  is  below  about  85  km  and  the  calculation  point  is  within  the  debris 
\region,  the  metal  species  densities  are  computed  prior  to  calling  routine 
CHEMIR.  Also,  if  the  fireball  temperature  is  below  2000  K  the  air  density  and 
entrainment  fraction  at  the  time  the  air  wis  at  2000  K  the  air  are  computed. 

For  fireballs  produced  by  high-altitude  bursts  and  calculation  times 
after  the  time  to  reach  apogee  the  fireball  ionization  is  calculated  as  the  sum 
of  ballistic  ionization  and  diffusion  ionization.  At  the  first  call  to  routine 
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Figure  13.  Flowchart  for  subroutine  FBEN. 


FBEN  after  t  the  apogee  fireball  ionization  along  the  longitudinal  axis 
is  fit  with  two  exponential  terms.  The  exponential  parameters  are  determined 
from  the  models  described  above.  Then  a  test  is  made  to  see  if  the  field  point 
is  within  the  ballistic  ionization  region.  If  it  is,  the  ballistic  ionization 
is  computed  by  finding  the  location  within  the  apogee  fireball  where  the  air 
at  the  calculation  point  came  from  and  computing  the  ionization  at  apogee  and 
the  change  in  ionization  since  apogee. 

Diffusion  ionization  parameters  are  calculated  once  for  each 
calculation  time.  Then,  the  diffusion  ionization  is  determined  and  combined 
with  the  ballistic  ionization.  For  field  points  near  the  top  of  the  fireball, 
tnc  electron  density  is  reduced  empirically  to  approximate  additional  volume 
expansion  not  modeled.  Also,  a  transverse  radial  ionization  distribution 
is  used  to  modify  the  ionization  in  terms  of  the  transverse  location  of  the 
field  point. 

After  obtaining  the  ionization  and  collision  frequencies,  the 
incremental  absorption  is  computed.  It  the  electron  density  is  large  enough 
that  the  plasma  frequency  is  greater  than  the  wave  frequency,  the  incremental 
absorption  is  set  equal  to  1  x  10*^  dB  km"*. 

When  footprint  calculations  are  made  the  electron  density  and 
temperature  are  computed  at  several  points  along  the  longitudinal  axis  and 
stored  in  common  block  CHEMFB  as  part  of  preliminary  calculations.  Then, 
these  data  are  used  to  fit  the  electron  density  and  temperature  at  points 
within  the  fireball  that  are  along  propagation  paths  in  each  scan.  The 
logical  flag  LFPCS  is  set  to  an  integer  greater  than  0  when  the  interpolation 
calculations  are  to  be  used  in  routine  FBEN. 
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SUBROUTINE  FBI NT 

Routine  FBINT  is  a  driver  routine  for  the  module  that  determines 
environment  and  propagation  quantities  for  propagation  through  fireball  regions. 
A  simplified  flowchart  for  the  routine  is  shown  in  Figure  14. 

After  integral  values  are  initialized  to  zero  for  each  fireball  a 
fireball  loop  is  set  up  and  routine  RINTER  is  called  to  determine  whether  the 
propagation  path  intersects  the  fireball.  If  an  intersection  exists,  fireball 
quantities  are  stored  in  labeled  common  block  CHEMB.  For  low-altitude  bursts 
the  metal  species  densities  (aluminum  and  uranium)  arc  computed  by  adding  the 
metal  densities  from  each  debris  region  associated  with  the  fireball.  Path 
intersections  with  debris  regions  are  also  found  for  low-altitude  bursts  for 
later  use  in  selecting  path  integration  segments.  A  beta  sheath  ionization 
flag  is  set  for  debris  regions  produced  by  the  same  burst  that  produced  the 
fireball.  Beta  sheath  ionization  from  other  debris  regions  is  ignored. 

Next,  path  integration  segments  arc  determined.  For  fireballs 
from  high-altitude  bursts  one  segment  is  used  for  burst  region  intersections 
and  another  for  conjugate  regions  intersections.  For  fireballs  from  low- 
altitude  bursts  up  to  eight  path  segments  can  be  defined  to  describe  different 
ionization  regions  within  the  fireball. 

After  the  integration  segments  have  been  defined  a  numerical 
integration  of  propagation  quantities  is  performed  for  each  segment.  A  running 
sum  over  the  segments  and  over  fireballs,  if  the  path  intersects  more  than  one 
fireball,  is  kept  and  added  to  the  path  propagation  quantity  terms. 
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FBINT 


Figure  14.  Flowchart  for  subroutine  FBINT. 
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SUBROUTINE  FBMHI 


This  routine  determines  time-independent  scaling  quantities  for 
fireballs  from  high-altitude  bursts. 

Routine  FBMil  is  called  for  each  burst  detonated  above  75  km.  The 
high-altitude  model  flag,  LMCF,  is  set  'to  zero  indicating  a  high-altitude  burst 
but  may  be  reset  to  one  (a  low-altitude  burst)  later  as  indicated  below. 

First,  the  debris  loss-cone  and  ion-leak  parameters  are  found  by 
calling  routines  LOSCON  and  IONLEK.  This  includes  the  amount  of  debris  kinetic 
energy  in  loss-cone  and  ion-leak  particles  which  is  not  available  to  produce 
the  fireball.  Next,  charge  exchange  and  UV  fractions  are  found  and  the  heave 
flag  is  set  (no  heave  for  bursts  detonated  below  120  km). 

Next,  the  initial  fireball  horizontal,  downward,  and  upward 
dimensions  are  determined.  Iteration  procedures  are  used  to  determine  the 
amount  of  X-ray  energy  contained  within  the  horizontal  radius  and  the  effect  of 
non-uniform  air  and  work  done  against  the  magnetic  field.  An  iteration  procedure 
is  also  used  to  find  the  downward  radius  in  a  non-uniform,  possibly  heaved 
atmosphere.  After  the  initial  dimensions  measured  from  the  burst  point  are 
found  the  location  and  dimensions  of  an  equivalent  spheroid  are  calculated. 

Next,  quantities  defining  the  fireball  rise  and  expansion  and  the 
initial  fireball  properties  (density,  scale  height,  mass,  and  volume)  are 
computed.  A  test  is  made  to  see  if  the  maximum  fireball  altitude  neglecting 
drag  forces  will  be  above  250  km.  If  it  is  not,  the  high-altitude  model  flag  is 
set  to  one  (low-altitude  burst)  and  control  is  returned  to  the  calling  routine. 

The  location  of  the  loss-cone  reference  point  (used  in  computing 
loss-cone  ionization)  is  determined  by  calling  routine  MAGFLD  to  determine  the 
location  of  the  point  on  the  field  line  through  the  burst  point  that  is  at  100  km 
altitude.  Also,  routine  EUXFIT  is  called  to  determine  the  fraction  of  ionizing 
UV  energy.  The  Lagrangian  coordinates  of  the  bottom  and  top  of  the  initial 
fireball  (used  in  the  fireball  wake  heave  modell  are  found  by  calling  routine 
LAGRAN. 

The  last  calculations  performed  are  to  determine  fireball  properties 
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when  the  fireball  reaches  apogee  and  at  the  magnetic  equilibrium  time.  The 
apogee  quantities  are  found  by  calling  routine  FBMHT  and  the  equilibrium 
properties  (mass  density)  by  calling  routine  FBITMD. 
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SUBROUTINE  FBMHT 


This  routine  determines  time-dependent  quantities  for  high-altitude 

bursts . 

The  routine  is  divided  into  two  major  computational  blocks;  one  for 
times  less  than  the  time  for  the  fireball  to  reach  apogee  (determined  in 
routine  FBMHI)  and  one  for  later  times.  An  input  logical  flag  (MODE)  is  used 
to  determine  apogee  fireball  quantities  (MODE  =  0)  or  fireball  quantities  at  a 
calculation  time  (MODE  =  1). 

l:or  MODE  =  0  or  for  times  less  than  the  apogee  time  for  MODE  =1, 
the  fireball  geometric  and  intrinsic  properties  are  determined  from  time- 
independent  quantities  calculated  in  routine  FBMHI  and  scaling  relations  defined 
for  several  time  periods  after  burst.  For  MODE  =  0  the  fireball  quantities 
(apogee  quantities)  are  stored  in  labeled  common  block  BREG. 

For  MODE  =  1  and  times  after  the  time  to  reach  apogee  the  location 
of  ionized  air  marker  parcels  are  determined  for  use  in  defining  ballistic  and 
diffusion  ionization.  Twenty  parcels  are  equally  spaced  along  the  longitudinal 
axis  of  the  apogee  fireball.  At  later  times  the  locations  of  the  highest  and 
lowest  parcels  in  the  burst  hemisphere  and  the  highest  and  lowest  parcels  in  the 
conjugate  hemisphere  (if  any  parcels  cross  the  magnetic  equator)  are  found.  The 
parcel  locations  are  used  to  define  coefficients  that  define  the  location  of  a 
parcel  in  the  apogee  fireball  in  terms  of  its  location  at  a  later  time. 

As  a  part  of  the  preliminary  calculations  the  time  (TMAX)  for  all  the 
parcels  to  fall  below  a  stopping  altitude  and  the  location  (XEMAX)  in  the  apogee 
fireball  separating  parcels  that  remain  in  the  burst  region  from  those  crossing 
the  equator  are  found.  Ionization  associated  with  parcels  that  are  above  the 
stopping  altitude  is  termed  ballistic  ionization.  Ionization  associated  with 
parcels  that  fall  below  the  stopping  altitude  is  called  diffusion  ionization. 

If  parcels  cross  the  equator,  parcel  location  scaling  coefficients  and  fireball 
dimensions  are  separately  defined  for  the  conjugate  region. 

The  parcel  locations  are  determined  by  calling  routine  BTRACE.  Inputs 
and  outputs  to  routine  BTRACE  arc  through  labeled  common  block  TRACE.  The 


138 


quantity  TPHS  is  the  time  for  the  parcel  to  fall  below  the  stopping  altitude, 
the  quantity  TPAP  is  the  time  for  the  parcel  to  reach  apogee  and  the  quantity 
CONJ  is  a  flag  that  is  set  to  a  negative  number  if  the  parcel  crosses  the 
equator. 
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SUBROUTINE  FBTINT 


This  routine  determines  the  theimal  emission  of  a  fireball  in  a 
specified  direction. 

The  thermal  emission  is  specified  in  terms  of  an  effective  fireball 
temperature  that  accounts  for  fireball  emissivity  and  temperature  along  a  path 
through  the  fireball  in  the  specified  direction.  For  low-altitude  bursts  the 
fireball  is  modeled  as  an  isothermal  region  and  the  fireball  temperature  and 
emissivity  at  the  fireball  edge  are  used  to  compute  the  effective  fireball 
temperature. 

For  high-altitude  bursts  both  the  emissivity  and  the  fireball 
temperature  are  functions  of  position  within  the  fireball.  The  effective 
fireball  temperature  is  computed  from 


where 


R 


(r)T(r)e 


-jj  a(x)dx 


dr 


a  *  incremental  absorption  within  the  fireball  (nepers) 
T  *  temperature  within  the  fireball  (K) 
r  *  the  distance  from  the  edge  of  the  fireball  along 
the  specified  path 

R  =  distance  through  the  fireball  along  the  specified 
path  . 

For  emissions  from  the  bottom  of  the  fireball  (paths  from  the  antenna  that  have  a 
positive  elevation  angle  in  the  fireball),  the  integration  is  performed  along 
the  path  from  the  receiver  according  to  the  above  equation.  Routine  FBEN  is 
called  to  obtain  a  and  T  as  a  function  of  position  within  the  fireball. 

For  emissions  from  the  top  of  the  fireball  a  direct  integration  of  the 
above  equation  could  require  an  excessive  number  of  integration  steps  since  most 
of  the  contribution  occurs  near  the  bottom  of  the  fireball.  In  order  to  limit 
the  number  of  steps,  an  integration  of  fireball  absorption  in  the  upward 
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direction  through  the  fireball  is  performed  first  and  emissivity  quantities 
at  each  integration  step  saved.  Then  an  integral  in  the  reverse  direction  is 
performed  using  the  integration  steps  determined  from  the  absorption  integral. 
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SUBROUTINE  FOOTPC 

This  routine  is  the  driver  routine  for  determining  fireball  region 
footprint  plots.  A  simplified  flow  chart  for  the  routine  is  shown  in  Figure 
15. 

Footprint  plots  are  not  prepared  for  low-altitude  bursts.  Absorption 
and  the  integral  of  electron  density  for  a  path  from  the  satellite  through  the 
fireball  center  are  given  as  output  along  with  the  location  of  the  corresponding 
ground  terminal. 

For  high-altitude  bursts  the  fireball  ionization  and  temperature  are 
fit  and  the  fit  parameters  stored  in  labeled  common  block  CHEMFB.  Then  a  scan 
loop  is  started  where  a  scan  consists  of  a  series  of  propagation  paths  from  the 
satellite  through  the  fireball.  The  scans  are  defined  by  points  along  the 
longitudinal  axis  of  the  fireball,  starting  at  the  fireball  bottom.  The 
propagation  paths  in  each  scan  are  chosen  to  be  in  a  plane  normal  to  the  fire¬ 
ball  longitudinal  axis  and  the  local  vertical  at  the  point  on  the  axis  defining 
the  scan.  An  attempt  is  made  to  define  five  paths  that  intersect  the  fireball 
(one  through  the  longitudinal  axis  and  two  on  each  side  of  the  axis). 

For  each  path  the  location  of  the  ground  terminal  is  computed  and 
propagation  quantities  determined  from  routine  FBINT.  Routine  FPLNK  is  called 
to  determine  link  quantities  including  the  effects  of  noise.  Two  additional 
paths  with  ambient  propagation  quantities  are  defined  for  each  scan  to  locate 
the  fireball  edges.  Then  routine  FPINT  is  called  to  obtain  parametric  values 
of  the  quantities  to  be  plotted  as  a  function  of  ground  terminal  location. 

After  all  the  scans  have  been  completed  routine  FPPLOT  is  called  to  prepare 
printer  page  plots. 
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Figure  15.  Flowchart  for  subroutine  FOOTPC. 


SUBROUTINE  FPCIR 

This  routine  determines  signal  processing  quantities  for  single  link 
circuits  used  for  footprint  plots. 

For  FM  systems  the  signal-to-noise  ratio  is  calculated  including  the 
effect  of  intermodulation  noise.  For  digital  systems  the  bit  energy-to-noise 
ratio  is  computed  and  depending  on  the  modulation  system  either  analytic  models 
are  used  to  determine  the  bit  error  rate  or  data  base  models  are  used  to 
determine  the  message  error  rate. 
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SUBROUTINE  FPINT 


This  routine  determines  parametric  values  of  footprint  plot 
quantities  as  a  function  of  ground  terminal  location. 

Calculations  are  performed  for  each  plot  quantity  and  stored  in  file 
LFIDO.  First,  the  maximum  value  of  the  plot  quantity  is  found  and  then 
parametric  values  between  the  minimum  value  (stored  in  a  data  statement)  and  the 
maximum  value  are  determined. 

Next,  the  location,  in  terms  of  distance  from  the  longitudinal  axis 
in  the  scan  plane,  of  paths  that  will  have  tne  parametric  values  of  the  plot 
quantity  are  determined  by  interpolation.  A  linear  interpolation  is  us- J  at  the 
edges  and  an  exponential  interpolation  elsewhere. 

After  the  location  of  the  paths  corresponding  to  the  parametric  values 
have  been  determined,  the  ground  terminal  location  for  each  path  is  computed  and 
the  results  stored  for  use  by  the  plot  routine. 


SUBROUTINE  FPLNK 


This  routine  determines  link  propagation  quantities  for  footprint 

plots. 

The  transmitter  and  receiver  terminal  locations  are  defined  from  the 
footprint  calculation  option  (downlink  for  KALCFP  =  1  and  uplink  for  KALCFP  -  2) 
and  the  terminal  systems  data  defined  from  the  systems  data  stored  for  index  1. 
Routines  GAIN  and  SLOBF,  are  used  to  determine  the  antenna  main  beam  and  side 
lobe  gains. 

If  a  system's  noise  temperature  was  not  specified  as  input,  the 
ambient  noise  level  at  the  receiver  is  computed.  Then  a  test  is  made  to  see  if 
noise  due  to  fireball  thermal  emissions  should  be  calculated.  The  fireball  noise 
is  always  calculated  for  p?rTs  in  the  first  scan  (see  description  of  routine 
FOOTPC)  and  the  maximum  temperature  saved.  For  paths  in  the  second  or  subsequent 
scans  the  fireball  noise  is  only  calculated  when  the  maximum  temperature  for 
the  first  scan  is  greater  than  250  degrees  Kelvin.  If  fireball  thermal  noise  is 
to  be  calculated,  routine  NIN'T  is  called  to  determine  a  normalized  integral 
over  the  antenna  pattern  and  routine  FBTINT  is  called  to  calculate  the  thermal 
emission  from  the  fireball. 

After  the  noise  has  been  determined  the  signal  power  noise  density, 
and  scintillation  quantities  are  stored  in  labeled  common  block  FTCIRC  for 
plotting  and  for  use  by  routine  FPCIR  to  determine  circuit  performance.  If  there 
are  angular  fluctuations  due  to  scintillation,  a  beam  spread  loss  is  computed 
and  the  signal  strength  reduced. 
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SUBROUriNE  FPPLOT 


This  routine  prepares  footprint  printer  page  plots  of  environment, 
propagation,  and  systems  performance  quantities. 

The  data  to  be  plotted  are  on  file  LFIDO  which  is  rewound  at  the 
start  of  each  plot.  A  plot  grid  is  determined  and  grid  points  located  at 
10  degree  intervals  in  latitude  and  longitude  for  an  elliptical  projection. 

Then,  for  each  scan  (see  description  of  routire  FOOTPC  for  description  of  a 
scan)  the  latitude-longitude  data  defining  contour  points  are  obtained  from 
fixe  LFIDO  and  plot  locations  and  plot  symbo' s  determined  for  each  contour 
point.  The  contour  points  always  start  at  che  minimum  value  (QTMIN  defined  in 
a  data  statement).  If  two  or  more  contour  points  are  at  the  same  plot  location, 
the  svmbol  for  the  largest  contour  value  is  used  to  define  the  point. 

After  the  plot  data  have  been  computed  the  printer  page  plot  is 
prepared  on  the  summary  output  file  (file  LFNO) . 

An  external  plot  file  can  be  prepared  by  setting  ID1  (dummy  control 
integer  in  the  OPTION  common  block)  to  an  integer  greater  than  zero.  This 
can  be  done  in  the  routine  INPUT  if  it  is  desired  to  control  preparation  of 
the  external  file  as  part  of  the  input  specification.  Data  are  written  on 
file  LFPLOT  which  is  set  to  7  in  the  control  routine.  The  following  describes 
the  data  written  on  the  plot  file. 

For  each  calculation  time  and  each  fireball  the  ordinate  and  abscissa 
scales  (applicable  for  each  plot  for  that  calculation  time  and  fireball)  and 
the  number  of  plots  are  written  on  the  file  as  follows: 


where 


TLAMAX,  TLAMIN,  T LOMAX,  TLOMIN,  NPLOT 

TLAMAX,  TLAMIN  =  maximum  and  minimum  values  of 

ordinate  scale  (north  latitude  in 
degrees) 

TLOMAX,  TLOMIN  =  maximum  and  minimum  values  of  abscissa 
scale  (east  longitude  in  degrees) 


NPLOT  =  number  of  plots  for  the  calculation 

time  and  fireball. 

Then  for  each  plot  two  records  are  written  on  the  file,  the  first  record 
describes  the  plot  and  the  grid  as  follows: 


where 


IQT,KFMFP,QTMINX,IXL,KK,((Y1(I),((X1(J=I), 
KK, 1*1 , IXL) 


IQT 


QTMINX  = 


IXL 

KK 

Y1(I),X!(J,I)  = 


plot  number 

1  =  absorption  (dB) 

_2 

2  =  electron  density  integral  (cm  ) 

3  *  reciprocal  of  correlation  time  (s"1) 

4  *  scintillation  index 

5  =  phase  standard  deviation  (radians) 

6  *  bit  error  rate  or  message  error 

rate  if  KFMFP=0  (s-1) 

*  reciprocal  of  signal -to-noise 
ratio  if  KFMFP=1 
minimum  value  of  contour. 

Contours  are  given  at  even  decades 
and  at  3  times  the  even  decade  values 
(eg,  10,30,100). 

number  of  latitudes  used  to  define  grid 
number  of  longitudes  used  to  define  grid 
coordinates  (in  terms  of  ordinate  and 
abscissa  scales)  of  points  showing  the 
location  of  longitude  traces.  The  points 
are  given  for  each  10  degrees  in  latitude 
and  longitude. 


I 
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The  second  record  describes  the  contours  as  follows: 


where 


NTnST,  N0C0N1,  MAXPTS,  ( (X2 ( I , J) ,Y2 (IJ) , 

1=1 , MAXPTS)  ,J--1 , NOCONI) 

NTEST  =  flag  to  determine  if  two  sets  of  contours 
will  be  prepared.* 

NCONI  =  number  of  contours 
MAXPTS  =  number  of  points  defining  contour 
X2(I ,J) ,Y2 (I , J)  =  coordinates  (in  terms  of  ordinate  and 
absscissa  scales)  of  points  defining 
contours. 


Note  that  if  NTEST=1,  the  next  data  on  the  file  are  the  data  for 
the  second  set  of  contours  (same  format  as  for  the  first  set  of  contours.) 
If  NTEST=0,  the  next  data  on  the  file  are  for  the  next  plot. 


*  As  described  in  Volume  1  (section  3. 2. 3. 4),  the  contour  plots  can  be  com¬ 
plicated  when  the  projection  of  the  longitudinal  axis  doubles  back  on  itself 
(two  points  on  the  projection  having  the  same  latitude  or  longitude).  For  these 
cases  two  separate  plots  are  prepared  for  the  external  plot  file  that  should  be 
plotted  on  the  same  grid.  The  actual  contours  would  be  continuous,  but  the 
separate  contours  illustrate  the  extent  of  effects  and  are  relatively  simple  to 
prepare.  If  two  sets  of  contours  are  prepared  NTEST  is  set  to  1  for  the  first 
set  and  to  0  for  the  second  set. 
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SUBROUTINE  GWPLME 

This  routine  determines  the  motion  of  a  high-altitude  fireball  after 
reaching  apogee  due  to  gravity  and  atmospheric  wind  forces.  A  simplified  flow 
chart  for  the  routine  is  shown  in  Figure  16. 

The  fireball  motion  is  controlled  by  gravity  and  wind  forces  acting 
near  the  altitude  of  peak  ionization.  While  this  altitude  is  calculated  in 
routine  FBEN,  it  is  simply  chosen  as  350  km  for  purposes  of  estimating 
fireball  motion.  The  location  of  the  point  along  the  central  field  line  in  the 
fireball  that  is  at  350  km  altitude  is  found  by  calling  routine  MAGFLD.  The 
number  of  time  steps  to  be  used  is  established  by  dividing  the  interval  between 
the  starting  time  and  the  final  time  into  the  minimum  jvimher  of  intervals  that 
are  not  over  1  hour  each.  The  date  (year,  JKMith,  day)  and  local  time  at  the 
starting  location  and  time  are  ?  "*■'  —  ,1  1  ->m  routines  ZTTOUT  and  DATE. 

Routine  DATE  is  called  to  le'"  «;e  »  time  at  the  time  of  interest  and 

then  routine  ZTTOdT  is  r  led  Xj  ue  i.e  the  local  time  at  the  location  of 
peak  ionization. 

For  each  time  step  the  local  time  at  the  location  of  peak  ionization 
is  found  and  the  east-west  and  north-south  velocities  due  to  winds  found 
from  routine  UVWD.  Then  the  velocities  due  to  gravity  are  computed  and  added 
to  the  wind  velocities.  The  effects  of  atmospheric  winds  will  not  be  included 
if  the  wind  option  (given  as  input  to  the  WF.SCOM  code)  is  exercised. 

After  the  velocities  are  determined,  the  location  of  the  peak 
ionization  point  at  the  end  of  the  time  step  is  found  by  assuming  a  constant 
velocity  during  the  time  interval.  rf  the  parameter  MODE  given  in  the  call 
statement  is  zero,  the  location  of  the  peak  ionization  at  the  end  of  the  last 
time  step  is  used  to  find  the  location  of  the  point  at  altitude  HSETL  (chosen  as 
150  km)  that  is  on  the  same  field  line  as  the  peak  ionization.  This  location 
is  used  as  the  debris  location. 

If  the  parameter  MODE  is  greater  than  zero,  the  location  and  velocity 
of  the  peak  ionization  point  is  found  at  the  end  of  the  last  time  step.  The 
reference  altitude  of  the  fireball  is  then  found  so  that  it  moves  from  the 
apogee  altitude  specified  in  routine  i  BMHT  to  the  altitude  of  peak  ionization 
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Figure  16.  Flowchart  for  subroutine  GWPLME. 


by  the  time  TMAX.  The  fireball  transverse  radius  and  tilt  angle  are  also 
modified  as  the  fireball  moves  across  field  lines. 
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SUBROUTINE  HPCHEM 

This  routine  partitions  the  total  energy  lost  by  loss  cone,  ion  leak,  and 

charge  exchange  particles  undergoing  inelastic  and  elastic  collisions.  The 

energy  is  partitioned  into  heavy  particle  thermal,  electron  thermal,  dissociation, 

ionization,  excitation,  and  radiation.  A  model  of  heavy  particle  prompt  chemistry 

4  2 

is  also  used  to  compute  the  change  in  number  densities  of  N2,  Oj,  N(  S),  N(  D), 

N*»  0*,  N+,  0*  and  electrons.  The  routine  is  essentially  the  same  as  routine 
HPCHEM  developed  for  R0SC0E. 
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SUBROUTINE  INITAL 

This  routine  determines  the  neutral  species  concentration  at  the 
end  of  phase  2  following  a  burst.  The  routine  is  similar  to  routine  INITAL 
developed  for  ROSCOE. 


First,  the  species  concentrations  after  the  prompt  energy  deposition 
are  obtained.  Then  a  test  is  made  to  see  if  the  energy  deposited  by  the 
burst  increases  the  gas  temperature.  The  temperature  is  computed  from 


where 

p  =  pressure  (dynes  cm 

k  =  Boltzmann's  constant  (1.38044  x  10"^  ergs  deg~* ) 
n  =  total  particle  density  including  charged  particles. 
If  the  temperature  has  changed  by  more  than  10  percent,  routine 
DRATF.  is  called  to  recompute  the  reaction  rate  coefficients. 
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SUBROUTINE  INPUT 


This  routine  reads  and  checks  input  data. 

On  the  first  call  to  routine  INPUT  routines  BLKDTA  and  IDEFLT  are 
called  to  obtain  constants,  parameters,  and  default  input  data.  Input  data  to 
be  read  in  for  the  current  problem  case  are  obtained  by  calling  routine  CARDIN 
which  stores  the  data  in  AFILE  when  the  ENCODE/DECODE  statements  are  used.  If 
these  statements  cannot  be  used,  the  input  data  can  be  stored  in  file  LFIDO  as 
indicated  by  comment  statements  in  routines  CARDIN  and  INPUT. 

For  each  numbered  input  block  the  input  data  read  in  on  the  numbered 
input  card  are  stored  in  common.  If  additional  input  data  are  to  be  read  in  as 
part  of  the  input  block,  routine  CARDIN  is  called  and  data  obtained  from  input 
cards  with  a  letter  in  column  1  as  described  in  Volume  1  (Section  3.1.1).  Tests 
are  made  to  insure  that  sufficient  data  are  given  and  that  the  data  are  in  the 
proper  order.  In  some  cases  routine  TEST  is  called  to  check  the  magnitude  of  an 
input  quantity.  If  errors  are  found,  error  messages  are  prepared. 

If  an  orbiting  satellite  is  specified  in  Input  Block  3,  routine  ORBIT 
is  called  to  obtain  orbit  parameters  that  can  be  used  to  determine  the  satellite 
position  at  each  calculation  time.  In  processing  Input  Block  6  data,  routine 
SPCTRM  is  called  to  obtain  weapon  device  data  and  radiation  spectrums.  Routines 
W0G1 ,  W0N1 ,  and  W0X1  are  called  to  determine  parametric  energy  deposition  data 
for  later  use  in  computing  energy  deposition  at  specific  locations. 

If,  after  reading  the  input  data  for  a  problem  case,  there  are  no 
errors  detected,  routine  INPUTP  is  called  to  perform  preliminary  processing  of 
the  data.  Calculations  are  stopped  by  calling  library  routine  EXIT  when  an  input 

card  with  a  twelve  in  columns  1  and  2  is  given.  If  errors  are  detected  in  the 

input  data  for  a  problem  case,  the  rest  of  the  input  cards  for  that  and  sub¬ 
sequent  problem  cases  are  checked  for  additional  errors  and  library  routine  EXIT 
called  when  the  input  card  with  a  twelve  in  columns  1  and  2  is  read. 

Additional  control  options  can  be  read  in  as  part  of  Input  Block  1  and 

stored  as  one  of  the  dummy  variables  in  OPTION  common  (see  Table  5).  This  may 
be  useful,  for  example,  to  control  preparation  of  user  plot  file. 


SUBROUTINE  INPUTP 

This  routine  performs  preliminary  processing  of  input  data  and 
prepares  summary  output  describing  input  data. 

The  major  data  processing  performed  is  geometry  coordinate  conversion. 
Routine  CORDTF  is  called  to  convert  input  coordinates  to  an  earth-centered  vector 
coordinate  system.  Tangent  plane  coordinates  and  geographic  coordinates  are 
also  determined  for  summary  output. 

Data  processing  is  performed  when  FM  modulation  is  specified.  Routines 
BB  and  RRMSTP  are  called  to  obtain  the  maximum  baseband  frequency  and  rms  and 
peak  frequency  deviations.  Calculations  are  also  made  to  convert  system  noise 
figure  to  noise  temperature. 


SUBROUTINE  MANPTS 


This  routine  determines  mandatory  integration  points  related  to  beta 
tubes,  low-altitude  ionization  sources,  and  heavy  particle  ionization  sources. 

The  low-altitude  ionization  source  points  are  located  near  the  closest 
point  of  approach  of  the  path  and  the  ionization  source.  The  beta  particle 
mandatory  points  are  located  at  the  points  where  the  ray  path  intersects  the 
beta  ionization  region  where  the  region  is  taken  as  large  or  larger  than  the 
actual  region.  If  detailed  propagation  output  is  requested  (input  option)  the 
beta  region  intersection  altitudes  are  written  out  on  a  detailed  output  file 
(file  LFDO) .  Note  that  because  the  beta  ionization  region  used  to  determine 
the  mandatory  points  is  chosen  conservatively  (larger  than  the  actual  region) 
the  actual  beta  particle  ionization  may  occur  over  a  smaller  altitude  interval 
along  the  ray  path  than  indicated  by  the  mandatory  altitudes. 

The  ion-leak,  loss-cone,  and  charge-exchange  heavy  particles 
produce  ionization  within  geomagnetic  tubes  similar  to  the  beta  tube.  It  would 
be  desirable  to  calculate  mandatory  points  for  the  heavy  particle  ionization 
regions  to  insure  calculations  where  effects  are  large  and  to  eliminate 
calculations  in  regions  where  effects  would  be  minimal.  A  special  problem 
arises  with  the  heavy  particle  mandatory  points  since  they  refer  to  the  time 
when  the  energy  is  deposited  (burst  times)  and  not  to  the  calculation  time.  Thus, 
the  location  of  the  air  ionized  by  ion-leak,  loss-cone,  and  charge-exchange 
particles  can  move  due  to  heave  such  that  it  is  not  between  mandatory  points 
determined  at  burst  time.  However,  since  only  vertical  air  motion  is  allowed 
in  the  current  heave  model,  heavy  particle  ionization  calculations  can  be 
minimized  by  only  making  calculations  for  propagation  paths  which  are  in  or 
above  the  initial  ionization  tubes.  The  location  of  the  propagation  path 
relative  to  the  heavy  particle  ionization  tubes  is  determined  by  comparing 
their  geomagnetic  coordinates  (latitude  and  longitude).  If  the  portion  of  the 
propagation  path  above  100  km  is  in  or  above  the  ionization  tubes,  logical 
calculation  flags  are  set  to  unity  (ionization  calculations  to  be  performed). 
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SUBROUTINE  NINT 

This  routine  determines  a  normalized  integral  over  the  antenna  pattern 
used  in  computing  the  receiver  antenna  noise  temperature  due  to  fireball  thermal 
radiation. 

The  noise  temperature  is  computed  assuming  that  the  radiation  is 
independent  of  position  within  the  antenna  beam.  Then  the  integral  of  receiver 
antenna  gain  over  the  fireball  can  be  computed  separately  from  fireball  temper- 
atir .»  and  emissivity  calculations. 

The  method  of  evaluation  depends  upon  the  location  and  size  of  the 
fireball  with  respect  to  the  antenna  pattern.  If  the  fireball  is  entirely  in 
the  sidelobe,  or  if  it  is  very  small  relative  to  the  main  lobe,  the  integral 
is  approximated  as 

w  0 

4  it  FB 


angle  off  beam  maximum  to  fireball  center 

solid  angle  subtended  by  fireball 
receiver  gain. 

If  the  fireball  extends  into  the  main  lobe  and  is  not  small  relative  to  the  main 
lobe,*  a  numerical  integration  procedure  is  employed.  The  main  lobe  is 
partitioned  into  zones  by  dividing  the  angular  extent  off  the  beam  maximum  into 
four  equal  rings  and  cutting  the  360-degree  extent  of  each  ring  around  the  beam 
maximum  into  six  equal  sections,  except  the  zone  which  is  circular.  This  yields 
a  total  of  19  zones,  A  ray  through  the  center  of  each  zone  is  checked  to 


♦The  main  lobe  is  considered  to  extend  one  receiver  antenna  beamwidth  off  the 
beam  maximum,  at  which  point,  the  (sin  x/x)^  function  is  down  about  18  dB  for  a 
mechanically  steerable  or  phased  array  type  antenna. 
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determine  if  it  intersects  the  fireball.  If  it  does,  the  fireball  is  assumed 
to  fill  the  zone.  The  resulting  numerical  approximation  to  the  desired  integrals 
may  be  written  as 


where 


^  d«  »  £  Gr(*j) 


4tt 


<p j  =  angle  off  beam  maximum  to  the  ray  through  the 
center  of  the  jth  zone  which  intersects  the 
fireball 

flj  =  solid  angle  covered  by  the  jth  zone. 


In  performing  the  integration,  tests  are  performed  to  make  sure  that  the 
integrals  do  not  exceed  their  theoretical  upper  limits  and  that  the  total  solid 
angle  included  in  the  integration  does  not  exceed  the  solid  angle  subtended  by 
the  fireball.  Once  all  zones  have  been  checked,  the  contribution  of  the 
remaining  portion  of  the  fireball,  if  any,  is  computed  by  assuming  that  the 
remaining  portion  is  in  the  side  lobe.  This  integration  procedure  allows  the 
several  fireball  geometries  defined  by  the  phenomenology  model  to  be  handled  in 
the  same  basic  manner,  with  the  only  major  differences  occurring  in  the  relations 
needed  to  determine  whether  a  ray  through  the  center  of  a  zone  intersects  the 
fireball.  This  determination  is  performed  by  the  intersection  subroutine  RINTER. 
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SUBROUTINE  PHENOM 

This  routine  is  a  driver  routine  for  the  phenomenology  models,  A 
simplified  flow  diagram  for  the  routine  is  shown  in  Figure  17. 

The  PHENOM  routine  is  called  for  a  specified  time  at  which  fireball, 
debris,  and  dust  region  parameters  are  desired.  A  test  is  made  to  determine 
whether  any  burst  times  are  between  the  current  calculation  time  and  the  last 
calculation  time  for  which  routine  PHENOM  was  called.  If  there  are  bursts  in 
the  last  time  interval,  time -independent  scaling  for  each  such  burst  is 
performed. 

The  first  step  in  the  time-independent  block  is  to  determine  the 
atmospheric  properties  including  those  of  a  dipole  geomagnetic  field  at  the  burst 
point.  If  the  modeling  option  parameter  MAMB  is  zero,  the  atmospheric  properties 
are  determined  for  conditions  applicable  at  the  origin.  If  MAMB  is  greater  than 
zero,  atmospheric  quantities  are  determined  for  the  particular  burst  location 
and  time.  Next,  the  time  independent  fireball  quantities  are  found.  For  burst 
altitudes  less  than  or  equal  to  75  km  the  low-altitude  time- independent  fire¬ 
ball  model  (routine  FBMLI),  the  radiation  output  model  (routine  KADOUT) ,  and 
the  dust  model  (routine  DUSTMI)  are  called,  Roulire  RADOUT  provides  the  total 
photon  emission  (assumed  to  occur  instantaneously)  causing  photodissociation  of 
ozone.  For  burst  altitudes  above  75-km  routine  FBMHI  is  called  to  determine  the 
high-altitude  time  independent  quantities.  However,  if  the  fireball  does  not 
rise  above  250  km,  a  logical  flag  LMCF  is  set  to  1  and  the  low-altitude  model  is 
used.  Thus,  the  "low-altitude"  model  is  used  for  all  bursts  below  75  km  and 
for  bursts  above  75  km  where  the  fireball  c'oes  not  rise  above  250  km. 

Three  labeled  common  blocks  are  used  to  store  quantities  calculated 
in  the  time-independent  section  of  routine  PHENOM.  Quantities  used  only  by 
routines  called  from  routine  PHENOM  are  stored  in  the  labeled  common  block  PHEN. 
Quantities  used  by  routines  other  than  those  called  by  routine  PHENOM  are  stored 
in  labeled  common  blocks  BREG  and  DCREG.  The  quantities  stored  in  block  OCR EG 
describe  dust  clouds. 

After  the  time-independent  quantities  have  been  calculated  for  bursts 
occurring  since  the  last  calculation  time,  the  time-independent  test  on  LMCF  is 
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Figure  17.  Flowchart  for  subroutine  PHENOM. 
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made  and  either  the  low-  or  high-altitude  models  used.  Routine  RADOUT  is  called 
in  the  time-dependent  section  to  d,cerminc  the  time-dependent  photon  emission 
causing  photodetachment  and  photodis  ociation  of  negative  ions.  Routines  DBMLT 
anu  D8MHT  provide  time-dependent  debris  region  scaling  for  low-  and  high-altitude 
bursts  respectively.  Routine  GWPLME  determines  the  effects  of  gravity  and 
atmospheric  winds  on  the  late-time  fireball  plume. 

* .  ae-independent  data  are  stored  in  two  labeled  common  blocks: 
block  FBREG  for  fireball  quantities  and  block  DBREG  for  debris  quantities. 

The  fireball  quantities  arc  stored  in  terms  of  the  fireball  index  (NEB)  which 
is  the  same  as  the  burst  index  (one  fireball  per  burst).  The  debris  region 
properties  are  indexed  with  a  separate  debris  index  (NDB) .  Up  to  four  debris 
regions  can  be  used  to  represent  the  debris  distributions  for  a  burst.  The 
number  of  debris  regions  is  redetermined  at  each  calculation  time  and  may 
change  for  a  given  burst.  Logical  indexes  are  stored  that  define  the  fireball, 
burst,  and  weapon  device  type  for  each  debris  region. 

After  the  time-independent  quantities  have  been  computed  for  each 
fireball  and  debris  region,  the  input  parameter  MWIND  is  checked.  If  MW I ND  =  1, 
routine  DBSTAM  is  called  to  determine  the  effect  of  atmospheric  winds  on  the 
late-time  debris  region  location  and  geometry. 
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SUBROUTINE  PIONF 


This  routine  determines  the  species  densities,  pressure,  and  temper¬ 
ature  at  a  specified  location  above  100  km  following  the  energy  deposition  due 
to  radiation  from  a  specified  burst.  A  flowchart  for  the  routine  is  shown  in 
Figure  18. 

The  effective  slant  range  from  the  burst  to  the  deposition  point  is 
computed  first.  If  the  specified  deposition  point  is  within  the  initial  fire¬ 
ball,  its  effective  location  is  assumed  to  be  at  the  edge  of  the  fireball  in 
the  same  direction  from  the  burst  point. 

Next,  the  UV  source  strength  for  each  energy  group  is  computed  in  the 
direction  of  the  deposition  point.  Parameters  describing  the  extent  of  the  UV 
"burnout"  region  are  obtained  from  BREG  common.  The  location  and  atmospheric 
properties  of  the  edge  of  the  burnout  region  in  the  direction  of  the  deposition 
point  are  obtained  for  use  in  the  air  mass-penetrated  calculation  for  attenu¬ 
ation  of  the  UV  flux. 

Next  the  properties  of  the  air  at  the  fireball  edge  are  found  from 
subroutine  ATMOSU  or  ATMOSF,  and  the  calculation  of  mass  penetrated  and  number- 
density  line  integrals  from  the  fireball  edge  to  the  deposition  point  is  begun. 
If  any  point  on  the  path  is  at  or  below  100  km,  the  effects  of  heave  on  the  air 
density  along  the  path  are  neglected  and  the  mass  penetrated  is  determined  with 
function  PMASS.  If  the  path  is  entirely  above  100  km,  the  mass  penetrated  is 
obtained  by  dividing  the  path  into  segments,  assuming  the  mass  density  varies 
exponentially  with  slant  range.  In  addition,  neutral-particl f  composition 
terms  are  determined  to  convert  mass  penetrated  to  the  line  integrals  of 
individual  neutral  particles. 

For  computing  the  mass  penetrated  in  the  heaved  air,  the  path  is 
divided  into  N  segments,  where  N  is  defined  as  follows: 


where 


N  =  Minimum[K,  Int(S/DMAX)  +  1] 

K  =  maximum  number  of  path  segments 


163 


EFFECT IVt  SLANT  RANGE  TO  DEPOSITION  POINT 


PROPERTIES  OF  AIR  AT  FB  EDGE 


path  eno  pomr 

^ BELOW  100  «A  . 


'point  on  path 

LOWER  THAN  ENP 
.  POINTS 


LOWEST  POINT 
BELOW  IOO  KN 


164 


=  N1  if  lowest  point  is  end  point 
=  N2  if  lowest  point  is  between  end  points 

Int  t X)  «  largest  integer  smaller  than  X 

S  =  path  length  (km) 

DMAX  =  minimum  length  of  path  segments  that  can  be  further 
divided  (km). 

The  parameter  values  used  are  N1  =  3,  N2  =  3,  and  DMAX  =  100  k^.  The  properties 
of  the  air  at  the  end  points  of  each  path  segment  are  obtained  from  subroutine 
ATMOSF,  and  the  mass  penetrated  along  each  segment  is  determined  assuming  the 
mass  density  varies  exponentially  with  slant  range.  If  the  initial  path  from 
the  fireball  edge  to  the  deposition  point  has  a  lowest  point,  then  each  seg¬ 
ment  may  be  further  subdivided  as  described  above. 

The  energy  deposited  by  heavy  particles  is  found  by  calling  routine 
DLRRTS.  Routine  HPCHEM  is  called  to  determine  the  change  in  species  concentra¬ 
tions  and  air  properties  caused  by  the  heavy  particle  energy  deposition.  The 
X-ray  energy  deposition  coefficient  is  found  by  calling  routine  WOXP  and  used 
to  determine  the  energy  depositied  by  X-rays. 

The  UV  flux  in  five  energy  groups  is  determined.  The  species  line 
integrals  arc  scaled  from  the  mass  penetrated  and  the  composition  terms.  The 
UV  energy  available  for  hydro  is  added  to  the  energy  deposition  due  to  X-rays  to 
compute  a  modified  pressure.  Then  routine  PCHEM  is  called  to  determine  the 
species  densities  and  air  properties  at  the  end  of  the  prompt  energy  deposition 
phase. 


i 
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SUBROUTINE  PROPEN 

This  routine  is  the  driver  routine  for  the  propagation  module.  The 
module  determines  environmental  quantities  and  propagation  effects  along  a 
straight-line  path  (ray  path)  between  two  specified  locations.  A  flowchart 
for  the  routine  is  shown  in  Figure  19. 

The  propagation  module  can  be  called  with  one  of  two  calculation 
options  determined  by  the  logical  option  MODE.  For  MODE  =  1  only  absorption 
between  the  terminals  is  calculated.  This  option  is  used  when  calling  routine 
PROPEN  from  the  noise  module.  For  MODE  =  2  all  of  the  propagation  effects 
requested  as  part  of  the  input  specification  are  calculated.  This  option 
is  used  for  the  main  ray  path  specified  as  input. 

First,  preliminary  geometry  calculations  are  made  and  integral 
quantities  initialized.  The  geometry  calculations  include  determination  of 
minimum  and  maximum  mandatory  slant  ranges  from  the  transmitter  terminal.  These 

arc  chosen  as  50  and  101  km  for  MODE  =  1  and  as  50  and  250  km  for  MODE  =  2  if  the 

ray  path  extends  over  these  altitudes.  If  not,  the  ray  path  terminals  are  used 

for  the  mandatory  points.  An  F-region  altitude  (IIFTEST)  is  chosen  for  use  in 

editing  bursts  as  F-region  ionization  sources.  This  altitude  is  chosen  as 
200  km  or  the  minimum  altitude  on  the  path  if  the  ray  path  is  above  200  km.  If 
the  ray  path  is  entirely  below  200  km,  HFTEST  is  chosen  as  the  highest  point  on 
the  ray  path. 

Next,  routine  AMBABS  is  called  to  determine  the  attenuation  due  to 
water  vapor  and  molecular  oxygen  in  the  ambient  atmosphere.  This  calculation 
is  not  made  for  MODE  =  1. 

Next,  mandatory  altitudes  (altitudes  at  which  environment 
calculations  must  be  made)  for  beta  tubes,  low-altitude  ionization  sources, 
and  heavy  particle  ionization  sources  are  obtained  by  calling  routine  MANPTS. 

After  the  mandatory  points  are  determined  a  test  is  made  on  HFTEST. 

If  HFTEST  is  greater  than  100  km,  routine  CMFEDT  is  called  to  select  bursts  to 
be  used  as  F-region  ionization  sources.  Routine  FBINT  is  called  to  determine 
fireball  intersections  and  propagation  effects  and  routine  DSTINT  is  called  to 
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Figure  19.  Flowchart  for  subroutine  PROPEN. 
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determine  dust  cloud  intersections  and  propagation  effects. 

Routine  START  prepares  an  ordered  set  of  mandatory  points  from  those 
determined  above.  The  mandatory  points  are  used  in  routine  PTHINT  called  to 
obtain  propagation  effects  along  the  path  between  terminal  points.  If 
detailed  output  was  temporarily  stored  in  file  LFIDO,  the  output  is  trans¬ 
ferred  to  file  LFDO. 
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SUBROUTINE  PSDM 

This  routine  is  the  driver  routine  for  the  post-stabilization  debris 
module  for  a  specified  debris  region.  Figure  20  shows  a  flowchart  for  the 
routine. 

First  the  number  and  location  of  marker  particle  altitudes  and  the 
number  and  times  of  evaluation  times  are  determined.  The  initial  location  of 
the  marker  particles  are  established.  The  initial  location  is  taken  as  the 
latitude  and  longitude  of  the  debris  center  at  stabilization  time  and  the 
altitudes  determined  above. 

Next,  a  time  loop  and  an  altitude  loop  are  started.  For  each  particle 
the  particle  velocity  is  determined  by  calling  routine  UVWD. 

After  vhe  marker  particle  velocities  have  been  determined,  the 
velocities  are  averaged  over  a  mixing  length  by  calling  routine  UVWAV.  Then 
the  velocities  are  used  to  determine  the  location  of  the  marker  particles  at  the 
end  of  the  time  interval. 

At  the  end  of  the  last  time  interval  (time  equal  to  input  calculation 
time),  routine  ARRLIM  is  called  to  determine  the  debris  location  and  geometry 
from  the  marker  particle  locations. 
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Figure  20.  Flowchart  for  subroutine  PSDM. 
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SUBROUTINE  PTH I NT 


This  routine  determines  ionization  and  propagation  effects  along  a 
specified  ray  path.  A  flowchart  for  the  routine  is  shown  in  Figure  21. 

Calculations  are  started  at  the  first  mandatory  point  (mandatory  point 
closest  to  transmitter  terminal)  and  integral  quantities  are  determined  for  the 
region  between  the  first  mandatory  point  and  the  transmitter.  At  each  calcula¬ 
tion  point  routine  ELDEN  is  called  to  determine  ionization  and  collision 
frequency  terms.  Then  the  incremental  absorption  and  incremental  phase  term 

N  where 
er 


O) 


are  calculated.  For  altitudes  below  100  km,  a  nominal  step  size  AS 
corresponding  to  a  2-1/2-km  altitude  increment  is  chosen  (but  not  more  than 
10  km),  with  a  minimum  step  size  AS  ^  equal  to  half  that  value  and  a 
maximum  step  size  AS  equal  to  twice  that  value.  For  altitudes  above  100  km, 
the  minimum  step  size  is  increased  to  2  AS  and  the  maximum  step  size  is 
increased  to  8  AS.  The  step  size  starts  out  at  AS  .  and  is  doubled  (subject 
to  ASmax)  each  time  the  phase  effects  term  and  incremental  absorption  vary  by 
less  than  a  factor  or  two.  If  the  electron  density  is  decreasing  with  position, 
the  maximum  step  size  is  allowed  to  continue  to  increase. 

The  phase  effects  term  and  the  incremental  absorption  are  assumed 
to  vary  linearly  between  two  field  points  if  the  change  in  magnitude  is  less 
than  a  factor  of  2;  otherwise  an  exponential  variation  is  assumed. 

The  integration  from  the  first  mandatory  point  to  the  transmitter  is 
stopped  when  all  of  the  following  checks  are  satisfied:  (1)  the  incremental 
absorption  is  decreasing;  (2)  the  incremental  absorption  multiplied  by  the 
remaining  path  length  is  less  than  0.1  dB;  and,  for  MODE  =  2  only,  (3)  the  phase 
effect  term  is  decreasing;  (4)  the  phase  effect  term  multiplied  by  the  remaining 
path  length  is  less  than  0.1  times  the  current  phase  integral. 


171 


Figure  21.  Flowchart  for  subroutine  PTHINT. 
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The  remaining  path  length  in  the  above  is  computed  as  the  minimum  of  the 
remaining  path  length  and  the  exponential  absorption  and  phase  effect  scale  lengths 

After  the  integral  quantities  between  the  first  mandatory  point  and 
the  transmitter  have  been  calculated  the  integration  between  the  first  mandatory 
point  and  the  receiver  terminal  is  started.  The  procedure  is  the  same  as 
previously  described  except  that  the  integration  is  not  allowed  to  stop  as  long 
as  there  are  any  mandatory  points  left.  At  each  mandatory  point  the  step  size 
is  set  to  the  minimum  step  size.  Tests  are  made  at  each  mandatory  point  to  see 
if  the  mandatory  point  is  associated  with  a  fireball.  For  low-altitude  fireballs 
and  the  fireball  from  the  last  high-altitude  burst,  the  region  along  the  path 
occupied  by  the  fireball  is  skipped  since  all  ionization  sources  are  included 
in  the  fireball  calculation.  If  there  is  more  than  one  high-altitude  burst, 
the  ionization  and  propagation  effects  for  the  region  occupied  by  all  but  the 
last  fireball  are  calculated  as  if  the  earlier  fireballs  were  not  there  since 
prompt  and  delayed  radiation  sources  are  not  included  in  the  region  calculations. 
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SUBROUTINE  RAOOUT 

This  routine  calculates  either  of  two  types  of  output  for  a  specified 
burst.  The  first  type  of  output  approximates  the  radiated  power  in  each  of  six 
frequency  bands  covering  IR  through  far  UV  at  a  specified  calculation  time.  The 
second  type  of  output  approximates  the  radiated  energy  resulting  from  the 
integration  of  power  in  the  sixth  frequency  band  over  a  time  span  from  10  6  s 
to  approximately  ten  times  thermal  maximum,,  The  routine  is  a  modified  version 
of  routine  RADOUT  developed  for  R0SC0E. 

Calculations  are  made  for  bursts  below  100  km;  however,  bursts 
between  50  and  100  km  are  treated  as  if  they  were  detonated  at  50  km.  The  power 
and  energy  outputs  are  determined  from  empirical  fits  to  RADFLO  data.  The 
power-time  curve  is  divided  into  six  time  segments  and  empirical  fits  used  for 
each  time  segment.  The  radiated  power  and  •  v  arc  modified  if  there  has 
been  a  radiation  merge  to  account  for  the  eftoc*  -r  reduced  air  density  caused 
by  the  first  fireball. 
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SUBROUTINE  RINTER 


This  routine  computes  the  intersection  points  between  a  ray  path  and 
a  fireball,  debris,  or  beta  tube  region.  The  region  shapes  included  are  the 
right-circular  or  skewed  spheroid,  cones,  magnetically  defined  tubes,  and  the 
torus.  The  routine  is  a  modified  version  of  routines  INTERS  developed  for 
WOE  code. 

The  intersection  with  the  spheroid  or  cone  is  obtained  analytically. 

•> 

Let  S  be  the  vector  from  the  transmitter  to  a  point  a  distance  s  up  the 
ray  path: 


S  =  sS  (1) 

A  "►  -> 

where  S  is  the  unit  vector  in  the  direction  of  S  .  Let  F  be  a  vector  from 
the  transmitter  to  the  center  of  the  spheroid  or  cone  and  let  T  a  a  vector 
from  the  center  of  the  spheroid  or  cone  to  a  point  on  its  symmetry  axis  a 
distance  t  from  the  center: 


T  = 


tT 


(2) 


where  T  is  again  a  unit  vector. 

The  vector  from  point  s  to  point  t  is 

D(s,t)  =  F  ♦  T(t)  -  S(s)  .  (3) 

The  distance  from  the  point  s  to  point  t  is  the  absolute  value  of  D  : 

D2(s,t)  =  |D(s,t)|2  (4) 

=  j  F J 2  +  t2  +  s2  +  2t (F  •  T) 

4  A  A  A 

-  2s(F  •  S)  -  2ts (T  •  S)  . 

For  any  distance  s  up  the  ray  path,  the  distance  from  s  to  the 

closest  point  on  the  symmetry  axis  of  the  spheroid  or  cone  can  be  found  by 
2 

minimizing  D  (s,t)  with  respect  to  t: 
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0  =  Jt  I)2^s*t) 


A  A 


=  2t  +  2(F  •  T)  -  2s(T  •  5) 


Thus 


A  A 


t  =  (T  •  S)s  -  F  •  T  =  as  +  b 


(3) 


Substituting  t  back  into  Equation  4  gives  the  shortest  distance  from  s  to 
the  axis: 


D2  =  As2  +  Bs  ♦  C 


(6) 


where 


A  A 


A  A 


A  =  1  +  a  -  2a (T  •  S) 

B  =  2 [ab  -  F  •  S  ♦  a(F  .  T)  -  b(T  •  S) 
C  =  | F | 2  ♦  b2  +  2b (F  •  T) 


>]  ■ 


The  distance  R  from  the  surface  to  the  axis  of  a  right-circular  spheroid 
or  cone  is  just: 


R2=At2+Bt+C 

q  q  q 


(7) 


where  for  the  spheroid 


A  =  -R^/R2l  . 


B  =  0 

q 


Cq  = 


R^  =  semi-axis  along  symmetry  axis. 


R^,  =  semi-axis  transverse  to  symmetry  axis. 
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and  for  the  cone 


A 

*  2 
tan 

9  . 

q 

B 

q 

2  R 

0 

2 

tan  6  , 

c 

q 

R 

0 

> 

R 

cone 

radius  at  center 

0 

0  = 

cone 

half-angle 

Writing  R  as  a  function  of  s  through  Equation  5  and  subtracting  it  from 
2 

P  gives  a  quadratic  equation  in  s: 

2  2  *  2  »  t 

Dt  -  R  +  A  s  ♦  B  s  ♦  C  =  0  ,  (8) 


where 


a'  =  A  -  a2A  , 

q 

B*  =  B  -  2  abA  -  aB  , 

q  q 

c'  =  C  -  b2A  -  bB  -  C 

q  q  q 

The  solutions  for  the  intersection  points  are  then 

S  =  -(b'/2a')±  /(B'/2A')2  -  c'/a'  ,  (9) 

where  an  intersection  exists  only  if  the  radical  is  real.  Two  special  cases  may 
arise.  If  A*  =  0,  which  corresponds  to  the  ray  path  being  parallel  to  part  of 
a  cone  surface,  then  there  is  only  one  intersection  given  by 

S  =  -c'/B'  (10) 

j 
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If  B*  is  positive  the  intersection  is  an  exit  point,  and  if  is  negative, 
the  intersection  is  an  entry  point.  If  B*  and  Af  are  zero,  the  ray  path  is 
parallel  to  a  cylinder  and  is  wholly  within  or  without  depending  on  whether  Cf 
is  negative  or  positive,  respectively. 

To  calculate  the  intersection  of  the  ray  path  with  a  skewed 
spheroid  or  cone,  Equation  5  must  be  replaced  with  a  relation  between  s  and 
t  such  that  point  s  is  horizontal  to  t  .  This  is  given  by 


s 


s 


p 


+ 


COS  T 

sin  E 

s 


(11) 


where 


s 


P 


slant  range  up  path  to  horizontal  plane 
through  region  center 


T  =  region  symmetry  axis  tilt  angle  from  vertical 


Eg  =  angle  between  ray  path  and  horizontal  plane 
through  region  center. 


Thus  a  and  b 


in  Equation  5  are  replaced  by 
sin  E 


_ s 

COS  T 


> 


b 


sin  E 

s 

s  - 

p  cos  T 


(12) 


and  the  rest  of  the  calculations  are  the  same  except  that  Equation  7  represents 
the  corresponding  skew  figures  (R  is  horizontal). 

The  equation  for  the  intersection  of  a  ray  path  and  a  torus  is  a 
quartic  expression  in  the  path  variable  and  is  tedious  to  solve.  It  is  simple, 
however,  to  determine  whether  or  not  a  given  point  is  inside  the  torus.  This 
allows  a  fast  numerical  procedure  for  finding  the  intersection  points.  The 
numerical  iteration  for  the  torus  intersection  points  is  carried  out  in 
Subroutine  TOROID  which  is  called  from  RINTER. 
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Intersections  with  magnetic  field  line  defined  tubes  are  also 
obtained  numerically  via  iteration  techniques.  The  function  TUBE  calculates 
the  intersection  points.  If  any  of  the  regions  are  truncated,  the  routine 
calculates  the  intersection  of  the  ray  path  with  the  bottom  and  top  altitudes 
and  adjusts  the  previously  obtained  intersection  points  if  necessary. 

Intersections  with  beta  tubes  are  found  by  first  determining  the 
closest  point  of  approach  of  the  ray  path  vector  and  a  vector  located  at  80  km 
on  the  geomagnetic  field  line  through  the  center  of  the  debris  region  having 
the  direction  of  the  geomagnetic  field.  The  closest  point  of  approach  is  used 
as  the  starting  point  in  routine  TUBE  to  determine  intersection  points.  First, 
routine  TUBE  is  called  with  a  single  point  calculation  option  to  determine 
whether  the  closest  point  of  approach  is  within  the  tube.  If  it  is  within  the 
tube  in  either  the  debris  region  or  the  conjugate  of  the  debris  region, 
routine  TUBE  is  called  again  with  an  option  to  allow  iterating  to  find  the 
ray  path  and  tube  boundary  intersections. 

Use  of  the  closest  point  of  approach  of  the  ray  path  with  a 
straight  line  approximation  to  the  geomagnetic  field  to  obtain  the  starting 
point  and  initial  test  of  whether  any  intersection  exists  is  a  reasonable 
approximation  for  beta  particle  ionization  where  the  interest  is  in  a  relatively 
small  altitude  region.  However,  for  field-aligned  fireballs  the  intersection 
of  interest  can  occur  over  a  large  altitude  region.  For  this  case,  the  region 
(field-aligned  fireball  region)  is  broken  up  into  five  altitude  segments  and 
tests  made  for  intersections  (in  both  the  local  and  conjugate  region)  for  each 
segment. 


179 


SUBROUTINE  SATC 

This  routine  is  the  driver  routine  for  satellite  communication 
calculations.  A  simplified  flowchart  for  the  routine  is  shown  in  Figure  22. 

First,  the  ground  terminal  locations  are  determined  from  either 
locations  given  as  input  or  from  routine  TORBIT  if  ground  terminal  motion  is 
specified.  For  KALCS=1  or  2  (multiple  link  circuits)  routine  SLINK  is  called 
to  determine  whether  propagation  and  noise  calculations  have  been  previously 
made  for  the  link.  If  they  have,  the  calculation  flags  LFLAG  and  NFLAG  are  set 
to  1.  For  KALCS-2  (down  link  broadcast  option)  each  link  is  different  and 
circuit  calculations  are  made  for  each  link. 

Next  path  geometry  quantities  are  defined  for  the  link  and  if 
MAMB  (input  option)  is  1  atmospheric  and  magnetic  field  routines  are  initialized 
at  the  location  of  the  path. 

Propagation  and  noise  quantities  are  calculated  and  stored  for  up 
to  41  links.  Link  indexes  LLINK  and  LNINK  are  used  to  store  and  retrieve 
propagation  and  noise  data.  If  LFLAG=1,  propagation  quantities  for  the  link 
have  been  previously  calculated  and  the  call  to  routine  PROPEN  is  skipped. 

If  NFLAG=1,  noise  calculations  have  been  made.  However,  if  NFLAG=0  further 
tests  (not  shown  in  Figure  22)  are  required  to  determine  if  noise  calculations 
are  required.  If  they  are,  routine  ANTND  is  called.  If  the  receiver  is  being 
jammed,  routine  JAMND  is  called  to  determine  the  jamming  noise  density.  Then, 
routine  SATLNK  is  called  to  define  link  performance  and  prepare  summary  output. 

After  calculations  for  all  the  links  in  a  particular  circuit  have 
been  performed,  routine  SATCIR  is  called  to  determine  circuit  performance 
quantities. 
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Figure  22.  Flowchart  for  subroutine  SATC. 
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SUBROUTINE  SATCIR 

This  routine  determines  performance  quantities  for  satellite  com¬ 
munication  circuits.  A  simplified  flow  diagram  for  the  routine  is  shown  in 
Figure  23. 

Circuit  performance  for  free  space  conditions  is  calculated  at  the 
first  calculation  time.  At  subsequent  calculation  times  circuit  performance 
is  only  calculated  for  disturbed  conditions.  For  either  free  space  (KAMBF=1) 
or  disturbed  (KAMBF=0)  conditions  a  link  loop  is  established  and  the  signal-to- 
noise  ratio  calculated.  If  there  is  no  signal  processing,  the  ratio  is  an 
effective  ratio  for  the  links  in  the  circuit.  If  there  is  signal  processing, 
the  signal-to-r»oise  ratio  and  bit  energy-to-noise  ratio  are  computed  for  each 
link.  For  FM  systems  an  FM  improvement  factor  is  computed  and  for  disturbed 
conditions  the  intermodulation  noise  due  to  dispersion  is  determined. 

Next,  the  bit  error  rate  (analytic  models)  or  message  error  rate 
(data  base  models)  are  calculated  for  digital  modulation  system  (see  Figure 
24).  If  there  is  signal  processing,  the  error  rates  are  computed  for  each 
link  and  then  combined  into  an  effective  error  rate  for  the  circuit.  After  the 
error  rates  are  computed  for  disturbed  conditions  summary  output  is  prepared 
and  written  on  file  LFNO. 


182 


MODELS  / 


Figure  24.  Detail  of  error  rate  calculation  in  Figure  23. 


SUBROUTINE  SATLNK 


This  routine  determines  link  propagation  quantities  for  satellite 
communication  circuits. 

First,  the  transmitter,  receiver,  and  antenna  parameters  are  ob¬ 
tained  for  the  link  and  quantities  required  to  compute  received  signal  and  noise 
are  determined.  Routines  GAIN  and  SLOBE  are  used  to  determine  the  antenna  main 
beam  and  side  lobe  gains. 

Next,  the  received  signal  and  noise  for  free  space  and  disturbed 
conditions  are  determined  and  stored  in  common  block  SIGCOM.  Scintillation 
quantities  for  the  link  are  also  defined  and  stored.  If  there  is  beam  spread 
loss  due  to  scintillation,  the  received  signal  for  disturbed  conditions  is 
reduced. 

Output  describing  propagation  quantities  for  the  link  is 
prepared  and  written  on  file  LFNO  if  level  1  summary  output  is  requested. 
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ROUTINE  SPECDP 


This  routine  determines  effects  of  prompt  radiation  on  neutral  species 
concentrations  at  a  specified  altitude  below  100  km  at  a  specified  calculation 
time.  The  routine  is  essentially  the  same  as  routine  SPECDP  developed  for 
ROSCOE. 

First  mean  values  for  0^  and  a  are  computed.  Then  the  neutral 
species  concentrations  are  computed. 

Next,  the  species  concentrations  are  limited  by  the  ambient  neutral 
species  values.  For  nighttime  conditions  this  only  involves  limiting  the 
values  of  0,  NO,  and  O^(^A).  For  da>time  a  new  ambient  model  consistent  with 
the  amount  of  odd  nitrogen  produced  is  first  determined.  This  requires  an 
iterative  procedure  which  has  been  chosen  to  minimize  the  number  of  iterations. 
Then  the  value  of  odd  oxygen  (0+0^)  is  allowed  to  decay  from  the  burst  produced 
value  toward  the  new,  modified  ambient  value.  The  value  of  odd  oxygen  at  the 
time  of  interest  is  used  to  determine  the  other  species  by  assuming  photo¬ 
chemical  equilibrium.  An  iterative  procedure  is  required  to  determine  the 
photoequilibrium  distribution  of  species  consistent  with  the  values  of  odd 
oxygen  and  odd  nitrogen.  The  decay  or  build  up  of  species  is  not  allowed  to 
surpass  the  mod.ficd  ambient  values. 
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SUBROUTINE  SLINK 


This  routine  determines  whether  propagation  or  noise  calculations 
have  been  previously  made  for  a  link. 

The  matrix  INDEX  stores  indexes  used  to  identify  propagation 
quantities.  Propagation  quantities  for  up  to  41  links  can  be  stored.  If 
propagation  quantities  have  not  been  calculated,  the  link  index  is  incremented 
and  a  propagation  calculation  flag  (LFLAG)  is  set  to  require  propagation 
calculations.  If  propagation  calculations  for  the  link  have  been  made,  the 
link  index  is  stored  in  INDEX.  Tests  are  also  made  to  see  if  noise  calcula¬ 
tions  have  beer,  made  and  a  noise  calculation  flag  is  set.  Additional  noise 
tests  and  a  noise  index  are  calculated  in  routine  SATC. 

If  L=1  (first  circuit),  propagation  and  noise  calculations  are 
required.  For  L>1  propagation  calculations  are  not  required  for  JLINK=2 
(satellite-to-satellite  link),  since  there  is  only  one  such  link.  For  all 
other  links,  tests  are  made  to  see  if  propagation  and  noise  calculations  have 

been  made.  First,  tests  arc  made  to  see  if  a  previous  link  has  the  same 
terminal  locations.  If  it  does  and  the  propagation  frequencies  are  the  same, 
propagation  calculations  are  the  same  for  both  links.  If  the  propagation 
frequencies  arc  the  same  and  the  receiver  antennas  are  the  same,  noise 
calculations  are  the  same  for  both  links  unless  one  link  is  an  up  link  and 
the  other  is  a  down  link. 
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SUBROUTINE  START 

This  routine  prepares  an  array  of  mandatory  points,  ordered  in 
increasing  slant  range  from  the  transmitter,  which  must  be  considered  when 
integrating  propagation  quantities  along  a  path. 

Mandatory  points  are  placed  at  the  minimum  and  maximum  slant  range 
as  computed  in  routine  PROPEN.  Mandatory  points  are  placed  at  the  entrance  and 
exit  points  of  the  ray  path  with  beta  tubes  and  fireballs  and  near  low-altitude 
gamma  sources.  In  preparing  the  mandatory  points  for  fireballs,  all  fireballs 
from  high-altitude  bursts  except  the  most  recent  burst  are  skipped.  The  mandatory 
points  are  contained  in  the  mandatory  point  array  SMP(I).  For  each  mandatory 
point  SMP(I)  there  are  two  integer  tags  in  the  array  INDMP  (2,100).  The  code 
for  INDMP  is 

INDMP(1,I)  =  0  Not  used. 

=  J  Fireball  index  associated  with 

the  intersection  points 

INDMP(2,I)  =  1000  Beta  tube,  low-altitude  gamma  points 

=  2000  Fireball  entry  point 

*  2002  Fireball  exit  point 
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SUBROUTINE  WOG1 

This  routine  determines  energy  deposition  coefficient  for  prompt  and 
delayed  gamma  rays.  The  routine  is  essentially  the  same  as  routine  W0G1 
developed  for  ROSCOE. 

The  effective  normalized  energy  deposition  coefficient  for  the 
spectrum  UBARGP(I,K)  is  calculated  and  stored  for  each  of  15  different  mass 
depths,  PMG(K).  Note  that  index  I  specifies  device  type,  index  J  specifies 
spectral  energy  group,  and  index  K  specifies  mass  depth.  The  double-dimensional 
array  UG(K,J)  stores  energy  deposition  coefficient  data  from  the  ATR  code  for 
monochromatic  energy  sources  at  photon  energy  corresponding  to  spectral  group 
J  and  at  mass  depths  PMG(K).  The  calculation  of  UBARGP  is  a  convolution  of  this 
monochromatic  data  with  the  energy  spectrum.  The  effective  normalized  energy 
deposition  coefficient  for  the  delayed  spectrum,  UBARGD(I,K),  is  calculated  and 
stored  for  each  of  the  15  mass  depths  in  the  same  manner  as  for  the  prompt 
energy  deposition  coefficient. 

The  parameters  FEDG(I)  and  FTHERM(I)  for  device  type  I  are  made 
available  to  the  other  phenomenology  models  through  labeled  common  block  DEVICE. 
The  tabular  data  UBARGP(I,K),  UBARGD(I,K),  and  PMG(K)  are  transmitted  to 
subroutines  WOGP  and  WOGD  through  labeled  common  block  WOG. 
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SUBROUTINE  WON1 

This  routine  determines  energy  deposition  coefficients  for  neutron 
energy  deposition.  The  routine  is  essentially  the  same  as  routine  W0N1 
developed  for  ROSCOE. 

First  the  energy  deposited  by  each  spectral  group,  the  total  energy 
deposited,  and  a  mean  neutron  cross  section  are  computed.  Then  the  yield 
fraction  in  deposited  neutron  energy  is  calculated  as  the  total  neutron  energy 
deposited  divided  by  the  total  weapon  yield.  Note  that  this  differs  from  the 
kinetic  energy  yield  fraction  given  as  input  in  that  it  includes  both  the 
kinetic  and  capture  energies. 

The  folded  energy  deposition  coefficients  for  prompt,  UBARNP, 
clastic  scatter,  UBARNE,  and  capture,  UBARNC,  phases  of  the  energy  deposition 
are  computed  for  each  mass  depth  PMN(K).  This  would  form  a  complete  set  of 
data  to  determine  the  prompt  and  delayed  energy  deposition.  However,  since 
the  computation  of  the  power  law  exponent  for  the  elastic  scatter  deposition 
phase  requires  a  time-consuming  iterative  solution,  we  make  that  calculation 
here  for  each  mass  depth  PMN(K) .  The  value  of  the  exponent  for  any  mass  depth 
is  found  later  by  interpolation.  The  parameters  UBAkWP,  UBARNC,  and  the 
exponent  F.XPQNE  then  form  a  complete  set  and  are  passed  to  subroutines  WONP 
and  WOND  through  WON  common,  along  with  the  mass  penetrated  array  PMN.  The 
yield  fraction  FEDN(I)  is  made  available  to  the  other  phenomenology  models 
through  DEVICE  common. 
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SUBROUTINE  WOX1 


This  routine  determines  energy  deposition  coefficients  for  X-rays. 
The  routine  is  essentially  the  same  as  routine  WOX1  developed  for  ROSCOE. 

In  addition  to  the  effective  energy  deposition  coefficient,  iT^  , 
an  X-ray  containment  factor,  f  ,  is  also  computed  for  use  in  the  fireball 
phenomenology  model  as 


f 

cx 


(M) 


y  (m)  dm 
x 


0) 


where  M  is  the  air  mass  penetrated  (gm/cm3).  The  function  fcx(M)  is  the 
fraction  of  the  X-ray  flux  that  would  be  absorbed  within  a  mass  depth  M. 

Because  the  X-ray  transport  may  vary  significantly  from  one  spectrum 
to  another,  the  deposition  coefficient,  URARX,  and  the  containment  factor, 
FCONTX,  are  precomputed  and  stored  for  a  different  set  of  air  masses  for  each 
device  type.  The  set  of  air  masses  for  device  type  I  are  given  by 
PMXBAR(I) *PMXR(K) ,  K  =  1,15.  The  normalizing  mass  PMXBAR(I)  is  defined  by 


18 

PMXBAR(I)  =  £  Oy.M  .  gm  cm'3  (2) 

j  =  l 


where  (F  ) .  is  the  energy  in  spectral  group  j 
^  J 

is  the  mass  normalization  for  spectral  group  j. 
defined  as  follows: 


for  device  type  I  and 
The  15  values  of  PMXR(K)  are 


K 

PMXR 

K 

PMXR 

K 

PMXR 

1 

0.001 

6 

0.4 

11 

2.0 

2 

0.01 

7 

0.6 

12 

3.0 

3 

0.05 

8 

0.8 

13 

4.0 

4 

0.1 

9 

1.0 

14 

5.0 

5 

0.2 

10 

1.5 

15 

10.0 

191 


I 


The  integral  of  the  X-ray  spectrum,  ETOTAL,  and  the  normalizing  inass 
PMXBAR(I)  are  computed  for  each  normalized  mass  depth.  Then  X-ray  energy 
deposition  coefficient,  UBARX,  and  the  containment  factor  FCONTX,  are  computed. 
(An  exponential  algorithm  is  employed  to  perform  the  integral  in  Equation  1.) 
The  normalizing  mass  PMXBAR,  and  the  tabulated  values  for  PMXR,  UBARX,  and 
FCONTX  are  placed  in  WOX  common  for  later  interpolation  by  routines  WOXP  and 
WOXC.  The  yield  fraction,  FEDX,  is  made  available  to  the  other  phenomenology 
routines  through  labeled  common  DEVICE. 
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ATTN:  TN 

Air  Force  Weapons  Laboratory 
Air  Force  Systems  Command 
ATTN:  NTYC 
ATTN:  SUL 
ATTN:  NTN 

Air  Force  Wright  Aeronautical  Lab 
ATTN  A.  Johnson 
ATTN:  W.  Hunt 

Air  Logistics  Command 
Department  of  the  Air  Force 

ATTN:  00-ALC/MM.  R.  Blackburn 

Air  University  Library 
Department  of  the  Air  Force 
ATTN:  AUL-LSE 

Air  Weather  Service.  MAC 
Department  of  the  Air  Force 
ATTN-  DNXP,  R  Babcock 

Assistant  Chief  of  Staff 
Intel  1 igence 

Department  of  the  Air  Force 
ATTN:  INED 

Assistant  Chief  of  Staff 
Studies  A  Analyses 
Department  of  the  Air  Force 

ATTN:  AF/SASC.  C.  Rfghtmejcr 
ATTN:  AF/SASC.  W.  Keaus 
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Ballistic  Missi le  Office 
Air  Force  Systems  Command 

ATTN  MNNXH,  J.  Allen 
ATTN:  MNNL ,  S  Kennedy 

Deputy  Chief  of  Staff 
Operations  Plans  and  Readiness 
Department  of  the  Air  Force 
ATTN  AFXOKT 
ATTN-  AFXOXFD 
ATTN  AFXOKCO 
ATTN:  AFXOKS 

Deputy  Chief  of  Staff 
Research,  Development,  &  Acq 
Department  of  the  Air  Force 
ATTN:  AFRDSP 
ATTN:  AFRDS 
ATTN:  AFRDSS 

Electronic  Systems  Division 
ATTN:  DCKC.  J.  Clark 

Electronic  Systems  Division 
ATTN:  XRW,  J.  Deas 

Electronic  Systems  Division 
ATTN:  YSM,  J.  Kobelski 
ATTN:  YSEA 

Foreign  Technology  Division 
Air  Force  Systems  Command 
ATTN:  NIIS  library 
ATTN:  TQTD.  B.  Ballard 

Headquarters  Space  Division 
Air  Force  Systems  Contnand 
ATTN:  SKA,  D.  Bolin 
ATTN:  SKY,  C.  Kennedy 

Headquarters  Space  Division 
Air  Force  Systems  Comnand 
ATTN:  YZJ ,  W.  Mercer 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN:  E.  Butt 

Rome  Air  Development  Center 
Air  Force  Systems  Command 
ATTN  TSID 
ATTN:  OCS.  V.  Coyne 

Rome  Air  Development  Center 
Air  Force  Systems  Command 
ATTN:  EEP 

Strategic  Air  Conmand 
Department  of  the  Air  Force 
ATTN:  ADWATE,  B.  Bauer 
ATTN:  OCX 
ATTN:  DCXT 
ATTN  NRT 

ATTN  DCXR,  T.  Jorgensen 
ATTN:  XPFS 

ATTN:  XPFS,  B.  Stephan 


OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSUR/NED 

Department  of  Commerce 

National  Burea  i  of  Standards 

ATTN:  Sec  Ofc  for  R  Moore 

Department  of  Commerce 

National  Oceanic  &  Atmospheric  Admin 
ATTN-  R.  Grubb 

Institute  for  Telecommunications  Sciences 
ATTN-  U.  Utlaut 
ATTN:  A.  Jean 
ATTN:  L  Berry 

U.S.  Coast  Guard 

Department  of  Transportation 

ATTN:  G-DOE-3/TP54.  B  Romine 

DEPARTMENT  OF  ENERCY  CONTRACTORS 

EGAG,  Inc 

ATTN:  D.  Wright 
ATTN-  J.  Colvin 

Lawrence  Livermore  National  Lab 

ATTN-  Tech  Info  Dept  Library 
ATTN:  L-389,  R.  Ott 
ATTN:  L-31 .  R.  Hager 

Los  Alamos  National  Scientific  Lab 
ATTN  MS  664.  J  Zinn 
ATTN:  J.  Wolcott 
ATTN-  E.  Jones 
ATTN:  D.  Westervelt 
ATTN-  R.  Taschek 
ATTN:  D.  Simons 
ATTN-  R.  Jeffries 
ATTN-  P.  Keaton 
ATTN-  MS  670.  J.  Hopkins 

Sandia  laboratories 

Livermore  Laboratory 


ATTN 

B. 

Murphey 

ATTN 

T. 

Cook 

Sandia  National 

Lab 

ATTN 

D. 

Thornbrough 

ATTN 

3141 

ATTN 

D. 

Dahlgren 

ATTN 

Space  Project  D:v 

ATTN 

Org  1250.  W  Brown 

ATTN 

Orq  4Z41,  T.  Wright 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN 

J, 

Straus 

ATTN 

S. 

Bower 

ATTN 

r. 

Salmi 

ATTN 

D. 

01  sen 

ATTN 

V. 

Josephson 

ATTN 

R. 

Slaughter 

ATTN 

I. 

Garfunkol 

ATTN 

N. 

Stockwell 
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University  of  Alaska 

ATTN:  Technical  Library 
AT'TI:  N.  Brown 
ATTN:  T.  Davis 

Analytical  Systems  Engineering  Corp 
ATTN-  Radio  Sciences 

Analytical  Systems  Engineering  Corp 
ATTN:  Security 

Barry  Research  Corporation 
ATTN-  J.  McLaughlin 

BOM  Corp 

ATTN:  T.  Neighbors 
ATTN:  L.  Jacobs 

Berkeley  Research  Associates.  Inc 


ATTN 

C.  Prettle 

ATTN: 

J.  Workman 

BETAC 

ATTN: 

J.  Hirsch 

Boeing  Co 

ATTN: 

G.  Hall 

ATTN: 

S.  Tashlrd 

ATTN: 

M/S  42-33,  J.  Kennedy 

Boelnq  Co 

ATTN: 

0.  Clauson 

Boot- Allen  A 

Hamilton.  Inc 

ATTN: 

B.  Wilkinson 

University  of  California  at  San  Diego 

ATTN- 

H.  Booker 

Charles  Stark  Draper  Lab,  Inc 

ATTN. 

J.  Gilmore 

ATTN: 

D.  Cox 

Communications  Satellite  Corp 

ATTN: 

D.  Fang 

Comsat  Labs 

ATTN: 

R.  Taur 

ATTN: 

G.  Hyde 

Cornell  University 

ATTN: 

D.  Farley,  Jr 

ATTN: 

M.  Keily 

Electrospace  Systems,  Inc 

ATTN: 

H.  Logston 

ATTN: 

P.  Phillips 

ESL.  Inc 

ATTN: 

J.  Marshall 

General  Electric  Co 

ATTN: 

A.  Harcar 

ATTN: 

M.  Bortner 

General  Electric  Co 

ATTN: 

A.  Steinmayer 

ATTN: 

C.  Zierdt 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  yContinuedj 

General  Electric  Co 

ATTN:  F.  Relbert 

General  Electric  Tech  Services  Co.  Inc 
ATTN:  G.  Millman 

General  Research  Corp 

ATTN:  J.  Ise,  Jr 
ATTN-  J.  Garbarino 

Hcrirons  Technology,  Inc 
ATTN:  R.  Kruger 

HSS.  Inc 

ATTN:  D.  Hansen 

IBM  Corp 

ATTN:  F.  Ricci 

University  of  Illinois 

ATTN:  Security  Supervisor  for  K.  Yen 

Institute  for  Defense  Analyses 


ATTN 

H  Wolfhard 

ATTN 

E.  Bauer 

ATTN 

J.  Bengston 

ATTN 

J.  Aein 

International  Tel  A  Telegraph 

ATTN 

Technical  Library 

ATTN 

G  Wetmore 

JAYCOR 

ATTN 

J.  Sperling 

JAYCOR 

ATTN 

J.  Doncarlos 

Johns  Hopkins  University 

ATTN 

J  Phillips 

ATTN 

P  Komiske 

AT  IN 

T  Potemra 

ATTN 

T.  Evans 

ATTN 

J.  Newland 

Kamun  TEMPO 

ATTN 

T.  Stephens 

ATTN 

M.  Stanton 

ATTN 

DAS  I AC 

ATTN 

W.  Knapp 

ATTN 

W.  McNamara 

Lfnkablt  Corp 

ATTN 

I .  Jacobs 

Litton  Systems.  Inc 

ATTN 

R.  Grasty 

Lockheed  Missiles  A  Space  Co, 

ATTN 

M.  Walt 

ATTN 

R  Johnson 

ATTN 

W.  Imhof 

Lockheed  Missiles  A  Space  Co. 

ATTN 

D.  Churchill 

ATTN 

Dept  60-12 
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M  I  T  Lincoln  Lab 

Rand  Corp 

ATTN:  D.  Towle 

ATTN 

E  Bedrozian 

ATTN:  L.  loughlin 

ATTN 

C.  Crain 

Martin  Marietta  Corp 

Riverside  Research  Institute 

ATTN.  R.  Heffner 

ATTN 

V  Trapani 

McDonn»U  Douglas  Corp 

Rockwell  International  Corp 

ATTN:  Tech  library  Services 

ATTN  U.  Olson 

AT1N 

R.  Buckner 

ATTN:  G.  Mroz 

Rockwell  International  Corp 

ATTN.  R.  Halprin 

ATTN:  N.  Harris 

ATTN 

S  Quillci 

ATTN:  J.  Moule 

Santa  Fe  Corp 

Meteor  Communications  Consultants 

ATTN 

D  Paolucci 

ATTN:  R.  Leader 

Science  Applications.  Inc 

ATTN 

C.  Smith 

Mission  Research  Corp 

ATTN 

D  Hamlin 

ATTN  R.  Bogusch 

ATTN 

L.  Llnson 

ATTN:  R.  Hendrick 

ATTN-  R.  Ki lb 

ATTN 

E.  Straker 

ATTN  S  Gutsche 

Science  Appl ications.  Inc 

ATTN:  F.  Fajen 

ATTN:  Tech  Library 

ATTN 

S  l 

ATTN  D.  Sappenf ield 

Science  Applications,  Inc 

Mitre  Corp 

ATTN 

J.  Cockayne 

ATTN:  A  Kymmel 

SRI  International 

ATTN:  C.  uallhan 

ATTN 

W  Jaye 

ATTN:  B.  Adams 

ATTN 

W.  Chesnut 

ATTN:  G.  Harding 

ATTN 

R.  Leadab’-and 

ATTN 

R  Hake,  Jr 

Mitre  Corp 

ATTN 

R  Tsunoda 

ATTN:  M.  Horrocks 

ATTN 

A  Bums 

AIN:  J.  Wheeler 

ATTN 

G.  Price 

A1FN:  U  Foster 

ATTN 

G.  Smith 

ATTN:  W.  Hall 

ATTN 

V.  Gonzales 

ATTN 

D.  Nell  son 

Pacific-Sierra  Research  Corp 

ATTN 

M  Baron 

ATTN:  E  Field.  Jr 

ATTN 

J  Petrlckes 

ATTN:  H.  Brode 

ATTN 

R.  Livingston 

ATTN:  f.  Thomas 

ATTN 

C  Rino 

Pennsylvania  State  University 

ATTN 

D.  McDaniels 

ATTN:  Ionospheric  Research  Lab 

Syi vanla  Systems  Group 

Photometries,  Inc 

ATTN 

M.  Cross 

ATTN:  I.  Kofsky 

Technology 

International  Corp 

Physical  Dynamics,  Inc 

ATTN 

W.  Boquist 

ATTN:  E.  Fremouw 

Trl-Com,  Inc 

R  &  D  Associates 

ATTN 

D  Murray 

ATTN:  R.  Lelevier 

TRW  Defense  &  Space  Sys  Group 

ATTN:  W.  Wright 

ATTN 

R  Plebuch 

ATTN  W.  Karzas 

ATTN  F.  Gilmore 

ATTN 

D.  Dee 

ATTN:  G.  Gabbard 

Utah  State 

University 

ATTN-  C  Greifinger 

ATTN 

J.  Dupnik 

ATTN :  H  Ory 

ATTN 

l.  Jensen 

ATTN.  R.  Turco 

ATTN-  M  Gantsweg 

ATTN 

K.  Baker 

ATTN  P.  Haas 

VisIDyne.  Inc 

ATTN 

C.  Humphrey 

PAD  Associates 

ATTN:  L.  Delaney 

ATTN:  B.  Yoon 

ATTN 

J.  Carpenter 
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